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Kurzfassung 


Nachlaufströmungen existieren hinter fluidumströmten Körpern und beeinflus- 
sen die auf den Körper wirkenden aerodynamischen Kräfte. Für Überschall- 
flugkörper kann sich eine Veränderung der Nachlaufströmung beispielsweise 
maßgeblich auf den Strömungswiderstand auswirken. Da die Auswirkungen 
von strömungsmechanischen Einflüssen auf Nachlaufströmungen bis heute 
nicht vollständig verstanden werden, ist die Untersuchung von Nachlaufströ- 
mungen nicht nur für technische Anwendungen im Bereich der Fahrzeug- und 
Flugkörperaerodynamik relevant, sondern auch aus wissenschaftlicher Sicht 
von grundlegender Bedeutung. 


Im Rahmen der vorliegenden Arbeit wurde der Einfluss eines vor der Modell- 
hinterkante in die Strömung eingebrachten Dralls auf den Überschallnachlauf 
eines längsangeströmten zylindrischen Körpers untersucht. 


In der Literatur ist mittels numerischer Untersuchungen bereits eine Verände- 
rung der Nachlaufstruktur aufgrund eines eingebrachten Dralls gezeigt worden. 
Experimentell konnte diese Veränderung der Nachlaufstruktur jedoch bisher 
nicht nachgewiesen werden. 


Der erste Teil der vorliegenden Arbeit umfasst experimentelle Untersuchungen 
zu rotierenden Überschallnachlaufströmungen. Zunächst wurden experimen- 
telle Untersuchungen mit rotierenden Strömungskörpern durchgeführt, die ge- 
genüber nichtrotierenden Nachlaufströmungen einen Druckabfall im Zentrum 
der Modellbasis erzeugten. Eine strukturelle Veränderung des Nachlaufs konn- 
te jedoch nicht gezeigt werden. Deshalb wurden anschließend Windkanalver- 
suche durchgeführt, bei denen festmontierte Strömungskörper mit angestellten 
Leitflächen höhere Drallraten im Nachlauf erzeugten. Durch die Vermessung 
des Geschwindigkeitsfeldes hinter den unterschiedlichen Strömungskörpern 
sowie Druckmessungen an der Basis der Strömungskörper konnte experimen- 
tell nachgewiesen werden, dass es aufgrund des eingebrachten Dralls zu zwei 
grundsätzlich verschiedenen, von der Drallrate abhängigen Veränderungen der 
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Kurzfassung 


Nachlaufstruktur kommt. Für einen Anstellwinkel der Leitflächen von 2 = 16° 
bildete sich eine rotierende Nachlaufströmung mit einem zentralen, stromab- 
gerichteten Wirbelschlauch aus, wie sie bisher nur durch numerische Simu- 
lationen vorhergesagt worden war. Ein höherer Leitflächenanstellwinkel von 
A = 32° führte zu einer bisher nicht beschriebenen Strukturänderung der 
Nachlaufströmung. In der Nähe der Strömungsachse bildete sich eine strom- 
aufgerichtete Strömung aus, die vom fernen Nachlauf bis zur Basis des Strö- 
mungskörpers reichte. 


Um die strömungsmechanischen Zusammenhänge sowie die Ursachen für die 
experimentell beobachteten Veränderungen der Strömungsstruktur beschrei- 
ben zu können, wurden im zweiten Teil dieser Arbeit numerische Simulationen 
durchgeführt. Mit Hilfe turbulenzauflösender Grobstruktursimulationen konn- 
ten die experimentellen Strömungsfelder qualitativ gut modelliert werden. Die 
Simulationsergebnisse für A = 16° zeigen, dass aufgrund des in den Nach- 
lauf eingebrachten Dralls in einem entlang der Bahnlinie eines Fluidelemen- 
tes mitgeführten und mitrotierenden Bezugssystems Zentrifugalkräfte erzeugt 
wurden. Die Zentrifugalkräfte resultierten in einem gegenläufig rotierenden, 
konzentrischen und torusförmigen Wirbelpaar an der Basis des Strömungskör- 
pers und führten zur Ausbildung des experimentell beobachteten, zentralen, 
stromabgerichteten Wirbelschlauches. Für A = 32° verhinderten die Zentrifu- 
galkräfte eine Ablenkung der an der Modellhinterkante ablösenden Strömung 
in Richtung der Zentralachse, sodass sich ein zylinderförmiges Totwasserge- 
biet ausbildete. Begrenzt wurde das Totwassergebiet im fernen Nachlauf durch 
die an den Windkanalwänden reflektierten Verdichtungsstöße, die an den Leit- 
flächenvorderkanten entstanden. Ein Teil der Strömung konnte den durch die 
reflektierten Verdichtungsstöße erzeugten axialen Druckgradienten nicht über- 
winden und wurde in Richtung der Basis umgelenkt. 


Die experimentellen und numerischen Ergebnisse zeigen außerdem, dass die 
Veränderungen der Nachlaufstruktur mit einer Verringerung des Basisdru- 
ckes um bis zu 80 % gegenüber einem nichtrotierenden Nachlauf einhergehen. 
Der eingebrachte Drall resultierte somit zusätzlich zu den strukturellen Ver- 
änderungen in einer maßgeblichen Erhöhung des Strömungswiderstandes in 
Abhängigkeit von der Drallrate der Nachlaufströmung. 
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Abstract 


Wake flows exist behind bodies in a fluid flow and affect the aerodynamic 
forces acting on these bodies. Small changes in the wake flow properties can 
for example significantly alter the total drag of supersonic flight vehicles. 
There are, however, effects on the wake that are not yet fully understood. 
Wake flows are therefore not only relevant for technical applications as vehicle 
aerodynamics, but also of a fundamental interest in scientific research. 


The present work investigates the effect of swirl introduced upstream of the 
base corner on the wake of axisymmetric afterbodies in a supersonic flow. 


Previous numerical studies in the literature indicated a change in the wake flow 
structure caused by spinning finned afterbodies, however, no experimental 
proof was found. 


Therefore, wind tunnel experiments on swirling supersonic wake flows were 
conducted in the first part of the present work. The first study with physically 
spinning afterbodies shows a decrease of the central base pressure compared 
to non-swirling wake flows. A structural change of the wake, however, was not 
observed. Hence, wind-tunnel experiments on models having stationary but 
canted fins introducing a higher swirling motion to the wake were conducted. 
Dependent on the swirl rate introduced by the canted fins, measurements of 
the velocity field behind the models and the pressure distribution at the model 
base show evidence of two distinctively different wake flow structures. For a fin 
cant angle of A = 16°, a rotating wake flow with a central downstream vortex 
tube formed, which was previously described only by numerical simulations. 
A higher fin cant angle of 2 = 32° resulted in a change of the wake structure, 
which was observed for the first time. Near the central axis of the flow field an 
upstream flow was established, extending from the far wake up to the model 
base. 


In order to describe the fluid-dynamical processes and the origins for the 
experimentally observed structural changes of the rotating wakes, numerical 


Abstract 


simulations were performed in the second part of this work. The results of large- 
scale turbulence-resolving simulations are qualitatively in good agreement with 
the measured flow fields. The numerical results for the model with 2 = 16° 
show that the swirling flow motion in the wake introduced centrifugal forces 
in the non-inertial Lagrangian reference frame following the pathline of an 
individual fluid parcel. The centrifugal forces led to a concentric, counter- 
rotating, toric vortex pair near the model base and resulted in the formation 
of the experimentally-observed central downstream vortex tube. For A = 32°, 
the centrifugal forces acting on the swirling fluid prevented the flow detaching 
at the base corner from being deflected toward the central axis. This resulted 
in a cylindrical dead air region behind the model. The dead air region was 
bounded in the far wake by the shocks caused by the fin leading edges, which 
were reflected from the wind tunnel walls. Parts of the flow were not capable 
of overcoming the adverse pressure gradient caused by the reflected shocks 
and thus were directed back to the base. 


Additionally, the experimental and numerical results show that the changes 
in the wake flow structure go along with a decrease in base pressure of up 
to 80 % in comparison to a non-rotating wake. The swirling motion resulted 
therefore not only in two distinct changes of the wake flow structure but also 
in a significant increase in drag dependent on the swirl rates in the wake. 
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417 äußerer Wirbelring im Nachlauf 
2 innerer Wirbelring im Nachlauf 
© zentraler, stromabgerichteter Wirbelschlauch im Nachlauf 


Lateinische Indizes 


b an der Modellbasis 

be Randbedingung unmittelbar hinter der Modellhinterkante 

B im Ballon 

e unmittelbar vor der Modellhinterkante 

i Zählvariable 

J Versuchsnummer 
bezogen auf das Rechengitter mit dem Index m 

min Uf bezogen auf die maximale relative 
Rückströmgeschwindigkeit Uf 

n normal zum Verdichtungsstoß 

NLTV aufgrund von Nichtlinearitäten, Hysteresen und 
Zeitabhängigkeiten 

Poo bezogen auf den statischen Druck im Freistrom pæ 

Po bezogen auf den Basisdruck pp 

Pe bezogen auf den relativen Basisdruck p% 

Pe bezogen auf den mittleren relativen Basisdruck p* 

Dr bezogen auf die Summe des wahren Wertes der Messgröße 
und der systematischen Abweichungen p% 

PIV zeitlicher Mittelwert der PIV-Messungen 

r Komponente in radialer Richtung r 

ref. Referenzmessung für PSP-Methode 
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rot. mit Modellrotation 
Sim. zeitlicher und azimutaler Mittelwert der 
numerischen Simulation 
TO aufgrund einer temperaturbedingten Nullpunktverschiebung 
Tl aufgrund einer temperaturbedingten Sensitivitätsveränderung 
VB im Vakuumbehälter 
x Komponente in axialer Richtung x 
X|minU* bezogen auf die axiale Position der maximalen relativen 
Rückströmgeschwindigkeit min U* 
XSP bezogen auf die axiale Position des hinteren 


Stagnationspunktes xsp 


Griechische Indizes 
y Komponente in azimutaler Richtung y 


Numerische Indizes 


0 Ruhebedingung, Totalgröße 
1 vor dem Verdichtungsstoß 

2 hinter dem Verdichtungsstoß 
o0 im Freistrom 


Bezeichnungsgrundsätze 


X Summe aus dem wahren Wert der Messgröße (X),.; und den 
systematischen Messabweichungen >); Bx i 

(X) real wahrer Wert der Messgröße X 

X” normiert mit Freistrombedingung X% 

X* komplex Konjugierte von X 

X Mittelwert von X 

X’ zeitlicher Schwankungswert der Größe X 

ox partielle Ableitung von X in radialer Richtung r 

Konstanten 


Ruut = 287,1 J/(kgK) spezifische Gaskonstante für trockene Luft 
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CCD Charge-Coupled-Device-Sensor 

CFL Courant-Friedrichs-Lewy 

CI Konfidenzintervall 

DES Detached-Eddy Simulation 

DNS Direkte Numerische Simulation 

GCI Grid Convergence Index, Gitterkonvergenzindex 

k-w-SST Wirbelviskositätsmodell nach Menter (1994) 

LDA Laser-Doppler-Anemometrie 

LES Large-Eddy Simulation 

LF Leitfläche 

MILES Monotone-Integrated Large-Eddy Simulation 

Nd:YAG Neodym-dotierter Yttrium-Aluminium-Granat-Kristall 

PIV Particle Image Velocimetry 

PSP drucksensitive Farbe (pressure-sensitive paint) 

PtTFPP Platin(I)-5,10,15,20-Tetrakis-(2,3,4,5,6-Pentaflurophenyl)- 
Porphyrin 

RANS Reynolds-gemittelte Navier-Stokes-Gleichungen 

RMS quadratischer Mittelwert (root mean square) 

SA Wirbelviskositätsmodell nach Spalart und Allmaras (1992) 

SGS Sub-Grid Scale 

SP hinterer Stagnationspunkt 

SST Shear-Stress Transport 
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1 Einleitung 


1.1 Motivation 


Als Nachlauf wird das Strömungsgebiet bezeichnet, welches sich unmittelbar 
hinter dem von einer Strömung umschlossenen Körper ausbildet. Strömungen 
dieser Art sind sowohl aus grundlagenorientierter, wissenschaftlicher Sicht 
als auch bei der Entwicklung von Flugkörpern und Fahrzeugen relevant, da 
sie die aerodynamischen Kräfte des Strömungskörpers beeinflussen können. 
Ein hoher Einfluss ergibt sich zum Beispiel auf den Strömungswiderstand von 
Überschallflugkörpern mit einem flachen Heck (Lamb und Oberkampf, 1995). 
Ähnliche Strömungen werden auch zur grundlagenorientierten Untersuchung 
von turbulenten, kompressiblen Scherschichten (Demetriades, 1968, 1976) 
verwendet. 


In Abbildung 1.1 ist der klassische Überschallnachlauf hinter einem achsen- 
symmetrischen Strömungskörper skizziert, wie er bereits von Chapman (1950) 
gezeigt wurde. Die Strömung expandiert an der Hinterkante des Körpers. Der 
maximale Winkel der Strömungsumlenkung durch die Expansion ist begrenzt, 
sodass die Grenzschicht an der Hinterkante des Körpers ablöst. Aus der Grenz- 
schicht entwickelt sich zwischen der äußeren und der inneren Strömung des 
nahen Nachlaufgebietes eine in Richtung der zentralen Achse des Strömungs- 
feldes abgelenkte freie Scherschicht. Das Fluid wird in der Nähe der Achse 
erneut umgelenkt, sodass sich in der äußeren Überschallströmung ein Rekom- 
pressionsstoß ausbildet. Durch die Kompression entsteht ein ansteigender 
Druckgradient und ein Teil der inneren Strömung wird aufgrund des gerin- 
geren Totaldrucks gegenüber der äußeren Strömung in Richtung des Körpers 
umgelenkt, sodass ein Rezirkulationsgebiet entsteht (Korst, 1956). Der an 
der zentralen Achse entstehende hintere Stagnationspunkt (SP) unterteilt das 
Strömungsgebiet in den nahen und fernen Nachlauf. Die Isotache mit der 
Axialgeschwindigkeit U, = 0 grenzt die stromauf- und stromabgerichteten 
Strömungsbereiche voneinander ab (Lamb und Oberkampf, 1995). 
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Abbildung 1.1: Skizze einer zeitlich gemittelten, nichtrotierenden Uberschallnachlaufstromung 
nach Chapman (1950) hinter einem achsensymmetrischen Strömungskörper (Weidner et al., 2017, 
2019b). 


Der Zustand der Strömung im Rezirkulationsgebiet ist von vielen Faktoren 
abhängig. In der Literatur konnte bereits der Einfluss verschiedener Mach- 
und Reynoldszahlen (Sieling und Page, 1970, Cope, 1953), Grenzschicht- 
eigenschaften (Augenstein et al., 1999, Durgesh et al., 2013) und geometri- 
scher Veränderungen an der Modellhinterkante (Sahu et al., 1985) auf die 
Eigenschaften von Nachlaufströmungen nachgewiesen werden. 


Hruschka und Leopold (2015) untersuchten, welchen Einfluss die Rotation 
von Strömungskörpern auf deren Überschallnachlauf hat. Mittels numeri- 
scher Simulationen konnte eine Veränderung der Strömungsstruktur im 
Nachlaufgebiet und eine Abnahme des Basisdrucks gezeigt werden. Stabilitäts- 
theoretische Untersuchungen für eine laminare, inkompressible Nachlauf- 
strömung lassen ebenfalls eine durch Drall hervorgerufene Veränderung der 
Strömungsstruktur erwarten (Jiménez-Gonzalez et al., 2014). Allerdings sind 
experimentelle Nachweise in der Literatur bisher weder für die von Hruschka 
und Leopold (2015) simulierten noch für die von Jimenez-Gonzälez et al. 
(2014) theoretisch vorhergesagten Änderungen der Nachlaufstruktur zu finden. 
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1.2 Zielsetzung 


In dieser Arbeit soll der Einfluss des Dralls auf die Nachlaufströmung hinter 
einem achsensymmetrischen Körper mit stumpfer Basis untersucht werden. 
Insbesondere sollen die von Hruschka und Leopold (2015) mittels numerischer 
Simulationen gezeigte Änderung der Strömungsstruktur eines Überschall- 
nachlaufs experimentell überprüft und die zugrundeliegenden strömungs- 
mechanischen Prozesse erläutert werden. 


1.3 Vorgehensweise 


Zunächst wurden in einem Überschallwindkanal experimentelle Unter- 
suchungen zu rotierenden Nachlaufströmungen durchgeführt. In einer ersten 
Studie wurde der zentrale Basisdruck rotierender Strömungskörper unter- 
sucht (Hruschka und Leopold, 2015, Weidner et al., 2019b). Anschließend 
wurden festmontierte Strömungskörper verwendet, die mittels vor der Modell- 
hinterkante angebrachter, gegenüber der Anströmung angestellter Leitflächen 
einen höheren Drall in der Strömung erzeugten, als dies aufgrund technischer 
Limitierungen mit rotierenden, geraden Leitflächen möglich war (Hruschka 
und Leopold, 2015, Weidner et al., 2019b). Um die von Hruschka und 
Leopold (2015) mittels numerischer Simulationen gezeigte Veränderung der 
Strömungsstruktur aufgrund eines eingebrachten Dralls experimentell zu 
überprüfen, wurden Messungen des Basisdruckes sowie der Strömungsge- 
schwindigkeiten im Nachlauf der nichtrotierenden Strömungsmodelle mit 
unterschiedlichen Leitflächenanstellwinkeln vorgenommen. 


Zur Beschreibung der strömungsmechanischen Prozesse rotierender Nachlauf- 
strömungen sind zusätzlich zu den experimentell untersuchten Größen die 
Azimutalgeschwindigkeit der Strömung sowie deren Einfluss auf die Druck- 
verteilung des Nachlaufs von Interesse. Diese beide Größen konnten im 
Rahmen der vorliegenden experimentellen Untersuchungen jedoch nicht 
bestimmt werden, sodass diese mit Hilfe von numerischen Simulationen 
untersucht wurden. Die abgelöste Strömung hinter dem Modell stellt für 
numerische Simulationen allerdings eine Herausforderung dar. Daher wurde 
die Gültigkeit der Simulationsergebnisse zunächst durch einen Vergleich mit 
den experimentellen Ergebnissen sichergestellt. 
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1.4 Aufbau der Arbeit 


Kapitel 2 gibt einen Überblick über die wichtigsten vorhergehenden Arbeiten 
über Nachlaufströmungen im Überschall. Hier werden sowohl analytische 
Modelle als auch empirische Korrelationen beruhend auf experimentellen 
Untersuchungen beschrieben. Anschließend werden die wichtigsten numeri- 
schen Simulationsergebnisse aus der Literatur vorgestellt und die Eignung der 
unterschiedlichen numerischen Methoden zur Berechnung von Überschall- 
nachlaufströmungen diskutiert. 


In der Literatur sind nur wenige Untersuchungen zu rotierenden Nachlauf- 
strömungen zu finden. Deshalb werden in Kapitel 2 mögliche Effekte eines 
Dralls auf den Nachlauf anhand von rotierenden Jetströmungen beschrieben. 


Kapitel 3 beschreibt den Versuchsaufbau, die Strömungsbedingungen, die 
eingesetzten Modellgeometrien sowie die durchgeführten Messungen des 
zentralen Basisdrucks für die Studie mit den rotierenden Strömungskörpern. 
In Kapitel 3 wird auch eine stark vereinfachte Modellbeschreibung zur 
Abschätzung der Basisdruckänderungen aufgrund der Modellrotation vor- 
gestellt und mit den experimentellen Ergebnissen verglichen. Außerdem be- 
inhaltet Kapitel 3 den Vergleich zwischen rotierenden Strömungskörpern mit 
geraden Leitflächen und nichtrotierenden Strömungskörpern mit angestellten 
Leitflächen als alternative Möglichkeit zur Drallerzeugung innerhalb einer 
Nachlaufströmung, deren Äquivalenz für geringe Drallraten ebenfalls in 
Kapitel 3 gezeigt wird. 


Kapitel 4 beschreibt die Verwendung der nichtrotierenden Strömungskörper 
mit angestellten Leitflächen zur Erzeugung rotierender Nachlaufströmungen, 
deren Drallraten aufgrund technischer Limitierungen mit den zuvor unter- 
suchten rotierenden Strömungskörpern nicht möglich waren. Außerdem wird 
in Kapitel 4 die Umsetzung der verwendeten Messmethoden zur Untersuchung 
rotierender Überschallnachlaufströmungen anhand eines längsangeströmten 
Zylinders ohne eingebrachte Rotation beschrieben und die zu erwartenden 
Messunsicherheiten analysiert. Die Messergebnisse für den längsangeströmten 
Zylinder wurden zusätzlich mit den vergleichbaren Experimenten von Leopold 
(1993) validiert. 


Die experimentellen Ergebnisse für rotierende Überschallnachlaufströmungen, 
die mit den nichtrotierenden, angestellten Leitflächen erzeugt wurden, werden 
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in Kapitel 5 beschrieben und auf Strukturänderungen in der Nachlaufströmung 
untersucht. Im Speziellen wurde die von Hruschka und Leopold (2015) 
simulierte strukturelle Veränderung des Nachlaufes experimentell überprüft. 


Da die Messergebnisse aus Kapitel 5 nur einen kleinen Ausschnitt des 
Strömungsfeldes umfassen und somit keine vollständige Beschreibung aller 
relevanten strömungsmechanischen Prozesse zulassen, wurden numerische 
Simulationen des gesamten Versuchsaufbaus durchgeführt, wie sie in Kapitel 6 
beschrieben werden. 


Die Simulationsergebnisse wurden in Kapitel 7 zunächst mittels der Mess- 
ergebnisse aus Kapitel 5 validiert. Anschließend ermöglichten die validierten 
Simulationsergebnisse einen Einblick in sämtliche Strömungsgrößen des 
Strömungsfeldes. Das auf diese Weise zugängliche Gesamtbild der rotie- 
renden Nachlaufströmungen wurde in Kapitel 7 anschließend zur Erläuterung 
der strömungsmechanischen Prozesse verwendet, die zu den unterschiedlichen, 
von der Drallrate abhängigen Strömungsstrukturen führten. 


Kapitel 8 fasst die Ergebnisse dieser Arbeit zusammen und gibt einen Ausblick 
über zukünftige Untersuchungen, die sich im Anschluss an diese Arbeit 
ergeben. 
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2.1 Achsensymmetrischer Überschallnachlauf 


2.1.1 Empirische Korrelationen 


Nachdem Mach und Salcher (1887) die ersten Strömungsvisualisierungen von 
Überschallflugkörpern veröffentlicht hatten, wurden von Lorenz (1917, 1928), 
Gabeaud (1931) sowie von Kärmän und Moore (1932) empirische Ansätze 
entwickelt, um den Basisdruck berechnen zu können. Die eingeschränkte 
Datengrundlage sowie die Beschränkung auf einzelne Abhängigkeitsparameter 
begrenzen jedoch die Anwendbarkeit dieser Korrelationen. 


Mit der Verfügbarkeit einer Vielzahl an Druckmessdaten aus Freiflug- und 
Windkanalversuchen, wie sie beispielsweise von Hill und Alpher (1949) 
zusammengefasst worden waren, konnten detailliertere empirische Modelle 
entwickelt werden. Diese beschreiben beispielsweise den Einfluss der Grenz- 
schichtdicke auf den Basisdruck als Funktion der von der Länge des Strömungs- 
körpers abhängigen Reynoldszahl Re, (Cope, 1953). Andere Korrelationen für 
den Basisdruck, wie die von Hoerner (1950) und Kurzweg (1951), beinhalten 
zusätzlich zur Machzahl der Anströmung und der Reynoldszahl Re; die 
Einflüsse, die sich aus dem Verhältnis von Grenzschichtdicke zu Modelldurch- 
messer, der Oberflächenrauigkeit der Modelle und der Temperatur an der Basis 
der Strömungskörper ergeben (Murthy und Osborn, 1976). 


2.1.2 Theoretische Modelle 


Aus der Annahme, dass sich abgelöste Überschallströmungen in die in Ab- 
bildung 1.1 gekennzeichneten Strömungsbereiche unterteilen lassen, hatte 
Chapman (1950) geschlussfolgert, dass die viskosen Mischungsprozesse in 
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der Scherschicht essenziell für die Bestimmung des Basisdruckes sind. Unter 
der Annahme einer laminaren Grenzschicht mit vernachlässigbarer Dicke, 
einer Prandtl-Meyer-Expansion an der Modellhinterkante und einer isentropen 
Rekompression der Strömung wurde von Chapman (1950) ein Modell zur 
Beschreibung von Überschallnachlaufströmungen entwickelt. Eine ähnliche 
Modellbeschreibung, basierend auf der Annahme einer turbulenten Grenz- 
schicht mit einer ebenfalls vernachlässigbaren Dicke und einer nichtisentropen 
Rekompression über einen Verdichtungsstoß, wurde von Korst (1956) entwi- 
ckelt. Beide Modelle treffen die Annahme, dass der Totaldruck entlang der 
Stromlinie, die das Rezirkulationsgebiet von der äußeren Strömung abtrennt, 
dem Totaldruck stromab der Rekompression entspricht. Auf diese Weise 
lassen sich die Stromungsgeschwindigkeiten in der Scherschicht bestimmen, 
sodass sich der Basisdruck ausschließlich abhängig von der Machzahl vor der 
Strömungsablösung ergibt (Rom, 1965). Aufgrund der getroffenen Vereinfa- 
chungen sind die analytischen Lösungen nicht allgemeingültig und somit nicht 
zur quantitativen Beschreibung von Strömungsgrößen im Rezirkulationsgebiet 
technisch relevanter Nachlaufströmungen verwendbar (Murthy und Osborn, 
1976). 


Deshalb wurden beispielsweise von Lamb und Oberkampf (1995) semi- 
empirische Korrelationen für den Basisdruck und die Temperatur im Rezir- 
kulationsgebiet entwickelt. Als empirische Datengrundlage hatten Lamb und 
Oberkampf (1995) größtenteils von Sieling und Page (1970) zusammengefasste 
Messergebnisse aus Freiflugversuchen sowie Windkanalexperimenten gedient, 
bei denen keine Halterungsstrukturen im Nachlauf der Modelle verwendet 
worden waren. Lamb und Oberkampf (1995) hatten sich wegen der technischen 
Anwendungen, für welche die Korrelation entwickelt wurde, auf Messergeb- 
nisse beschränkt, bei denen das Verhältnis von Grenzschichtdicke zu Modell- 
durchmesser gering ist. Auf diese Weise kann der Basisdruck für Nachlauf- 
strömungen mit einer turbulenten Grenzschicht ausschließlich abhängig von 
der lokalen Machzahl der Strömung vor der Ablösung an der Modellhinterkante 
näherungsweise beschrieben werden (Lamb und Oberkampf, 1995). 
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2.1.3 Experimentelle Analysen und 
numerische Simulationen 


Die Verwendung von Methoden zur Vermessung des Geschwindigkeitsfeldes 
mit geringer Beeinflussung der Strömung wie der Laser-Doppler-Anemometrie 
(LDA) durch Leopold (1993) und Herrin und Dutton (1994) sowie der Particle- 
Image-Velocimetry-Methode (PIV) durch Scarano und van Oudheusden 
(2003) hatten zu präzisen Messergebnissen für die Geschwindigkeitsfelder 
kompressibler Nachlaufströmungen geführt. Die somit verfügbaren quantita- 
tiven Messfelder wurden zur Validierung numerischer Simulationen benötigt 
und werden bis heute zu diesem Zweck verwendet (Fureby et al., 1999, Baurle 
et al., 2001, Forsythe et al., 2002a,b, Fureby und Kupiainen, 2003, Kawai und 
Fujii, 2005, Barone und Roy, 2006, Simon et al., 2006, Sandberg und Fasel, 
2006a,b,c, Sandberg, 2012). 


Putnam und Bissinger (1985) fassten die Simulationen aus den späten 70ern 
und frühen 80ern zusammen. Sie schlussfolgerten, dass die damals verwen- 
deten Simulationsmethoden nicht in der Lage waren, den Druck an der Basis 
ausreichend genau zu berechnen. Die Simulationsmethoden basierten auf 
den Reynolds-gemittelten Navier-Stokes-Gleichungen (RANS) und model- 
lierten sämtliche turbulenten Skalen. Durch die Verbesserung der Modelle 
zur Beschreibung des Einflusses der Turbulenz auf die Strömung wurde eine 
genauere Berechnung des mittleren Basisdrucks ermöglicht, jedoch ergaben 
sich noch immer Abweichungen im Geschwindigkeitsfeld des Nachlaufes. 
Einen Überblick darüber gibt die Zusammenfassung von Dutton et al. (1995), 
in welcher die Verwendung von Large-Eddy-Simulationen (LES) oder Direkten 
Numerischen Simulationen (DNS) vorgeschlagen wird, um die Probleme bei 
der Modellierung der Turbulenz zu umgehen. 


Sandberg und Fasel (2006a,c) hatten den achsensymmetrischen Nachlauf 
der experimentellen Untersuchungen von Herrin und Dutton (1994) für eine 
Anströmung mit einer Machzahl M von 2,46 mittels DNS-Methode berechnet, 
welche sämtliche Wirbelskalen auflöst und somit keine Turbulenzmodellierung 
benötigt. Um die hohen Rechenzeiten der DNS-Methode zu reduzieren, hatten 
Sandberg und Fasel (2006a,c) eine im Vergleich zu den experimentellen 
Untersuchungen von Herrin und Dutton (1994) um den Faktor 33 reduzierte 
Reynoldszahl Rep, basierend auf dem Modelldurchmesser D, verwendet. 
Außerdem wurde im Gegensatz zu den experimentellen Bedingungen eine 
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laminare Grenzschicht simuliert. Trotz der unterschiedlichen Strömungs- 
bedingungen zeigen die Simulationen von Sandberg und Fasel (2006a,c) 
für einen Halbzylinder mit einer Reynoldszahl Rep von 1-10% ähnliche 
Strömungsmuster, wie sie von Bourdon und Dutton (1999) bei Rep = 3,3 - 10° 
und mit einer turbulenten Grenzschicht experimentell gezeigt worden waren. 
Insgesamt hatten sich jedoch deutliche quantitative Unterschiede zu dem von 
Herrin und Dutton (1994) vermessenen Stromungsfeld ergeben. Abbildung 2.1 
zeigt den Basisdruck p» relativ zum statischen Druck der Anströmung p. in 
Abhängigkeit von der mit dem Modellradius R = D/2 entdimensionalisierten 
radialen Position r. Es ist zu erkennen, dass die DNS-Simulation von Sandberg 
und Fasel (2006a) den Basisdruck im Vergleich zu den experimentellen Ergeb- 
nissen von Herrin und Dutton (1994) um 30 % überschätzte. Abbildung 2.2 
zeigt den Verlauf der mit der Anströmgeschwindigkeit U” normierten Axial- 
geschwindigkeit U, entlang der Zentralachse x. Es ist zu erkennen, dass 
sich in der Simulation von Sandberg und Fasel (2006a) ein im Vergleich zu 
den Messergebnissen von Herrin und Dutton (1994) etwa doppelt so langes 
Rezirkulationsgebiet ergeben hatte. Später hatte Sandberg (2012) den gleichen 
Referenzfall von Herrin und Dutton (1994) als Vollzylinder mit einer turbu- 
lenten Grenzschicht am Einlass des Simulationsgebietes simuliert und konnte 
auf diese Weise die Übereinstimmung mit den experimentellen Ergebnissen 
verbessern, obwohl auch hier die gegenüber den Experimenten reduzierte 
Reynoldszahl Rep von 1 - 10° verwendet worden war. 


Um Simulationen bei den hohen Reynoldszahlen der experimentellen Unter- 
suchungen von Herrin und Dutton (1994) oder auch Leopold (1993) durch- 
führen zu können, müssen trotz der gestiegenen Rechenleistungen bis heute 
Teile der Turbulenz modelliert werden (Spalart, 2000, Spalart und Venkata- 
krishnan, 2016). Die Abhängigkeit des Rechenaufwands vom Modellierungs- 
grad häufig verwendeter Methoden zur Turbulenzbeschreibung ist in Abbil- 
dung 2.3 schematisch dargestellt. 


Fureby et al. (1999) hatten eine LES-Methode auf den von Herrin und Dutton 
(1994) untersuchten achsensymmetrischen Überschallnachlauf angewendet. 
Bei einer LES-Methode werden die großskaligen Wirbel explizit berechnet 
und nur die kleinsten turbulenten Skalen werden mit einem Turbulenz- 
modell beschrieben (Sagaut, 2002). Die LES-Simulation von Fureby et al. 
(1999) hatte den gemessenen Basisdruck um etwa 6% und die Länge des 
Rezirkulationsgebietes um 5 % unterschätzt. Fureby et al. (1999) hatten die 
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Abbildung 2.1: Vergleich des radialen Basisdruckverlaufes unterschiedlicher numerischer Strö- 
mungssimulationen mit den experimentellen Ergebnissen von Herrin und Dutton (1994) für 


M = 2,46. 
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—— —— DNS, laminar, Halbzylinder (Sandberg und Fasel, 2006a) 
DNS, turbulent, Vollzylinder (Sandberg, 2012) 

—-—-— LES, grobes Gitter (Fureby et al., 1999) 

oo LES, feines Gitter (Fureby und Kupiainen, 2003) 

— - - - — MILES (Forsythe et al., 2002b) 

— — — DES, SA (Forsythe et al., 2002a) 

— — — DES, k-w-SST (Forsythe et al., 2002a) 

teens RANS, SA (Forsythe et al., 2002a) 

rennen RANS, k-w-SST (Forsythe et al., 2002a) 

_ Hybrid-RANS/LES, fein (Kawai und Fujii, 2005) 


Abbildung 2.2: Vergleich der Axialgeschwindigkeiten entlang der Zentralachse des Strömungs- 
feldes unterschiedlicher numerischer Simulationen mit den experimentellen Ergebnissen von 


Herrin und Dutton (1994) für M = 2,46. 
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DNS 


LES und 
MILES 


Hybride RANS/LES- 
Methoden (z.B. DES) 


RANS 


Abbildung 2.3: Schematischer Zusammenhang zwischen Rechenaufwand und Modellierungsgrad 
verschiedener Turbulenzmodellierungsverfahren. 


Unterschiede zwischen Experiment und Simulation mit der im Vergleich zum 
Experiment von Herrin und Dutton (1994) zu geringen Grenzschichtdicke 
vor der Strömungsablösung an der Modellhinterkante erklärt. Forsythe et al. 
(2000) hatten die zu geringe Grenzschichtdicke wiederum mit der für eine LES- 
Methode nicht ausreichenden Auflösung des Rechengitters in der Grenzschicht 
erklärt. Später hatten Fureby und Kupiainen (2003) auf einem Rechengitter 
mit der vierfachen Anzahl an Gitterzellen einen im Vergleich zum Experiment 
von Herrin und Dutton (1994) um 4% erhöhten Basisdruck berechnet und 
auf diese Weise die Übereinstimmung der berechneten Strömungsgeschwin- 
digkeiten entlang der zentralen Achse mit den experimentellen Ergebnissen 
verbessert, wie in Abbildung 2.2 gezeigt. Die Ergebnisse der zusätzlich mittels 
Monotone-Integrated-LES-Methode (MILES) durchgeführten Simulationen, 
bei der anstelle eines Turbulenzmodells die dissipativen Eigenschaften der 
numerischen Diskretisierungsmethode zur Modellierung der kleinsten Turbu- 
lenzskalen verwendet werden (Boris, 1990, Grinstein et al., 2011), zeigen nur 
geringfügige Unterschiede zu den LES-Simulationen. Daraus hatten Fureby 
et al. (1999) sowie Fureby und Kupiainen (2003) geschlussfolgert, dass die 
Turbulenzmodellierung auf den kleinsten Skalen nur einen geringen Einfluss 
auf die Simulationsergebnisse hat. Baurle et al. (2001) hatten den Referenzfall 
von Herrin und Dutton (1994) ebenfalls mit einer MILES-Methode berechnet. 
Im Gegensatz zu den Ergebnissen von Fureby und Kupiainen (2003) unter- 
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schätzen die Ergebnisse von Baurle et al. (2001) den von Herrin und Dutton 
(1994) gemessenen Basisdruck im Mittel um 5%. Die Simulationen von 
Forsythe et al. (2002b), die ebenfalls mittels MILES-Methode berechnet 
worden waren, zeigen sowohl für den Basisdruck als auch für die Geschwin- 
digkeitsverteilung entlang der Zentralachse eine sehr gute Übereinstimmung 
mit den experimentellen Ergebnissen. Forsythe et al. (2002b) hatten dafür 
jedoch ein feineres Rechengitter als Fureby und Kupiainen (2003) und Baurle 
et al. (2001) verwendet. 


Forsythe et al. (2002a,b) hatten außerdem eine Detached-Eddy-Simulations- 
methode (DES) angewendet, bei der im Vergleich zur LES-Methode nicht nur 
die kleinsten turbulenten Skalen modelliert werden, sondern auch die gesamte 
Turbulenz in der Nähe von festen Oberflächen (Spalart, 2001). Deshalb kann im 
Allgemeinen ein gröberes Rechengitter verwendet werden, sodass der Rechen- 
aufwand gegenüber einer LES-Methode reduziert wird. Die Ergebnisse der 
DES-Simulationen hatten Forsythe et al. (2002a,b) mit RANS- und MILES- 
Ergebnissen verglichen. Wie in Abbildung 2.1 gezeigt, waren die RANS- 
Simulationen nicht in der Lage, den von Herrin und Dutton (1994) experi- 
mentell bestimmten Basisdruckverlauf wiederzugeben (Fasel und Sandberg, 
2006). Auch im Geschwindigkeitsfeld entlang der zentralen Achse ergaben 
sich deutliche Unterschiede. Im Vergleich zu den experimentellen Ergebnissen 
von Herrin und Dutton (1994) ist die maximale Rückströmgeschwindigkeit 
etwa doppelt so hoch und weist einen geringeren Abstand von der Basis auf. 
Im Gegensatz zu den RANS-Simulationen waren die Simulationen mittels 
DES-Methode in der Lage, den Basisdruck qualitativ korrekt wiederzugeben 
und weichen um weniger als 6 % von den experimentellen Ergebnissen ab. Die 
DES-Simulation mit dem Wirbelviskositätsmodell von Spalart und Allmaras 
(1992) (SA) zeigt im Vergleich zu der DES-Simulation mit dem k-w-Shear- 
Stress-Transport-(SST)-Modell nach Menter (1994) (k-w-SST) eine geringere 
Abweichung von weniger als 3 % gegenüber den experimentellen Ergebnissen 
von Herrin und Dutton (1994). Die Abhängigkeit der DES-Ergebnisse von 
der Gitterauflösung war weniger ausgeprägt als bei den Simulationen mittels 
MILES-Methode. Aufgrund der geringeren Abhängigkeit der Ergebnisse von 
der Gitterauflösung und des verminderten Rechenaufwandes durch die Turbu- 
lenzmodellierung in der Grenzschicht hat sich die DES-Methode gegenüber 
LES- und MILES-Methoden als effizientere Methode zur Simulation von 
Überschallnachlaufströmungen erwiesen (Forsythe et al., 2002b). Kawai und 
Fujii (2005) hatten ebenfalls eine hybride RANS/LES-Methode auf den 
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Referenzfall von Herrin und Dutton (1994) angewendet und konnten die gute 
Übereinstimmung des berechneten Basisdruckes mit den experimentellen 
Ergebnissen bei einem im Vergleich zu LES- und MILES-Methoden gerin- 
geren Rechenaufwand bestätigen. Allerdings weichen die numerischen Ergeb- 
nisse auf dem feinsten Gitter (superfein) deutlicher von den experimentellen 
Ergebnissen ab als auf dem gröberen Gitter (fein). Die Ergebnisse der unter- 
schiedlichen hybriden RANS/LES-Methoden (Forsythe et al., 2002a, Kawai 
und Fujii, 2005) zeigen in Abbildung 2.2 eine systematische Abweichung der 
Länge des Rezirkulationsgebietes von mehr als 10 % gegenüber den experi- 
mentellen Ergebnissen von Herrin und Dutton (1994). Baurle et al. (2001) 
sowie Forsythe et al. (2002a) hatten die Abweichungen in der Länge des Rezir- 
kulationsgebietes auf den Übergang zwischen der modellierten Turbulenz der 
Grenzschicht und der aufgelösten Turbulenz in der Scherschicht zurückge- 
führt. Da die bereits in der Grenzschicht existierenden Wirbelstrukturen bei 
hybriden RANS/LES-Simulationen erst im LES-Bereich aufgelöst werden 
können (Baurle et al., 2001), werden die Instabilitäten der Scherschicht unmit- 
telbar hinter der Modellhinterkante von hybriden RANS/LES-Simulationen im 
Allgemeinen nicht korrekt wiedergegeben (Forsythe et al., 2002a). Dies führt 
zu einem veränderten Wachstum der Scherschicht (Forsythe et al., 2002a) und 
resultierte in dem in Abbildung 2.2 zu erkennenden Überschätzen der Länge 
des Rezirkulationsgebietes im Vergleich zu den experimentellen Ergebnissen 
von Herrin und Dutton (1994). Eine Übersicht über numerische Simulations- 
ergebnisse für Überschallnachlaufströmungen mit unterschiedlichen Simulati- 
onsmethoden zeigen auch Simon et al. (2006). 


Insgesamt konnte gezeigt werden, dass die teilweise turbulenzauflösenden 
Simulationsverfahren (LES, MILES, DES, Hybrid-RANS/LES) in der Lage 
sind, die Messergebnisse von Herrin und Dutton (1994) qualitativ gut wieder- 
zugeben. Gegenüber den Messergebnissen ergeben sich jedoch die in Ab- 
bildung 2.1 gekennzeichneten quantitativen Abweichungen des Basisdrucks 
von +6%, verursacht durch die unterschiedlichen Parameter der turbulenz- 
auflösenden numerischen Methoden. Die DES-Methode sowie vergleichbare 
hybride RANS/LES-Methoden überschätzen dabei die Länge des Rezirkulati- 
onsgebietes systematisch. 


2 Stand der Forschung 


2.1.4 Beeinflussung von Nachlaufströmungen 


Aufgrund des absolut instabilen Strömungsbereichs innerhalb des Rezirkulati- 
onsgebietes (Leopold, 1993) können selbst kleine Störungen in diesem Bereich 
das Strömungsfeld grundlegend beeinflussen. Dies ist für die Vermessung 
des Strömungsfeldes hinderlich, kann aber auch dazu genutzt werden, das 
Strömungsfeld und dessen Eigenschaften gezielt zu verändern. 


Eine Möglichkeit der Einflussnahme auf die Nachlaufströmung besteht darin, 
die Grenzschicht bereits vor der Ablösung zu verändern. Augenstein et al. 
(1999) sowie Durgesh et al. (2013) hatten zum Beispiel den Einfluss unter- 
schiedlicher Oberflächenrauigkeiten an der Modelloberfläche auf die turbu- 
lente Scherschicht gezeigt. Durgesh et al. (2013) hatten festgestellt, dass eine 
mit einer erhöhten Oberflächenrauigkeit einhergehende dickere Grenzschicht 
einen höheren Basisdruck zur Folge hat. Dies wurde auch von Lamb und 
Oberkampf (1995) beschrieben. 


Die experimentellen Studien von Bourdon und Dutton (2001, 2002) sowie 
Janssen und Dutton (2005) hatten den Einfluss von in der Grenzschicht 
erzeugten Wirbeln auf die turbulenten Strukturen in der resultierenden Scher- 
schicht und auf das Rezirkulationsgebiet untersucht. Es konnte eine Verän- 
derung der turbulenten Durchmischung in der Scherschicht festgestellt werden, 
die zu einem veränderten Druck im Rezirkulationsgebiet führte. Die verwen- 
deten Wirbelgeneratoren befanden sich vollständig innerhalb der Grenzschicht 
und waren teilweise in der Lage, den Basisdruck um bis zu 3 % (Bourdon und 
Dutton, 2002, Janssen und Dutton, 2005) zu erhöhen. 


Außerdem besteht die Möglichkeit, die Strömung durch eine veränderte 
geometrische Formgebung zu beeinflussen. Eine häufig genutzte geometrische 
Modifikation ist das Anbringen eines Kegelstumpfes (boat-tail) hinter dem 
zylindrischen Körper (Cortright und Schroeder, 1951b). Das Rezirkulations- 
gebiet kann auf diese Weise verkleinert und die Strömung bereits vor der 
Ablösung in Richtung der Achse umgelenkt werden. 


Gasausblasungen (base bleed) an der Modellbasis (Cortright und Schroeder, 
1951a) und Verbrennungsprozesse (base combustion) im Rezirkulationsgebiet 
(Clayden und Bowman, 1968) verändern die Strömungsstruktur des Nachlaufs 
grundlegend und beeinflussen somit ebenfalls den Basisdruck. Ward und 
Baltakis (1976) hatten den Einfluss rotierender Gasausblasungen auf dieses 
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veränderte Strömungsfeld untersucht und konnten eine Erhöhung des Basis- 
druckes in Abhängigkeit von der Rotationsrate des ausströmenden Gases 
messen, da sich der maximale Massenstrom in das Rezirkulationsgebiet durch 
die Rotation erhöht hatte. Zu dem gleichen Ergebnis kamen Kayser et al. 
(1988), die Gasausblasungen bei einem rotierenden Modell untersucht hatten. 


Numerische Untersuchungen von Sandberg und Fasel (2006a) zeigen, dass 
das Eliminieren azimutaler Wirbelmoden niedriger Ordnung im Nachlauf 
durch geometrische Beschränkung des Simulationsgebietes den Basisdruck 
erhöhen kann. Experimentell konnte dies durch die Verwendung von Trenn- 
elementen im Rezirkulationsgebiet nicht gezeigt werden, jedoch wurden die 
Druckschwankungen an der Basis durch die Trennelemente um bis zu 39 % 
vermindert (Reedy et al., 2012). 


Dupuis und Hathaway (1997) hatten den Einfluss unterschiedlich stark 
angestellter Leitflächen auf die aerodynamischen Kräfte untersucht, die auf 
ein Freiflugmodell wirken. Die Anströmung versetzte die Modelle je nach 
Anstellwinkel der Leitflächen in Rotation. Die Modelle rotierten anschließend 
nahezu mit Synchrondrehzahl, sodass nur ein geringer Drehimpuls in das 
Nachlaufgebiet eingebracht worden war. Deshalb konnten keine signifikanten 
Veränderungen des Nachlaufs beobachtet werden. Im Gegensatz dazu konnten 
Miller und Molnar (1986) für eine transsonische Anströmung eine geringfügige 
Veränderung des Basisdruckes eines Freiflugmodells aufgrund von Rotation 
zeigen. Silton und Weinacht (2008) hatten im Überschall einen geringen 
Anstieg des Gesamtströmungswiderstandes für rotierende Strömungskörper 
mit an der Oberfläche in axialer Richtung eingebrachten Einkerbungen gegen- 
über glatten Strömungskörpern gezeigt. Die Einkerbungen könnten gegenüber 
den glatten Modellen ein höheres Drehmoment an die Strömung übertragen 
haben, sodass die unterschiedlichen Strömungswiderstände auf einen Einfluss 
des Dralls zurückgeführt werden könnten. Silton und Weinacht (2008) hatten 
jedoch nicht gezeigt, ob sich der Stromungswiderstand aufgrund veränderter 
Reibungskräfte an der Modelloberfläche, eines an die Strömung übertragenen 
Dralls oder anderer Gründe erhöht hatte. 


Hruschka und Leopold (2015) hatten in einem Windkanalexperiment den 
Einfluss rotierender Leitflächen auf den Basisdruck und das Rezirkulations- 
gebiet längsangeströmter zylindrischer Strömungskörper untersucht. Während 
die unmittelbar vor der Modellhinterkante angebrachten Leitflächen mit einem 
eingebauten Antrieb in Rotation versetzt wurden, rotierte die Modellbasis nicht. 
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Bei einer Drehrate von 80 Hz hatte sich der zentrale Basisdruck gegenüber dem 
nichtrotierenden Strömungskörper um 1,6 % erhöht. Für die maximal erreichte 
Drehrate von 170 Hz hatten Hruschka und Leopold (2015) eine Abnahme des 
zentralen Basisdruckes gegenüber dem statischen Modell von 1,5 % festgestellt. 
Die Rauchvisualisierungen des Nachlaufs zeigten jedoch keine strukturellen 
Veränderungen. Die zusätzlich durchgeführten numerischen Simulationen, 
bei denen die Modellbasis ebenfalls nicht rotierte, konnten für 160 Hz den 
Druckabfall gegenüber dem nichtrotierenden Modell wiedergeben, zeigten 
jedoch für 80 Hz keinen Druckanstieg. Für Drehraten von 320 Hz und 480 Hz 
hatte sich in den numerischen Ergebnissen entlang der Zentralachse ein strom- 
abgerichteter Wirbelschlauch ausgebildet, der von der Modellbasis bis in den 
fernen Nachlauf reichte und zu einem Basisdruckabfall von bis zu 65 % führte. 


Die numerischen Simulationen von Hruschka und Leopold (2015) zeigen 
bisher als einzige eindeutig den Einfluss eines Dralls auf Basisdruck und 
Strömungsstruktur einer Überschallnachlaufströmung. Die stabilitätstheore- 
tische Untersuchung von Jim&nez-Gonzälez et al. (2014) für eine laminare, 
inkompressible Strömung lassen außerdem eine Veränderung im Nachlauf 
erwarten, die einem Wirbelaufplatzen (Peckham und Atkinson, 1957, Hall, 
1972) ähnlich ist. Experimentelle Nachweise für Änderungen der Strömungs- 
struktur kompressibler Nachlaufströmungen aufgrund eines eingebrachten 
Dralls wurden in der Literatur jedoch bisher nicht beschrieben. 


2.2 Drallstromungen 


2.2.1 Wechselwirkung von Wirbelschläuchen 
mit Verdichtungsstößen 


In der Literatur existieren allerdings Untersuchungen zu rotierenden Über- 
schallströmungen, die sich nicht unmittelbar hinter Strömungskörpern befin- 
den. Délery et al. (1984) hatten gezeigt, dass die Strömungsstruktur von Wirbel- 
schläuchen durch die Interaktion mit Verdichtungsstößen beeinflusst werden 
kann. Delery et al. (1984) konnten damit erstmals eine Art Wirbelaufplatzen in 
Überschallströmungen beobachten, was zuvor ausschließlich für inkompres- 
sible Strömungen beschrieben worden war (Peckham und Atkinson, 1957). 
Für inkompressible Strömungen hatte Hall (1972) Wirbelaufplatzen als eine 
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plötzliche Veränderung der Wirbelstruktur beschrieben, bei der die Strömung 
entlang der Wirbelachse stark verlangsamt und in radialer Richtung abgelenkt 
wird. Auch Rusak et al. (2015) hatten Wirbelaufplatzen als ein schlagartig 
auftretendes Phänomen bezeichnet, welches einen Wirbelschlauch in einen 
energetisch günstigeren Zustand einer um ein nahezu stagnierendes Totwasser- 
gebiet rotierenden Strömung überführt. In der Literatur sind mehrere mögliche 
Erklärungen für dieses Strömungsphänomen zu finden. Benjamin (1962) hatte 
Wirbelaufplatzen beispielsweise mittels stabilitätstheoretischer Überlegungen 
als eine sprunghafte Transition von einem überkritischen in einen unter- 
kritischen Zustand beschrieben. Auch in neueren Zusammenfassungen über 
Wirbelaufplatzen, wie der von Lucca-Negro und O’Doherty (2001), werden 
allerdings noch immer verschiedene Erklärungsversuche für das Auftreten 
dieses Phänomens dargestellt. In der Literatur ist bis heute selbst für inkom- 
pressible Strömungen keine konsistente Beschreibung dieses Phänomens zu 
finden (Rusak et al., 2015). 


Visbal und Gordnier (1995) hatten für kompressible Unterschallstromungen 
festgestellt, dass das Aufplatzen eines Wirbelschlauches bei höheren Mach- 
zahlen höhere Drallraten erfordert. Dies wurde auch durch die numeri- 
schen Simulationsergebnisse von Kuruvila und Salas (1990), Melville (1996), 
Herrada et al. (2003) bestätigt. Die Simulationsergebnisse von Luginsland und 
Kleiser (2015) zeigen hingegen den gegenteiligen Effekt, dass mit steigender 
Machzahl die kritische Drallrate abnimmt. Eine globale Stabilitätsanalyse 
von kompressiblen Unterschall-Wirbelströmungen von Rusak et al. (2015) 
zeigt wiederum einen Anstieg der kritischen Drallrate mit ansteigender axialer 
Machzahl. 


Metwally et al. (1989) hatten für Überschallströmungen zwischen einer schwa- 
chen und einer starken Wechselwirkung zwischen Verdichtungsstößen und 
Wirbelschläuchen differenziert, wobei nur bei der starken Wechselwirkung 
eine Art Wirbelaufplatzen zu beobachten ist. Schon Delery et al. (1984) hatten 
mit der Machzahl der äußeren Strömung normal zum Stoß und der Drehrate des 
Schlauchwirbels zwei Parameter zur Charakterisierung des Wirbelaufplatzens 
in Überschallströmungen beschrieben, die bis heute verwendet werden. Catta- 
festa und Settles (1992) hatten die von Delery et al. (1984) ermittelten kriti- 
schen Bedingungen für Wirbelaufplatzen in Überschallströmungen auf einen 
größeren Machzahlbereich normal zum Verdichtungsstoß erweitert. Settles 
und Cattafesta (1993) hatten anschließend die beiden Untersuchungen und 
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die daraus resultierenden Grenzbedingungen für Wirbelaufplatzen zusammen- 
gefasst. Sowohl Delery et al. (1984) als auch Settles und Cattafesta (1993) 
hatten für Überschall-Wirbelschläuche eine Abnahme der kritischen Drallrate 
mit steigender Machzahl M, normal zum Stoß festgestellt. Für einen unendlich 
schwachen Stoß mit M, = 1 hatten Settles und Cattafesta (1993) den kriti- 
schen Winkel zwischen Axial- und Azimutalgeschwindigkeit zu etwa 30° 
bestimmt. Für inkompressible Strömungen (M, — 0) findet Wirbelaufplatzen 
bei höheren Strömungswinkeln von etwa 40° statt (Hall, 1972). Ausgelöst wird 
das Wirbelaufplatzen in inkompressiblen Strömungen oftmals durch einen in 
Strömungsrichtung ansteigenden Druckgradienten. In Überschallströmungen 
wird Wirbelaufplatzen oft von Verdichtungsstößen ausgelöst, die ebenfalls 
einen ansteigenden Druckgradienten darstellen. Somit sind nicht nur die 
beobachteten Strömungsstrukturen sondern auch die auslösenden Faktoren für 
Wirbelaufplatzen in inkompressiblen und supersonischen Strömungen ähnlich. 


2.2.2 Einfluss des Dralls auf turbulente 
Mischungsprozesse 


In der Literatur wurde darüber hinaus der Einfluss des Dralls auf die turbu- 
lente Durchmischung zweier koaxialer Strömungen beschrieben. Charakte- 
ristisch für Jetströmungen, bei denen die Axialgeschwindigkeit der beiden 
koaxialen Strömungen im Allgemeinen unterschiedlich ist, ist die sich ausbil- 
dende Mischungsschicht zwischen den beiden Strömungen. Die strömungs- 
mechanischen Prozesse in den Mischungsschichten von Jet- und Nachlauf- 
strömungen sind miteinander vergleichbar, sodass Effekte, die für einen der 
beiden Strömungsfälle beschrieben werden, im Allgemeinen auch für den 
anderen beobachtbar sind (Demetriades, 1968, 1976). Beispielsweise treten 
wellenartige großskalige Turbulenzstrukturen innerhalb der Scherschicht 
sowohl in Jet- (Ponton und Seiner, 1995) als auch in Nachlaufströmungen 
(Bourdon und Dutton, 1999, 2000) auf und können dort Machwellen erzeugen 
(Tam, 2009). Die durch diese Turbulenzstrukturen erzeugten hochfrequenten 
Druckschwankungen sind akustisch wahrnehmbar (Ponton und Seiner, 1995, 
Tam, 2009), verändern die zeitlich gemittelte Strömung jedoch kaum (Bourdon 
und Dutton, 2000) und werden deshalb in der vorliegenden Arbeit nicht unter- 
sucht. 
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Im Gegensatz dazu hatten Swithenbank und Chigier (1969) anhand von 
Experimenten mit rotierenden Jetströmungen im Trans- und Unterschall 
(Moore, 1964, Chigier und Chervinsky, 1967) mittels theoretischer Über- 
legungen eine im zeitlichen Mittel erhöhte turbulente Durchmischung für 
rotierende Überschall-Jetströmungen geschlussfolgert. Dieser Effekt konnte 
jedoch weder durch die Experimente von Povinelli und Ehlers (1972) noch 
von Schetz und Swanson (1973) nachgewiesen werden. Mutmaßlich war der 
erzeugte Drall, der nicht näher spezifiziert worden ist, zu gering, um einen 
messbaren Einfluss auf die Durchmischung zu haben. Erst spätere Experimente 
von Swithenbank et al. (1989), Naughton et al. (1989), Tillman et al. (1991) 
und Naughton und Settles (1992) hatten die Theorie von Swithenbank und 
Chigier (1969) bestätigen können. 


Die Änderung der turbulenten Durchmischung kann mit veränderten turbu- 
lenten Strukturen in der Mischungsschicht (Naughton und Settles, 1992) 
und dem Stabilitätsverhalten der Strömung (Panda und McLaughlin, 1994) 
erklärt werden. Ein ähnliches Verhalten konnte von Bourdon und Dutton 
(2001, 2002) sowie Janssen und Dutton (2005) auch für Nachlaufströmungen 
beobachtet werden. Bourdon und Dutton (2001) hatten gezeigt, dass in 
Strömungsrichtung ausgerichtete turbulente Strukturen in der Scherschicht 
die Durchmischung der Strömung im Rezirkulationsgebiet mit der äußeren 
Strömung erhöhen. In der Folge wird im Vergleich zu einer unbeeinflussten 
Nachlaufströmung (Bourdon und Dutton, 1999) mehr Fluid aus dem Rezir- 
kulationsgebiet entnommen, sodass sich der Basisdruck vermindert (Bourdon 
und Dutton, 2001). Im Gegensatz dazu hatten Bourdon und Dutton (2002) für 
in azimutaler Richtung eingebrachte turbulente Strukturen gegenüber einer 
unbeeinflussten Nachlaufströmung ein geringeres Wachstum der Scherschicht- 
dicke zwischen der Ablösung und der Rekompression der Strömung feststellen 
können. Dies deutet auf eine geringere Durchmischung des Rezirkulations- 
gebietes mit der äußeren Strömung hin und führte zu einem um bis zu 3 % 
erhöhten Basisdruck. 
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3 Studie mit rotierenden 
Strömungskörpern 


Zunächst wurde eine Studie mit rotierenden Strömungskörpern durchge- 
führt. Technisch bedingt waren die Drehraten der im Folgenden beschrie- 
benen experimentellen Untersuchungen allerdings nicht ausreichend, um eine 
Änderung der Nachlaufstruktur hervorzurufen. Es konnte jedoch die von 
Hruschka und Leopold (2015) mittels numerischer Simulationen gezeigte Än- 
derung des zentralen Basisdrucks experimentell bestätigt werden. Außerdem 
diente diese Studie der Entwicklung einer alternativen Vorgehensweise zur 
Erzeugung höherer Drallraten im Nachlauf. 


3.1 Versuchsaufbau und Anströmbedingungen 


Modellhalterung 


Strömungskörper 


Befestigung 


-_— [22 


Pre 
Ruhekammer Lavaldüse Messkammer 


Abbildung 3.1: Skizze des verwendeten Versuchsaufbaus in einem Überschallwindkanal (Weidner 
et al., 2016, 2017). 
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| # Verdichtungsstoß 
am Ende der 
Lavaldüse 


Abbildung 3.2: Lavaldüse und Messkammer mit Modellhalterung, Strömungskörper und einge- 
zeichneten Verdichtungsstößen am Ende der Lavaldüse. 


Zur Untersuchung des Einflusses von Drall auf Überschallnachlaufströmungen 
wurde im Folgenden der in Abbildung 3.1 skizzierte Versuchsaufbau verwen- 
det. Ähnliche Versuchsaufbauten sind bereits in früheren Studien (Sieling und 
Page, 1970, Leopold, 1993, Herrin und Dutton, 1994, Augenstein et al., 1999) 
zur Untersuchung von Nachlaufstromungen ohne Drall eingesetzt worden. 


Aus der mit Drucktanks verbundenen Ruhekammer strömte trockene Luft mit 
einem nahezu konstanten Ruhedruck po von 4,8 - 10° Pa und einer mittleren 
Ruhetemperatur 79 von 295K aus. Aufgrund der überkritischen Durch- 
strömung der an die Ruhekammer anschließenden Lavaldüse wurde in der in 
Abbildung 3.2 gezeigten Messkammer, mit einem quadratischen Querschnitt 
von 0,2mx0,2m, eine Überschallströmung mit der Machzahl M» = 2 
erreicht. 


Die zu untersuchenden Strömungskörper wurden an einer zentral im Windkanal 
montierten, achsensymmetrischen Modellhalterung befestigt. Um Stöße in der 
Anströmung der Modelle zu vermeiden, wurde die Modellhalterung, wie in 
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Abbildung 3.1 skizziert, in der Ruhekammer befestigt. Auf diese Weise war 
das Nachlaufgebiet frei von Halterungsstrukturen und durch diese verursachte 
Verdichtungsstöße, die die Nachlaufströmung unmittelbar beeinflussen können 
(Sieling und Page, 1970). Außerdem wurde das Modell so im Windkanal 
montiert, dass die an dem nicht glatten Übergang zwischen Lavaldüse und 
Messkammer entstehenden und in Abbildung 3.2 eingezeichneten Verdich- 
tungsstöße den Nachlauf so gering wie möglich beeinflussten (Chapman, 
1950). 


3.1.1 Messung des statischen Drucks der Anströmung 


radiale Bohrung 
zur Messung des 
statischen Drucks 


Abbildung 3.3: Skizze zur Position der statischen Druckmessungen an der Oberfläche eines 
zylindrischen Modells. 


Der statische Druck im Freistrom wurde, wie in Abbildung 3.3 skizziert, 
38,5 mm stromauf der Modellhinterkante gemessen. Abbildung 3.4 zeigt die 
Messungen für vier azimutale Positionen y und vergleicht diese mit einer 
isentropen Expansion des Ruhedrucks po = 4,8: 10°Pa auf die nominelle 
Machzahl von M = 2. Die Fehlerbalken geben die statistische Messabwei- 
chung der verwendeten Druckmessgeber an. Der arithmetische Mittelwert 
der Druckmessungen von 0,605 - 10° Pa ist um 1,3 % geringer als der theore- 
tische Druck von 0,613 - 10° Pa für eine isentrope Expansion auf die nominelle 
Machzahl von M = 2. Aus dem Verhältnis des gemessenen statischen Drucks 
zum Ruhedruck po wurde die Freistrom-Machzahl experimentell zu M = 2,01 
bestimmt. 
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Abbildung 3.4: Messung des statischen Freistromdrucks an der Modelloberfläche in Abhängigkeit 
von der Azimutalkoordinate 4. 


3.1.2 Messung der Grenzschichtdicke 


Der Gesamtdruck stromauf der Modellhinterkante wurde mittels einer Pitot- 
Sonde vermessen. Die Pitot-Sonde wurde, wie in Abbildung 3.5 skizziert, 
38,5 mm vor der Modellhinterkante angebracht. Der Einlass der Sonde befand 
sich 45,5 mm vor der Modellhinterkante und der Abstand Ar des Einlasses von 
der Modelloberfläche konnte beliebig eingestellt werden. Unmittelbar vor dem 
Einlass der Pitot-Sonde, mit einem Außendurchmesser von 1,25 mm, bildete 
sich der in Abbildung 3.5 sichtbare Verdichtungsstoß aus, der im Bereich des 
Einlasses als nahezu senkrecht betrachtet werden kann. 


Abbildung 3.6 zeigt die Messergebnisse für den Pitotdruck ppitot, die mit 
dem 38,5 mm vor der Modellhinterkante an der Modelloberfläche gemessenen, 
mittleren statischen Druck von 0,61 - 10° Pa normiert wurden. Für Ar > 5mm 
ist der gemessene relative Pitotdruck um 3 % geringer als der theoretische Wert 
von 5,64, der sich nach der Rayleigh-Pitot-Formel (Anderson, 2001) für eine 
Strömung mit M = 2 ergibt. Die systematisch geringeren Messwerte ergeben 
sich unter anderem daraus, dass der Innendurchmesser der Pitotsonde mit 
0,8mm im Vergleich zum Außendurchmesser von 1,25 mm relativ groß war 
(Tropea et al., 2007). Aus dem experimentell bestimmten Pitotdruck wurde 
die Machzahl in der äußeren Strömung zu mindestens 1,97 bestimmt und war 
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Verdichtungsstoß 


Abbildung 3.5: Skizze zum Aufbau der Pitotdruck-Messungen in der Grenzschicht mit Schlieren- 
visualisierung der Strömung um die Pitot-Sonde. 


12 pq 
| ©. Pitotdruck-Messungen O | 
L | —<— Numerische Simulation J 
10 F4 - - - Rayleigh-Pitot-Formel für M = 2 (Anderson, 2001) J 
gT 8F 
8 H 
x L 
< L 
Ean 
= 
E L 
= 4j} 
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ox [IA RER Oe TE RR RR BR BR DR DR KR DR DR DR RR DR OOO A HER 
1,0 1,5 2,0 2,9 3,0 3,5 4,0 4,5 5,0 5,5 6,0 


Relativer Pitotdruck ppitot/ Poo 


Abbildung 3.6: Vergleich des gemessenen und simulierten Pitotdrucks in der Grenzschicht mit 
dem theoretischen Pitotdruck gemäß Rayleigh-Pitot-Formel fiir M = 2 (Weidner et al., 2019b). 
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somit um bis zu 1,5 % geringer als die nominelle Machzahl Mœ% = 2. Die 
Abweichungen des experimentell bestimmten Pitotdrucks und der Machzahl 
der freien Anströmung von den theoretisch abgeleiteten Werten liegen damit 
innerhalb der erwarteten experimentellen Unsicherheiten. 


In Abbildung 3.6 ist außerdem das Ergebnis einer RANS-Simulation (An- 
hang A) des in Abbildung 3.1 gezeigten Versuchsaufbaus dargestellt. Der 
relative Pitotdruck der Simulationen wurde nach Sauer (1951) aus der 
ungestörten Machzahl in der Grenzschicht bestimmt. Der relative Pitot- 
druck im Freistrom der Simulationsergebnisse entspricht einer Machzahl der 
Strömung von 2,02. 


Sowohl die Ergebnisse der Messung als auch der Simulationen zeigen eine 
mittlere Grenzschichtdicke 6 von 5mm mit einer Unsicherheit von +1 mm 
und entsprechen damit den Messungen von Leopold (1993), die mit dem 
gleichen Versuchsaufbau jedoch fiir einen Modelldurchmesser von 38,66 mm 
durchgeführt wurden. Für den Modelldurchmesser, den auch Leopold (1993) 
verwendet hatte, hatten Augenstein et al. (1999) die Grenzschichtdicke mit 
einem speziell entwickelten Laser-Doppler-Anemometer (Grenzschicht-LDA) 
zu 4 mm bestimmt. 


Mit den Geschwindigkeitsmessungen innerhalb der Grenzschicht konnten 
Augenstein und Leopold (1999) außerdem zeigen, dass die Grenzschicht 
ein turbulentes Grenzschichtprofil aufwies. Abbildung 3.7 vergleicht die von 
Augenstein und Leopold (1999) gemessenen und die numerisch simulierten 
dimensionslosen Strömungsgeschwindigkeiten u* (Schlichting und Gersten, 
2006) in Abhängigkeit vom dimensionslosen Wandabstand y* (Schlichting 
und Gersten, 2006). Sowohl die Messungen von Augenstein und Leopold 
(1999) als auch die Simulationsergebnisse zeigen eine gute Ubereinstimmung 
mit dem logarithmischen Uberlappungsgesetz einer turbulenten Grenzschicht 
nach Schlichting und Gersten (2006). 


Abbildung 3.7 zeigt auBerdem, dass der dimensionslose Wandabstand y* des 
wandnächsten Zellzentrums in etwa 1 beträgt, sodass die reinviskose Unter- 
schicht von der Simulation aufgelöst wurde und nicht mittels Wandfunktionen 
modelliert werden musste. 
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Abbildung 3.7: Dimensionsloses Geschwindigkeitsprofil der Grenzschicht. 


3.2 Modellgeometrien und 
Stromungsbedingungen im Nachlauf 


An der in Abbildung 3.1 gezeigten Modellhalterung wurden die in Tabelle 3.1 
aufgelisteten und in Abbildung 3.8 skizzierten Strömungskörper montiert 
(Weidner et al., 2019b). Die verwendeten Strömungskörper unterschieden 
sich in der Anzahl nyp und der Höhe hyp der Leitflächen. Es wurden sowohl 
Leitflächen verwendet, die vollständig innerhalb der Grenzschicht lagen, als 
auch Leitflächen, die mehr als doppelt so hoch wie die mittlere Grenzschicht- 
dicke 6 = 5mm waren. Alle Leitflächen hatten eine Dicke von 1 mm und 
reichten von der Modellhinterkante 20 mm stromauf. 


Alle Modelle hatten auf der in Abbildung 3.2 eingezeichneten Länge Lp von 
135 mm vor der Modellhinterkante einen konstanten Modelldurchmesser D von 
40 mm. Aufgrund der Reynoldszahl Rez 5 = PooUooLs/ [Hoo > 8 - 10° basierend 
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Tabelle 3.1: Konfiguration der Leitflächen, gemittelter zentraler Basisdruck ohne Rotation, 
Anzahl der Experimente, maximal erreichte Modelldrehrate und relative Standardabweichung 
der Messungen ohne Rotation für die rotierenden Modelle (Weidner et al., 2019b). 


Modell ne Aurlmm] PE p-o (O5%CH Nep. JMmax [Hz] me 
Zylinder - = 0,634 +0,009 10 720 22% 
RI 12 15 0,634 + 0,006 20 765 2,0% 
R2 12 30 0,634 + 0,003 21 602 1,0% 
R3 12 60 0,615 + 0,014 7 417 2,4% 
R4 12 125 0,589 + 0,007 7 183 12% 
R5 6 125 0,578 + 0,004 8 212 0,7% 


R6 3 12,5 0,612 + 0,005 22 299 1,8 % 


Abbildung 3.8: Skizze der verwendeten rotierenden Strömungskörper mit geraden Leitflächen 
(Weidner et al., 2019b). 
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auf der Lauflänge der Grenzschicht Ls > Lp, der Strömungsgeschwindigkeit 
der Anströmung U~ von 520 m/s, der Dichte der Anströmung 
Po 


ae = 1,28kg/m, (3.1) 
A Rıuft1o el 


die sich aus dem Druck pæ = 0,61 - 10° Pa und der Temperatur Tœ = 166K 
in der Anströmung sowie der spezifischen Gaskonstanten für trockene 
Luft Riut = 287,1 J/(kg K) ergibt, und der dynamischen Viskosität in der 
Anströmung nach Sutherland 


co = 17,89 - 1076 P 
H °°°7_+110K \288K 


= 113 - 1076 Pas (3.2) 


3 
288K+110K/ Ts V 
DERI | (Anderson, 2006) 


lassen sich die untersuchten Nachlaufströmungen dem turbulenten Bereich 
zuordnen (Chapman, 1950, Kurzweg, 1951). 


3.3 Messungen des zentralen Basisdruckes 


Die Modelle wurden als Hohlkörper konzipiert, sodass an der Basis Druck- 
messungen mit einem nichtrotierenden Messsystem durchgeführt werden 
konnten. Der Druck p, im Zentrum der Basis wurde mit dem simultan gemes- 
senen Freistromdruck pœ Zu px = Pv/P» normiert. Für jedes Modell sind in 
Tabelle 3.1 der Mittelwert Py fa=0 und das 95 % Konfidenzintervall (CI) fiir 


Nexp. statistisch unabhängige Messungen ohne Rotation angegeben (Weidner 
et al., 2019b). 


Die Messungen ohne Rotation zeigen, dass Leitflächenhöhen h; kleiner als die 
mittlere Grenzschichtdicke ö = 5 mm keinen messbaren Einfluss auf den Basis- 
druck hatten. Leitflächen, die höher als die Grenzschicht waren, verminderten 
hingegen den Basisdruck. Mit zunehmender Höhe der Leitflächen stellten 
diese ein größeres Hindernis für die Strömung dar. Bei ansonsten unver- 
änderten Strömungsbedingungen führte dies im Vergleich zum längsange- 
strömten Zylinder ohne Leitflächen zu einem geringeren Basisdruck (Weidner 
et al., 2019b). Des Weiteren wurde für das Modell R6 mit nur drei Leitflächen 
im Vergleich zu den Modellen R4 und R5 mit sechs beziehungsweise zwölf 
Leitflächen ohne Modellrotation ein höherer Basisdruck gemessen. Da sich die 
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Konfidenzintervalle für die Modelle R4 und R5 überschneiden, konnte nicht 
abschließend geklärt werden, ob der Basisdruck ohne Modellrotation für das 
Modell RS tatsächlich geringer ist als für das Modell R4, wie die gemessenen 
Mittelwerte für die beiden Modelle vermuten lassen. 


Um den Einfluss der Modelldrehrate fy = Q/(27) mit der Modellwinkel- 
geschwindigkeit Q auf den zentralen Basisdruck untersuchen zu können, 
versetzte ein elektrischer Antriebsmotor die in Abbildung 3.8 gezeigten 
Modell-Hohlkörper auf der Länge Lp von 135 mm in Rotation. Die Modell- 
basis und das System zur Druckmessung befanden sich jedoch jederzeit in 
Ruhe. Die maximal erreichte Modelldrehrate fM,max für jedes Modell ist 
ebenfalls in Tabelle 3.1 angegeben und reicht bis zu einer entdimensionali- 
sierten Modelldrehrate Pmax = 27 fm.maxR/Ux ~ 0,19 (Weidner et al., 2019b). 


Für jede Drehrate fm wurde während des gleichen Windkanalexperimentes 
eine Referenzmessung ohne Rotation durchgeführt, sodass sich der Einfluss 
des Dralls auf den Basisdruck aus dem Verhältnis ¢ = Pý ,ot./Pù, far 
ergibt. Die Wiederholbarkeit der einzelnen Messungen, basierend auf den 
zufalligen Abweichungen in den Versuchsbedingungen, wurde durch die 


relative Standardabweichung der Messungen ohne Rotation zu s, De oh Pr fir=0 


bestimmt und ist ebenfalls in Tabelle 3.1 angegeben (Weidner et al., 2019b). 
Die statistische relative Standardabweichung der Einzelmessungen von <0,1 % 
aufgrund von Messabweichungen der verwendeten Messgeber (GE Sensing, 
2007) ist gegenüber der Wiederholbarkeit der Einzelmessungen in der Größen- 
ordnung von +2 % vernachlässigbar gering (Coleman und Steele, 2009). 


Abbildung 3.9 zeigt das Druckverhältnis ø in Abhängigkeit von der Modell- 
drehrate fm beziehungsweise p für die Modelle mit zwölf Leitflächen unter- 
schiedlicher Höhe, sowie für das Modell ohne Leitflächen. Basierend auf 
der Wiederholbarkeit der Einzelmessungen wurden die Unsicherheiten der 
gemessenen Druckänderungen mit der doppelten relativen Standardabwei- 


chung s in beide Richtungen angegeben (Weidner et al., 2019b). 


PS fra ol Po f= =0 


Für den längsangeströmten Zylinder ohne Leitflächen wurde keine Verän- 
derung des Basisdruckes in Abhängigkeit von der Modelldrehrate gemessen. 
Die durch die rotierende Oberfläche erzeugten azimutalen Schubspan- 
nungen waren nicht in der Lage, ausreichend Drehimpuls an die Strömung 
zu übertragen, was in einer vernachlässigbaren Strömungswinkelgeschwin- 
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Abbildung 3.9: Basisdruckverhältnis in Abhängigkeit von der Drehrate für Modelle mit unter- 
schiedlichen Leitflächenhöhen hyp (Weidner et al., 2019b). 
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digkeit wrot. im Nachlauf resultierte (Weidner et al., 2019b). Demgegenüber 
konnten die Modelle mit Leitflächen, aufgrund der höheren Rolldämpfungs- 
raten (Bhagwandin, 2016, DeSpirito und Heavey, 2004), einen größeren 
Drehimpuls in die Strömung einbringen. Die verwendeten Modelle mit 
Leitflächen erzeugten durch die azimutale Ablenkung der Strömung eine 
Strömungswinkelgeschwindigkeit wrot. im Nachlauf in der Größenordnung der 
Modelldrehrate Q (Weidner et al., 2019b). 


Die Messergebnisse in Abbildung 3.9 zeigen für die Modelle mit Leitflächen 
bei geringen Drehraten zunächst einen Druckanstieg von bis zu 3%, der 
mutmaßlich auf Inhomogenitäten der Windkanalströmung zurückzuführen ist 
(Hruschka und Leopold, 2015), da dieser Effekt bei numerischen Simulationen 
mit homogener Anströmung (Hruschka und Leopold, 2015) sowie bei Expe- 
rimenten in einem anderen Windkanal nicht auftraten (Weidner et al., 2019a). 
Für höhere Drehraten ist in Abbildung 3.9 ein Abfall des zentralen Basis- 
druckes mit ansteigender Modelldrehrate zu erkennen. Die kleinsten Leit- 
flächen mit einer Höhe von 1,5mm verminderten den Basisdruck um bis 
zu 10 % (Weidner et al., 2019b). Aufgrund des geringeren Rolldämpfungs- 
momentes im Vergleich zu Modellen mit höheren Leitflächen wurden mit den 
kleinsten Leitflächen die höchsten Drehraten erreicht (Weidner et al., 2019b). 
Im Gegensatz dazu wurden für die Leitflächen mit hr = 12,5mm unter 
Verwendung des Modells mit nyf = 12 nur geringe Drehraten erreicht, sodass 
die gemessenen Druckänderungen im Bereich der Messunsicherheit liegen 
(Weidner et al., 2019b). 


Tabelle 3.1 zeigt, dass eine Verminderung der Leitflächenanzahl nyp die 
maximal erreichbare Drehrate fM,max erhöhte. Die Änderungsrate des Basis- 
druckes mit der Drehrate war hingegen nahezu unabhängig von der Anzahl 
der Leitflächen, wie in den in Abbildung 3.10 dargestellten Messergebnissen 
für die drei Modelle mit App = 12,5mm zu erkennen ist. Das Modell mit 
nur = 3 zeigte bei fm = 300 Hz einen Druckabfall von etwa 5 %, was für 
kleinere Leitflächenhöhen erst bei fm > 500 Hz erreicht wurde. Mit zuneh- 
mender Leitflächenhöhe ergab sich somit eine Zunahme der Änderungsrate 
des Basisdruckes bezogen auf die Modelldrehrate. 
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Entdimensionalisierte Drehrate p 
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Abbildung 3.10: Basisdruckverhältnis in Abhängigkeit von der Drehrate für Modelle mit unter- 
schiedlicher Leitflächenanzahl nyp (Weidner et al., 2019b). 
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3.4 Druckabschätzung für rotierende 
Nachlaufströmungen 


Die gemessenen Druckänderungen aufgrund des eingebrachten Dralls wurden 
näherungsweise mittels der einfachsten Modellbeschreibung für rotierende 
Strömungen, einer Starrkörperrotation mit konstanter Winkelgeschwindigkeit 
Wrot., abgeschätzt (Weidner et al., 2019b). Hruschka und Leopold (2015) 
zeigten mittels numerischer Simulationen, dass die axialen und radialen 
Strömungsgeschwindigkeiten im Nachlauf durch Drehraten, wie sie in den 
experimentellen Untersuchungen von Kapitel 3.3 erzeugt wurden, kaum beein- 
flusst werden. In einem mit dem Fluidelement dV mitrotierenden Koordina- 
tensystem muss sich folglich ein Gleichgewicht 


ô 

S dV = pwr dV (3.3) 
zwischen der Kraft aufgrund eines radialen Druckgradienten und der durch 
den Drall erzeugten Zentrifugalkraft mit der radialen Koordinate r und der 
Dichte p ergeben, wie es auch in Gaszentrifugen entsteht (Kemp, 2009). Für 
kompressible Strömungen beschreibt die Isentropenbeziehung 


i 


poe (2y (Anderson, 2001) (3.4) 
Pb 


C 


mittels des Isentropenkoeffizienten y und den am Radius rp. herrschenden 
Randbedingungen pp, und pp, den Zusammenhang zwischen der Dichte p und 
dem Druck p. Somit konnte Gleichung (3.3) analytisch gelöst werden, um die 
Druckverteilung 


P(T, wrot.) = (3.5) 


ya 1 2 2 
Eana = 
be y C Poc rot. 


in Abhängigkeit von r und wrot, zu erhalten. Gleichung (3.5) beschreibt 
die Druckverteilung an der Basis, wenn die Randbedingungen ppc und ppc 
den Werten an der äußeren Grenze des Rezirkulationsgebietes bei ru. = R 
entsprechen (Weidner et al., 2019b). 
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3.4 Druckabschätzung für rotierende Nachlaufstromungen 


Hruschka und Leopold (2015) hatten mittels numerischer Simulationen 
gezeigt, dass die Länge des Rezirkulationsgebietes und damit der mittlere 
Ablenkungswinkel der Scherschicht für die in Kapitel 3.3 experimentell 
untersuchten Drehraten nahezu konstant sind. Deshalb wurde das Verhältnis 
zwischen dem Druck an der inneren Grenze des Rezirkulationsgebietes 
Poc und dem Druck unmittelbar vor der Hinterkante des Körpers pe als 
unabhängig von der Drehrate angenommen (Weidner et al., 2019b). Im 
Folgenden wurden hierzu die in Tabelle 3.1 angegebenen Druckverhältnisse 


Po. fa=0 bei Nulldrehzahl verwendet. Die sich ergebende Druckrandbedingung 


Poe = Ph, g,=0 ` Pe (3.6) 


führte mittels thermischer Zustandsgleichung eines idealen Gases (Oswatitsch, 
1952) und der Ruhetemperatur 79 = 295K, die im Rezirkulationsgebiet 
näherungsweise herrscht (Lamb und Oberkampf, 1995), zu der zusätzlich 
benötigten Dichterandbedingung 


Poc 
Be (3.7) 
Pog RiuttTo 


Für das Modell ohne Leitflächen wurde pe ~ pœ angenommen, da die 
ausschließlich durch azimutale Schubspannungen erzeugte rotierende Schicht 
stromauf der Ablösekante vernachlässigbar dünn war (Weidner et al., 2019b). 
Für die Konfigurationen mit Leitflächen wurde pe jedoch abhängig von 
Leitflächenhöhe und Drehrate definiert, da die für diese Modelle auf- 
tretenden Zentrifugalkräfte den Druck bereits vor der Ablösung in einer 
nichtvernachlässigbaren Größenordnung veränderten (Weidner et al., 2019b). 
Hierzu wurde angenommen, dass sich die Druckverteilung im Bereich der 
Leitflächen ebenfalls mittels Starrkörperrotation beschreiben lässt. Dies resul- 
tiert in einem reduzierten Druck vor der Ablösekante von 


Pr. 
y-l 


Po 


(3.8) 


1 
v1 ST 1 \y 2Rhır + h? 
atmen [ol "a(n 


Poo = see rot. 2 


im Vergleich zu dem Druck ohne Modellrotation von pe(hrr, wrot. = 0) = Po 
(Weidner et al., 2019b). 
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Die resultierenden Abschätzungen für das zentrale Basisdruckverhältnis & 
als Funktion von Leitflächenhöhe und Drehrate sind in den Abbildungen 3.9 
und 3.10 enthalten. Im Gegensatz zu den Messergebnissen zeigt die Abschätz- 
ung für das Modell ohne Leitflächen einen Druckabfall von bis zu 5 % für die 
untersuchten Drehraten. Die Abweichung zwischen der Modellbeschreibung 
und den Messergebnissen ist durch den geringen Drehimpuls bedingt, der durch 
das Modell ohne Leitflächen an die Strömung übertragen wurde. Daraus resul- 
tierte eine Winkelgeschwindigkeit der Strömung im Nachlauf, die wesentlich 
geringer war als die Modelldrehrate (wrot. « Q). Die stark vereinfachte Modell- 
beschreibung zur Abschätzung des Einflusses der Rotation auf den Basisdruck 
ist somit für das Modell ohne Leitflächen nicht anwendbar (Weidner et al., 
2019b). 


Für Modelle mit Leitflächen ist die Modellbeschreibung in der Lage, die 
Größenordnung des Druckabfalls, der mit der Leitflächenhöhe zunimmt, sowie 
dessen näherungsweise quadratischen Anstieg mit der Modelldrehrate wieder- 
zugeben. Für die höchsten Drehraten unterschätzt die Modellbeschreibung 
den Abfall des zentralen Basisdruckes deutlich, da sich die Länge des Rezir- 
kulationsgebietes, wie von Hruschka und Leopold (2015) gezeigt, aufgrund 
des höheren Dralls nichtvernachlässigbar verändert haben könnte, sodass die 
Gleichungen (3.3) und (3.6) nicht mehr gültig sind (Weidner et al., 2019b). 
Außerdem stellt die Annahme einer Starrkörperrotation des Rezirkulations- 
gebietes eine zu starke Vereinfachung dar. Unter der Annahme eines freien 
Wirbels mit einem inneren Strömungskern, dessen Bewegung sich wiederum 
als Starrkörperrotation beschreiben lässt, ergeben sich systematisch größere 
Druckänderungen (Weidner et al., 2019a). Aufgrund der begrenzten Daten- 
grundlage und zu Gunsten einer möglichst einfachen Modellbeschreibung 
wurde in dieser Arbeit jedoch auf die Erweiterung des Modells um diese 
zusätzliche Annahme verzichtet. 


Der Einfluss der Kompressibilität, das heißt die Verwendung von Glei- 
chung (3.4) anstelle der Annahme einer konstanten Dichte, auf die Druck- 
änderung war für die untersuchten Drehraten und Modelle mit <0,12 % 
vernachlässigbar klein (Weidner et al., 2019b). 
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3.5 Vergleich von rotierenden, geraden mit 
nichtrotierenden, angestellten Leitflächen 


Ux(r) A Leitfliche Ux(r) un | 
Qr D Qr 
7 Ux 
Un 
(a) Anströmung rotierender Leitflächen (b) Anströmung angestellter Leitflächen 


Abbildung 3.11: Vergleich der Anströmung von rotierenden, geraden Leitflächen und nicht- 
rotierenden, angestellten Leitflächen. 


Die rotierenden Leitflächen können, wie in Abbildung 3.1la skizziert, 
näherungsweise als quer zur Anströmung U,(r) bewegte ebene Platten 
betrachtet werden. Aufgrund der Rotation ergibt sich in dem mit der Leitflä- 
chenvorderkante mitrotierenden Koordinatensystem der Anströmwinkel 


Arot. = arctan (3.9) 


Qr 
Ur)‘ 
Unter der Annahme einer vernachlässigbaren Abhängigkeit des Anström- 
winkels von der radialen Koordinate können nichtrotierende Leitflächen mit 
einem konstanten Anstellwinkel, wie in Abbildung 3.11b skizziert, als Ersatz- 
modell für rotierende Leitflächen verwendet werden. Der Anstellwinkel A 
wurde entsprechend dem Anströmwinkel derrotierenden Leitflächen auf halber 
Höhe der Leitflächen zu 


A= Avot (r = R + hyp/2) (3.10) 


gewählt. Auf diese Weise wurde jedem nichtrotierenden Modell mit konstantem 
Anstellwinkel eine Äquivalenzdrehrate 


U(r =R+ hır/?2) tan A 
a re ee 11 
Piqui. 2n(R + hır/2) Bah 


zugeordnet. 
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Für ein Modell mit der gleichen Leitflächenanzahl nır und der gleichen Leit- 
flächenhöhe hyp wie Modell R6 aus Tabelle 3.1, jedoch mit einem Anstell- 
winkel der Leitflächen von A = 4°, ergab sich mit U(r = R+hır/2) ~ Us eine 
Äquivalenzdrehrate von fiqui. ~ 220 Hz. Der gemessene zentrale Basisdruck 
für dieses Modell relativ zum Modell R6 ohne Rotation und ohne Anstell- 
winkel A liegt, wie in Abbildung 3.10 gezeigt, im Bereich des mit fu = 220 Hz 
rotierenden Modells R6 mit geraden Leitflächen. 


Auch ein Modell mit A = 8° und der gleichen Leitflächenanzahl und -höhe wie 
Modell R2 aus Tabelle 3.1 zeigt in Abbildung 3.9 eine gute Übereinstimmung 
mit den Messungen der rotierenden Leitflächen. Da die mittlere Höhe der 
Leitflächen innerhalb der Grenzschicht lag, wurde U,(r = R + hır/2) entspre- 
chend der numerischen Simulation in Kapitel 3.1.2 zu etwa 440 m/s gewählt, 
sodass sich eine Äquivalenzdrehrate von Fiqui. ~ 460 Hz ergibt. 


Ein Vergleich für höhere Anstellwinkel A war aufgrund der experimentell 
limitierten Drehraten nicht möglich. Für den Vergleich höherer Anstellwinkel 
müsste zusätzlich der Unterschied in den Anströmgeschwindigkeiten von 
Ua rot. für die rotierenden Modelle und U, für die Modelle mit angestellten 
Leitflächen berücksichtigt werden, der nur für kleine Anstellwinkel A vernach- 
lässigbar ist. Experimentell war die Axialgeschwindigkeit U, durch die 
Geometrie des Windkanals und die sich ausbildende Grenzschicht festgelegt 
und konnte nicht verändert werden, sodass sich für rotierende Leitflächen eine 
relative Abweichung der Anströmgeschwindigkeit von 

U A,rot. 7 U x 1 


= -1 12 
Ux cosA (ata 


im Vergleich zu den angestellten Leitflächen ergibt. 


Außerdem wird durch rotierende Modelle die Grenzschicht bereits vor den 
Leitflächen teilweise in Rotation versetzt, sodass rotierende und angestellte 
Leitflächen nur für geringe Drehraten und Anstellwinkel direkt vergleichbar 
sind. 
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3.6 Fazit zur Studie mit rotierenden 
Strömungskörpern 


Die Messergebnisse für rotierende Leitflächen in den Abbildungen 3.9 und 3.10 
zeigen einen Abfall des zentralen Basisdruckes in Abhängigkeit von Leit- 
flächenhöhe, -anzahl und Drehrate (Weidner et al., 2019b). Die Messungen 
bestätigen somit die numerischen Ergebnisse von Hruschka und Leopold 
(2015) für geringe Drehraten. Die Veränderung der Leitflächenanzahl hatte für 
die verwendeten Modelle keinen messbaren Einfluss auf die Änderungsrate 
des zentralen Basisdruckes gezeigt (Weidner et al., 2019b). Im Gegensatz dazu 
nahm die Änderungsrate des zentralen Basisdruckes mit der Leitflächenhöhe 
zu (Weidner et al., 2019b). 


Der Einfluss von Drehrate und Leitflächenhöhe wurde mit einer stark ver- 
einfachten Modellbeschreibung (Weidner et al., 2019b) für die experimentell 
untersuchten Drehraten abgeschätzt. Für die höchsten experimentell unter- 
suchten Drehraten unterschätzt die Modellbeschreibung (Weidner et al., 
2019b) den Druckabfall, da die Annahme einer von der Rotation unabhängigen 
Expansion an der Modellhinterkante für höhere Drehraten nicht mehr gültig ist 
(Hruschka und Leopold, 2015) und die Annahme einer Starrkörperrotation des 
Rezirkulationsgebietes im Allgemeinen eine zu starke Vereinfachung darstellt 
(Weidner et al., 2019a). 


Es wurde außerdem experimentell gezeigt, dass für geringe Drallraten rotie- 
rende, gerade und nichtrotierende, angestellte Leitflächen einen ähnlichen 
Einfluss auf den zentralen Basisdruck haben. Da Modelle mit angestellten 
Leitflächen keine technischen Limitierungen der Drallraten durch den verwen- 
deten Antriebsmotor aufweisen, wurden diese im Folgenden zur Untersuchung 
wesentlich höherer Drallraten genutzt. Für Modelle mit hohen Leitflächen- 
anstellwinkeln ist ein direkter Vergleich mit rotierenden Leitflächen jedoch 
nicht möglich, sodass im Rahmen dieser Arbeit auf die Bestimmung einer 
Äquivalenzdrehzahl verzichtet wurde. 
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Mit den in Kapitel 3 gezeigten rotierenden Strömungskörpern konnte keine 
strukturelle Veränderung der Nachlaufströmung beobachtet werden, da die 
erzeugten Drallraten technisch begrenzt waren. Um höhere Drallraten unter- 
suchen zu können, wurden im Folgenden nichtrotierende Strömungskörper 
mit angestellten Leitflächen als alternative Methode zur Drallerzeugung im 
Nachlauf verwendet. Der in Kapitel 3.1 beschriebene Versuchsaufbau und die 
Anströmbedingungen blieben hierfür unverändert. 


4.1 Konfigurationen mit angestellten Leitflachen 


Abbildung 4.1: Skizze der verwendeten Strömungskörper mit zwölf angestellten Leitflächen 
(Weidner et al., 2017). 
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Tabelle 4.1: Übersicht über die verwendeten nichtrotierenden Strömungskörper. 


Leitflächen (LF) 
Modell Anzahl Höhe hir Anstellwinkel A 


Zylinder - - - 

A=0° 12 12,5 mm 0° 
A= 16° 12 12,5 mm 16° 
A = 32° 12 12,5 mm 32° 


Die verwendeten nichtrotierenden Strömungskörper hatten, wie in Abbil- 
dung 4.1 skizziert, zwölf Leitflächen mit einer Höhe App von 12,5 mm. 
Wie bereits bei den rotierenden Strömungskörpern reichten die Leitflächen 
von der Modellhinterkante 20mm stromauf und hatten eine Dicke von 
1mm. In Tabelle 4.1 sind die im Folgenden verwendeten Modelle, die sich 
ausschließlich im Leitflächenanstellwinkel A unterschieden, sowie der Zylinder 
ohne Leitflächen als Referenzmodell aufgelistet. 


4.2 Messmethoden 


Um eine Veränderung der Nachlaufstruktur nachweisen und beschreiben zu 
können, wurden die verwendeten Messmethoden im Vergleich zu Kapitel 3 
erweitert. In den folgenden Kapiteln 4.2.1 bis 4.2.3 werden die Methoden zur 
qualitativen Bestimmung des wandnahen Strömungsfeldes, der axialen und 
radialen Strömungsgeschwindigkeiten im Nachlauf sowie zur Druckmessung 
an der Oberfläche und an verschiedenen radialen Positionen der Basis der 
Strömungskörper beschrieben. 


4.2.1 Ölanstrichverfahren 


Mittels eines Ölanstrichverfahrens, wie es auch von Leopold (1993) verwendet 
worden ist, wurde die Orientierung des wandnahen Strömungsfeldes sichtbar 
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gemacht. Hierzu wurden Mischungen aus verschiedenfarbigen Pulvern und 
Mineralöl auf die Strömungskörper aufgetragen. Die Strömung über die 
Modelloberfläche erzeugte Schubspannungen, die das Öl auf der Oberfläche 
in Bewegung versetzten. Auf diese Weise entstanden charakteristische Strö- 
mungsmuster auf der Oberfläche. Eine Kamera zeichnete das transiente 
Verhalten der induzierten Ölströmung mit einer Bildwiederholungsrate von 
24 Hz auf. 


Für einen Strömungskörper ohne Leitflächen wurde das in Abbildung 4.2 
gezeigte Strömungsmuster aufgezeichnet, welches den Aufnahmen von 
Leopold (1993) ähnlich ist. Statt eines rotationssymmetrischen Strömungs- 
feldes im Nachlauf des Zylinders mit einem einzelnen Stagnationspunktes im 
Zentrum der Basis ergaben sich, aufgrund der geringen Länge der Lavaldüse 
und deren rechteckigem Querschnitt, zwei Stagnationspunkte bei z/R = +0,6 
und ein periodisches Strömungsverhalten, bei dem die Strömung zwischen 


Abbildung 4.2: Momentaufnahme der Visualisierung der Orientierung des wandnahen Strö- 
mungsfeldes mit zwei Stagnationspunkten (x) für einen längsangeströmten Zylinder, sichtbar 
gemacht mittels Ölanstrichverfahren. 
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den beiden Stagnationspunkten (x) alternierte (Leopold, 1993). Wie jedoch in 
Kapitel 5.1 gezeigt wird, wirkte sich diese Strömungsasymmetrie nur gering- 
fügig auf den Nachlauf der Modelle mit Leitflächen aus und beeinflusste somit 
die Untersuchungen im Rahmen dieser Arbeit nicht. 


Mit dem Ölanstrichverfahren wurde für die Modelle mit angestellten Leit- 
flächen qualitativ überprüft, ob eine rotierende Strömungsbewegung in das 
nahe Nachlaufgebiet eingebracht werden konnte. Außerdem dienten die Visua- 
lisierungen der Orientierung des wandnahen Strömungsfeldes dazu, eine 
mögliche Umkehr der radialen Strömungsrichtung in der Nähe der Basis 
qualitativ nachzuweisen, wie sie von Hruschka und Leopold (2015) mittels 
numerischer Simulationen gezeigt worden war. 


4.2.2 Particle Image Velocimetry (PIV) 


Die Axial- und Radialgeschwindigkeiten der Strömung im Nachlauf wurden 
mittels Particle-Image-Velocimetry-Methode (PIV) bestimmt. Da Strömungen 
durch Messungen mittels PIV-Methode im Allgemeinen kaum beeinflusst 
werden (Tropea et al., 2007), ist die Methode besonders geeignet, um Nachlauf- 
strömungen mit hoher räumlicher Auflösung zu untersuchen (Scarano und van 
Oudheusden, 2003). 


Erzeugung und Einbringen der Partikel 


Um das gesamte Strömungsfeld mit Partikeln versehen zu können, wurden 
die für die Messmethode notwendigen Partikel an zwei verschiedenen Stellen 
mit jeweils unterschiedlichen Methoden in die Strömung eingebracht. Die 
äußere Strömung wurde vor der Ruhekammer mit einem Dampf aus Wasser, 
Propylen-Glykol und weiteren Glykolen versehen. Durch das Einbringen der 
Partikel vor der Ruhekammer wurde eine homogene Verteilung der Partikel in 
der freien Anströmung erreicht, ohne die Nachlaufströmung durch den einge- 
brachten Massenstrom zu beeinflussen (Anhang B). Das nahe Nachlaufgebiet 
wurde durch eine Öffnung mit einem Durchmesser von 2mm im Zentrum 
der Modellbasis mit Partikeln versetzt. Ein Schlauch wurde an der Öffnung 
befestigt und durch die Modellhalterung aus dem Windkanal hinausgeführt. 
Zur Partikelerzeugung wurden am Ende des Schlauches Zigaretten angebracht. 
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Die Differenz zwischen dem Umgebungsdruck und dem zentralen Basisdruck 
führte zu dem in Anhang B.1 gemessenen Massenstrom in das Rezirkula- 
tionsgebiet von weniger als 5,0-10->kg/s. Der Einfluss des verursachten 
Massenstroms auf den Druck im Rezirkulationsgebiet konnte mit der semi- 
empirischen Korrelation von Danberg (1990) zu weniger als 0,1 % bestimmt 
werden (Anhang B.2). Numerische Simulationsergebnisse in Anhang B.2 
zeigen außerdem, dass der eingebrachte Massenstrom lediglich einen lokalen 
Einfluss auf die Struktur der Nachlaufströmung hatte. 


Folgeverhalten der Partikel 


Da bei einer PIV-Methode nicht unmittelbar die Geschwindigkeit des Fluids, 
sondern die Geschwindigkeit der in die Strömung eingebrachten Partikel 
vermessen wird (Tropea et al., 2007), ist die Genauigkeit der Messung unter 
anderem davon abhängig wie gut die Partikel der Strömung folgen können 
(Samimy und Lele, 1991). Die in Gasströmungen üblicherweise verwen- 
deten Partikel mit Durchmessern von 1-3 um haben im Allgemeinen eine 
höhere Trägheit als das Strömungsmedium, sodass die Partikel verzögert auf 
Geschwindigkeitsgradienten reagieren (Tropea et al., 2007). Insbesondere 
Diskontinuitäten im Geschwindigkeitsfeld, wie sie durch Verdichtungsstöße 
hervorgerufen werden, können mit einer PIV-Methode nur ungenau vermessen 
werden (Scarano und van Oudheusden, 2003). Eine exakte Vermessung 
der Strömungsgeschwindigkeit im Bereich der stumpfen Vorderkanten der 
verwendeten Leitflächen war deshalb nicht möglich. Im Nachlauf waren außer 
dem Rekompressionsstoß, der gegenüber den abgelösten Verdichtungsstößen 
vor den Leitflächen vergleichsweise schwach ist, keine Verdichtungsstöße zu 
erwarten. Scarano und van Oudheusden (2003) hatten gezeigt, dass Partikel 
mit einem Durchmesser von 0,27 um, die Agglomerationen mit einem effek- 
tiven Durchmesser von 0,5 um und einer Dichte von 10° kg/m? bilden, zur 
Vermessung einer Überschallnachlaufströmung geeignet sind. Am kritischsten 
Punkt des Strömungsfeldes — in der Scherschicht unmittelbar nach der 
Modellhinterkante (Herrin und Dutton, 1995) — wurde die Stokes-Zahl, als 
Maß dafür wie gut die Partikel der Strömung folgen können (Tropea et al., 
2007), von Scarano und van Oudheusden (2003) zu 0,14 bestimmt, basierend 
auf der nach Samimy und Lele (1991) berechneten Zeitskala der Strömung. 
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Der Zigarettenrauch im nahen Nachlauf enthielt feste und flüssige Partikel 
mit Durchmessern von 0,1 —0,7 um, deren geometrisches Mittel etwa 0,3 um 
betrug (Liund Hopke, 1993). Die Partikel bestanden hauptsächlich aus Wasser, 
Glycerin, Propylen-Glykol und Nikotin (Rodgman und Perfetti, 2013), die eine 
mittlere Dichte von 10° kg/m? aufwiesen. Mit der dynamischen Viskosität im 
Freistrom Ho, Änderungen der Strömungsgeschwindigkeiten in der Größen- 
ordnung von U% über eine Länge von 1 mm (Scarano und van Oudheusden, 
2003) und einem Partikeldurchmesser von weniger als 0,7 um ergibt sich für 
die Partikel des Zigarettenrauchs in der Scherschicht eine Stokes-Zahl von 
weniger als 0,12. Die Partikel des Zigarettenrauchs konnten der Strömung 
somit ähnlich gut folgen wie die Partikel, die in den Experimenten von Scarano 
und van Oudheusden (2003) eingesetzt worden waren. Nach Samimy und Lele 
(1991) ergibt sich deshalb in Bereichen hoher Geschwindigkeitsgradienten 
zwischen den gemessenen Partikelgeschwindigkeiten und der tatsächlichen 
Strömungsgeschwindigkeit eine systematische Abweichung mit einem quadra- 
tischen Mittel (RMS) von etwa 1%. Abgesehen von den Verdichtungsstößen 
waren die Geschwindigkeitsgradienten in der äußeren Strömung geringer als 
im nahen Nachlauf und der angrenzenden Scherschicht. Deshalb waren die 
dort verwendeten Partikel mit einem durchschnittlichen Partikeldurchmesser 
von etwa lum (Haertig et al., 1996) und einer mittleren Dichte von etwa 
10° kg/m? trotzdem in der Lage, der Strömung zu folgen. 


Eine Abweichung zwischen Strömungs- und Partikelgeschwindigkeit ergibt 
sich auch in rotierenden Strömungsgebieten. Die Zentripetalkräfte, die nötig 
sind, um die radialen Trägheitskräfte auszugleichen und somit eine kreis- 
förmige Bewegung um eine Rotationsachse zu erzeugen, sind von der 
Dichte abhängig und deshalb im Allgemeinen für das Strömungsmedium 
und die eingebrachten Partikel unterschiedlich. Aufgrund der höheren radialen 
Trägheit weisen die Partikel eine radiale Geschwindigkeit relativ zur Strömung 
auf (Samimy und Lele, 1990). Die sich ergebende systematische Abwei- 
chung skaliert mit der Rotationsgeschwindigkeit und dem Dichteverhältnis 
zwischen dem Strömungsmedium und den Partikeln (Samimy und Lele, 1991). 
Außerdem wird die Partikelkonzentration in den Wirbelzentren vermindert 
(Samimy und Lele, 1990), sodass die Geschwindigkeiten dort nur ungenau 
vermessen werden können. Die Position des Wirbelzentrums wird durch die 
gegenüber der Strömung höhere Trägheit der Partikel jedoch nicht beeinflusst 
(Samimy und Lele, 1990), sodass das mittlere Strömungsfeld zwar quantitative 
Abweichungen aufweist, qualitativ jedoch kaum beeinflusst wird (Jäger, 1983, 
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Jäger und George, 1983). Insbesondere der qualitative Nachweis einer verän- 
derten Strömungsstruktur im Nachlauf wird von dieser Messabweichungen 
nicht beeinflusst. 


Beleuchtung der Partikel 


Die eingebrachten Partikel wurden mit einem Nd:YAG-Laser-Doppelpuls- 
System beleuchtet. Das System erzeugte mit einer Frequenz von 7Hz zwei 
Laserpulse mit einem zeitlichen Abstand von 1 us. Jeder Laserpuls hatte eine 
Dauer von 10ns, sodass instantane Partikelpositionen aufgenommen werden 
konnten. Die Strahlengänge der beiden Laserpulse wurden kombiniert und 
mittels einer Linsenkombination in einen Laserlichtschnitt umgeformt, der in 
die zentrale Ebene hinter dem Strömungskörper geleitet wurde. Der Laser- 
lichtschnitt hatte im Messbereich eine Dicke von etwa 500 um. Die Streuung 
des Lichts durch die Partikel wurde mit einer Kamera mit Charge-Coupled- 
Device-Sensor (CCD) aufgenommen, die für jeden Laserpuls ein separates 
Bild mit einer Auflösung von 2048 px x 2048 px aufzeichnete. 


Kreuzkorrelationsalgorithmus 


Mittels eines iterativen Mehrfenster-Kreuzkorrelations-Algorithmus (Haertig 
et al., 2008) wurden die radialen und axialen Stromungsgeschwindigkeiten 
für jedes Bildpaar ausgewertet. Hierzu wurde ein äquidistantes Auswertungs- 
gitter mit einer Kantenlänge von 16px und einer räumlichen Auflösung 
von 38 Messpunkten pro Basisdurchmesser D verwendet. Mit diesem 
Auswertungsgitter konnten turbulente Strukturen in der Größenordnung von 
0,1 R räumlich aufgelöst werden. Für jeden Gitterpunkt wurden zwei quadra- 
tische Auswertungsfenster mit einer Kantenlänge von 48 px und eines mit 
32 px verwendet. Größere Fenster enthalten mehr Partikel als kleinere, sodass 
das Signal-Rausch-Verhältnis für das Korrelationsmaximum im Allgemeinen 
höher ist (Tropea et al., 2007). Allerdings mitteln größere Fenster kleinskalige 
Geschwindigkeitsgradienten heraus, weshalb hochauflösende Messergebnisse 
nur mit den kleineren Fenstergrößen möglich sind. Als endgültige Strömungs- 
geschwindigkeit wurde der Wert gewählt, der mit dem kleinsten Fenster und 
einem Verhältnis des quadratischen Mittelwerts (RMS) zum Maximum der 
Korrelation von höchstens 0,2 bestimmt wurde. Aus der Genauigkeit mit der 
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die Position des Korrelationsmaximums bestimmt werden konnte (Haertig 
et al., 2008), ergab sich für jeden Messpunkt eines einzelnen Bildpaares die 
Abweichung der instantanen Stromungsgeschwindigkeit zu +7 m/s mit einem 
Konfidenzniveau von 95 %. 


Geschwindigkeitskomponente senkrecht zum Laserlichtschnitt 


Der mit den angestellten Leitflächen eingebrachte Drall resultierte in einer 
azimutalen Strömungsgeschwindigkeit senkrecht zur Ebene des Laserlicht- 
schnitts. Hohe Azimutalgeschwindigkeiten führen dazu, dass Partikel sich aus 
der Ebene des Laserlichtschnitts herausbewegen und somit nur von einem 
der beiden Laserpulse beleuchtet werden. Je höher die Azimutalgeschwin- 
digkeit ist, desto geringer wird somit das Signal-Rausch-Verhältnis der Kreuz- 
korrelation, bis schlussendlich kein eindeutiges Maximum mehr bestimmt 
werden kann (Tropea et al., 2007). Bei der verwendeten Dicke des Laser- 
lichtschnitts von 500 um und einer Verzögerung von | us zwischen den beiden 
Laserpulsen waren die Partikelpositionen der beiden Einzelbilder in jedem 
Falle unkorreliert, wenn die Strömungsgeschwindigkeit senkrecht zum Laser- 
lichtschnitt höher als 500 m/s war. Für die verwendeten Anstellwinkel der 
Leitflächen von A < 32° wurden im Nachlauf allerdings mittlere Azimutal- 
geschwindigkeiten von weniger als U~% sin A = 276 m/s erwartet, sodass ausrei- 
chend Partikel von beiden Laserpulsen beleuchtet werden konnten, um ein 
eindeutiges Korrelationsmaximum bestimmen zu können. 


Mittelungsverfahren 


Während jedem Windkanal-Experiment wurden 100 Bildpaare aufgezeichnet. 
Das Ergebnis für jedes Bildpaar wurde aus dem kartesischen Koordinaten- 
system (x,y) in das zylindrische Koordinatensystem (x,r = |y|) überführt, 
sodass aus jeweils vier Windkanalexperimenten, mitinsgesamt 400 Bildpaaren, 
ein Mittelwert von bis zu 800 Einzelmessungen für jeden Messpunkt des zylin- 
drischen Bezugssystems gebildet werden konnte. Die Messabweichungen, die 
sich aus der Unsicherheit in der Bestimmung des Korrelationsmaximums 
ergaben, wurden aufgrund des turbulenten Charakters der Strömung als 
zufällige Abweichungen betrachtet, sodass deren Einfluss auf den Mittelwert 
des Strömungsfeldes aufgrund der hohen Anzahl an Einzelmessungen vernach- 
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lässigbar wurde. Die Messpunkte, für die die Standardabweichung des 
Mittelwertes der Strömungsgeschwindigkeiten bezogen auf Uœ% geringer 
als 0,5% war, wurden zur Berechung des zeitlich gemittelten Strömungs- 
feldes verwendet. Die Standardabweichung des Mittelwertes wurde von drei 
Größen beeinflusst: den physikalischen Schwankung aufgrund des transienten 
Strömungsverhaltens, den zufälligen Messabweichungen und der Anzahl der 
erfolgreichen Einzelmessungen für jeden Messpunkt. Eine hohe Standardab- 
weichung des Mittelwertes wurde somit in den unstetigen Strömungsgebieten 
sowie in Bereichen mit geringer Partikelanzahl oder geringer Beleuchtungs- 
stärke erwartet. Außerdem reduzierte sich die Anzahl der Einzelmessungen 
aufgrund von Strömungsgeschwindigkeiten senkrecht zum Laserlichtschnitt, 
wie es für die Modelle mit angestellten Leitflächen erwartet wurde. Dies führte 
für einzelne Messpunkte dazu, dass die Strömungsgeschwindigkeiten nur für 
einen geringen Anteil der 800 Einzelmessungen bestimmt werden konnten. Der 
Mittelwert der Messergebnisse konnte dann von dem tatsächlichen zeitlichen 
Mittel des Strömungsfeldes abweichen. Um eine auf diese Weise entstehende 
systematische Abweichung der Messergebnisse vom wahren Mittelwert der 
Messgröße zu minimieren, wurden für das gemittelte Strömungsfeld nur 
Messpunkte verwendet für die mindestens 200 Einzelmessungen vorlagen. 
Die Messunsicherheit der PIV-Ergebnisse mit einem Konfidenzniveau von 
95% als Kombination der zufälligen Abweichungen und der systematischen 
Abweichung aufgrund der verzögerten Beschleunigung der Partikel ergab sich 
dann zu weniger als +8 m/s. 


Ergebnisse für einen längsangeströmten Zylinder 


Für einen längsangeströmten Kreiszylinder konnten auf diese Weise die in 
Abbildung 4.3 gezeigten Felder der zeitlich gemittelten, relativen Axial- und 
Radialgeschwindigkeiten 

Us pry = Ux/U%, (4.1) 


X 


UF piv = U, /Ux (4.2) 


y 


bestimmt werden. Die Position xsp des hinteren Stagnationspunktes ergab sich 
aus dem Schnittpunkt der Isotache U, = 0 mit der Achse bei x/R ~ 2,84 + 0,05. 
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(a) Relative Axialgeschwindigkei 
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(b) Relative Radialgeschwindigkeit 


Abbildung 4.3: Zeitlich gemitteltes, relatives Geschwindigkeitsfeld im Nachlauf eines längsan- 
geströmten Zylinders (PIV-Messung), eingezeichnet sind die zeitlich gemittelten Stromlinien, der 
Expansionsfächer an der Modellhinterkante (---), die Rekompression (— - —), der Verdich- 
tungsstoß am Ende der Lavaldiise (------ ) und die Isotache Ux = 0. 
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Tabelle 4.3: Vergleich der Position des hinteren Stagnationspunktes, der maximalen relativen 
Rückströmung und deren absoluter und relativer, axialer Position für experimentell vermessene 
Überschall-Zylindernachlaufströmungen mit Anström-Machzahlen von 2,0 und 2,46. 


Moo Xsp/R min UF X min ux/R X min us /Xsp 
Leopold (1993) 2,0 2,94 +0,06 -0,28 1,54 + 0,48 = 
Herrin und Dutton 
(1994) 2,46 2,65 -0,27 1,5 57 Io 
PIV Messung 2,0 2,84 +0,05 -0,26 + 0,02 1,67 +0,05 59 +2% 


Tabelle 4.3 zeigt, dass die relative Länge des Rezirkulationsgebietes xsp/R 
mit dem von Leopold (1993) mittels LDA bestimmten Messwert gut überein- 
stimmt. Die Ergebnisse von Herrin und Dutton (1994) zeigen hingegen 
ein kürzeres Rezirkulationsgebiet, da die Experimente bei einer höheren 
Machzahl durchgeführt worden waren (Murthy und Osborn, 1976). Die 
von Herrin und Dutton (1994) gemessene relative Position der maximalen 
Rückströmgeschwindigkeit x|minu*/xsp stimmt jedoch mit dem Ergebnis 
der PIV-Messungen gut überein. Außerdem zeigen die PIV-Messungen eine 
maximale relative Rückströmgeschwindigkeit min UX von —0,26 + 0,02 an der 
Achse bei x|minu+/R = 1,67 +0,05, die im Mittel um 7 % beziehungsweise 
4% geringer ist als bei Leopold (1993) und Herrin und Dutton (1994). 


4.2.3 Druckmessungen 
Drucksensitive Farbe (PSP-Methode) 


Das Druckfeld im Bereich der Leitflächen wurde mittels drucksensitiver 
Farbe (pressure-sensitive paint, PSP) bestimmt. Diese Methode ist in der 
Lage, den statischen Druck an der Oberfläche eines Körpers mit einer hohen 
räumlichen Auflösung zu bestimmen. Sie wurde verwendet, um Verände- 
rungen im Druckfeld vor der Ablösung der Strömung an der Hinterkante des 
Körpers qualitativ zu vermessen. 
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Für die in dieser Arbeit verwendete PSP-Methode wurden Porphyrin-Moleküle 
(PtTFPP) innerhalb einer sauerstoffdurchlässigen Matrix auf die Modell- 
oberfläche aufgetragen. PtTFPP-Molekiile fluoreszieren, wenn sie mit Licht 
der Wellenlängen 380-450nm beleuchtet werden. Für die Anregung der 
Moleküle wurde Licht der Wellenlänge (400 + 15)nm verwendet. Dies führte 
zu einer Transition der Elektronen in einen höheren Energiezustand. Angeregte 
Elektronen der PtTFPP-Molekiile können in ihren ursprünglichen Energiezu- 
stand zurückfallen, indem sie entweder Photonen mit 620- 750nm Wellen- 
länge aussenden oder indem sie die Energie an Sauerstoffmoleküle in der 
Umgebung übertragen. Diese konkurrierenden Prozesse bestimmen die Inten- 
sität des emittierten Lichts. Bei geringeren Drücken führt die geringere Sauer- 
stoffkonzentration zu einer erhöhten Fluoreszenz (Tropea et al., 2007). Das 
emittierte Licht wurde mittels einer monochromen Kamera aufgezeichnet. 
Vor der Kamera wurde ein Bandpassfilter mit einem Wellenlängenbereich 
von (650 + 10)nm angebracht, um das intensive Licht herauszufiltern, das 
zur Anregung der Moleküle benötigt wurde. Der an der Modelloberfläche 
herrschende Druck p wurde mittels der Stern-Volmer-Relation (Stern und 
Volmer, 1919) 


re = A+B (4.3) 


Pref. 


mit den temperaturabhängigen Koeffizienten A und B berechnet. Für die 
verwendete Farbe wurden die Koeffizienten für Temperaturen von (283 + 10)K 
sowie 0,1 - 10° Pa < p<s?2: 10° Pa zu A = 0,14 +0,01 und B = 0,86 + 0,04 
bestimmt. Mit Hilfe einer Referenzintensität Zef, unmittelbar nach jedem 
Windkanalversuch für einen Referenzdruck pref., der dem Umgebungsdruck 
entspricht, wurde aus der aufgezeichneten Intensität J des fluoreszierenden 
Lichts während des Versuches der lokal herrschende Druck p berechnet 
(Martinez, 2007). 


Die hohen Druckgradienten einer Überschallströmung, wie sie durch Verdich- 
tungsstöße hervorgerufen werden, sorgen für große Intensitätsunterschiede des 
emittierten Lichts. Allerdings verhinderten die parallel auftretenden Tempe- 
raturgradienten in Verbindung mit der Temperaturabhängigkeit der Koeffizi- 
enten A und B von 0,4 %/K ein exaktes Messergebnis (Tropea et al., 2007). 
Außerdem unterlagen die PtTFPP-Molekiile einem Alterungsprozess, da sie 
durch Feuchtigkeitsabsorption in einen nichtanregbaren Zustand überführt 
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wurden. Deshalb nahm die Intensität des emittierten Lichts mit der Zeit ab, 
sodass zusätzlich eine zeitabhängige Abweichung vom tatsächlichen Wert der 
Messgröße entstand (Tropea et al., 2007). 


Abbildung 4.4 zeigt das Ergebnis der Druckmessung mittels PSP-Methode 
an der Modelloberfläche des längsangeströmten Zylinders mit einem zeitlich 
gemittelten Druck von 0,66 - 10°Pa. Das Messergebnis der PSP-Methode 
weicht um etwa 8 % von dem mittels Druckgebern bei x/R = -1,925 gemes- 
senen mittleren statischen Druck von 0,61: 10°Pa ab, der dem theoretisch 
erwarteten Wert einer isentropen Expansion der Strömung mit einer Abwei- 
chung von lediglich 1,3 % entspricht (Kapitel 3.1.1). Obwohl die Ergebnisse 
der PSP-Methode quantitative Messunsicherheiten in der Größenordnung 
von 10% aufweisen, konnten die PSP-Ergebnisse für qualitative Vergleiche 
des statischen Druckes an der Modelloberfläche verwendet werden. Hierfür 
wurden die Messergebnisse der PSP-Methode im Folgenden mit dem mittleren 
statischen Druck vor den Leitflächen normiert. 


1,0 Pp [105 Pa] 
1,2 

0,5 
1,0 
< 0,0 0,8 

A 

0,6 

-0,5 
0,4 
-1,0 0,2 
0,0 


Abbildung 4.4: Druckmessung an der Modelloberfläche eines längsangeströmten Zylinders 
stromauf der Modellhinterkante mittels PSP-Methode, eingezeichnet sind die Positionen der 
Druckmessungen aus Kapitel 3.1.1 (x). 
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Druckmessung mittels Miniaturmessgebern 


Die Druckverteilung an der Modellbasis wurde mittels Druckgebern (Kulite, 
2014) und einem Messverstärker vermessen, die direkt in das Modell integriert 
wurden. Die Druckgeber wurden an der Basis an verschiedenen radialen 
Positionen angebracht. Ein Druckgeber wurde in die zentrale Bohrung 
integriert, die auch für das Einbringen der Partikel während der PIV-Messungen 
verwendet wurde. Konstruktionsbedingt wurden die weiteren Bohrungen zur 
Druckmessung über die gesamte Basisfläche verteilt, wie es in Abbildung 4.5 
dargestellt ist. In Kapitel 5.5 und Kapitel 7.1.5 wird gezeigt werden, dass der 
Basisdruck für die verwendeten Modellgeometrien und radialen Positionen 
zur Druckmessung nahezu unabhängig von der Azimutalposition war. 


Der Druck an der Basis pp wurde mit dem statischen Druck in der Anströmung 
Poo zu dem relativen Basisdruck px normiert, sodass der Druckanstieg in der 
Anströmung von bis zu 2% während eines Windkanalversuches von bis zu 
40s das Messergebnis nicht beeinflusste. Der während eines Experimentes 
ermittelte Messwert 


Ph Bi (Du) + ps Ppi +Ep> (4.4) 


Abbildung 4.5: Position der Druckbohrungen an der Modellbasis (Maße in mm). 
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ist eine Kombination aus dem wahren Wert der Messgröße (px Jessie den syste- 
matischen Abweichungen Ppt „i und den zufälligen Abweichungen Epr- Bei dem 


Wert p handelt es sich um die theoretisch messbare, für jede Konfiguration 
konstante Summe aus wahrem Wert der Messgröße und den systematischen 
Abweichungen Ppz.i- 


Die systematischen Abweichungen der Druckgeber Ppr „i konnten nicht exakt 
bestimmen werden, wurden jedoch mittels der Herstellerangaben als statis- 
tische Standardabweichungen bp, und b,,, approximiert (Coleman und Steele, 
2009). 


Die statistische Standardabweichung b,, ergibt sich aus den temperatur- 
bedingten Nullpunktverschiebungen bp, ro und Sensitivitätsveränderungen 
bp,,r1 sowie aufgrund von Nichtlinearitäten, Hysteresen und Zeitabhängig- 
keiten der Messgeber by, nery (Kulite, 2014) zu 


bp, = \ (bp,.70)” + (dm.rı) + (bm Nerv). (4.5) 


Die Messstrecke wurde bei Umgebungstemperatur kalibriert. Da wahrend der 
Versuche in der Messkammer des Windkanals tiefere Betriebstemperaturen 
vorlagen, entstand eine temperaturbedingte systematische Abweichung. Die 
maximale Temperaturdifferenz zwischen den Bedingungen während der Kali- 
brierung und den Versuchen wurde zu (25 + 5) K abgeschätzt, da die Tempe- 
ratur des Modells und damit der Messgeber tiber die Versuchszeit abfiel. 


Messgeber und Verstärker für den statischen Druck in der Anströmung 
waren außerhalb des Windkanals angebracht und unterlagen deshalb keinen 
Temperaturschwankungen. Der gemessene Wert po. weicht deshalb systema- 
tisch mit einer statistischen Standardabweichung 


bpo = bpo, NLTV (4.6) 


ausschließlich basierend auf Nichtlinearitäten, Hysteresen und Zeitabhängig- 
keiten der Messstrecke von dem wahren Wert der Messgröße (Poo) req, ab (GE 
Sensing, 2007). 
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Tabelle 4.5: Abschätzung der relativen systematischen Messabweichungen für po = 0,61 - 10° Pa 
und pp = 0,36 - 10° Pa. 


Dpy,To/ Pb 16 +03 % 
bp,.T1/Pb 03 +0,1 % 
bp,,NLTV [Pb 0,2 % 
bps,NLTV/Po <0,1 % 


In Tabelle 4.5 wurden die relativen Standardabweichungen der einzelnen syste- 
matischen Messabweichungen für pə = 0,61 - 10° Pa und einem angenom- 
menen Basisdruck von pp = 0,36 - 10° Pa abgeschätzt. 


Aufgrund der Fehlerfortpflanzung ergibt sich für den relativen Basisdruck p< 
eine relative systematische Standardabweichung 


2 2 2 2 
Pre Z (22) E (m) f (“a Pr (et) (47) 
Po Pb Pb Pb Poo 
als Kombination der relativen systematischen Standardabweichungen fiir pp 
und po (Coleman und Steele, 2009). Gegenüber der relativen Standard- 


abweichung b,,,ro/p» sind die anderen in Tabelle 4.5 angegebenen systemati- 
schen Messabweichungen vernachlassigbar klein. 


Fiir den theoretisch messbaren Wert pe ergab sich somit eine Wahrscheinlich- 
keitsdichteverteilung T basierend auf den von den Herstellern der Messgeber 


spezifizierten, E auftretenden systematischen Messabweichungen. Mit 


einer Wahrscheinlichkeit von 95% lag der theoretisch messbare Wert p< 


+ 196b p|: Die 


real ~~ 


innerhalb des in Abbildung 4.6a gezeigten Intervalls [o ) 
Summe der Abweichungen %; Bpx i war fiir alle Versuche mit dem gleichen 
Modell identisch und führte somit zu einem konstanten Wert De Dieser theore- 
tisch messbare Wert př wurde, wie in Abbildung 4.6b skizziert, von zufälligen 


Abweichungen Ep* überlagert. Die zufälligen Abweichungen Eps resultierten in 


einer r Normalverteilung der Messwerte mit der Wahrscheinlichkeitsdichte Y, * 


um pe mit der Standardabweichung s—. Mit der Anzahl der Versuche Ngxp. 
b pP p 
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Wahrscheinlichkeits- 
dichte Y 


* 
(pi la 


(a) Wahrscheinlichkeitsdichteverteilung von pe aufgrund systematischer Messabweichungen 


Wahrscheinlichkeits- 
dichte Y 


* 
D 
(pi Ves i b 


(b) Wahrscheinlichkeitsdichteverteilung von Pe aufgrund zufälliger Messabweichungen 


A 


PAM Exp.) 
Py 


Up» ‚95 To 


Wahrscheinlichkeitsdichte Y 


* 7k Dk 
(pi = Po Po 


(c) Wahrscheinlichkeitsdichteverteilungen von Pe und DE für NExp. < œ 


Abbildung 4.6: Schematische Darstellung systematischer und zufälliger Messabweichungen und 
deren Einfluss auf die Abweichung des Messwertes von dem wahren Wert der Messgröße. 
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wurde, wie in Abbildung 4.6c gezeigt, der Mittelwert pe der Einzelmessungen 
Pej zu 


zZ 1 Nexp. 
tai * 
a (4.8) 
Po Nexp. 2 Pb, 


und die Standardabweichung des Mittelwertes Sox zu 
b 


Nexp. 9 
> (rt, - 7%) (49) 
j=l 


1 
Nexp.(NExp. u 1) 
bestimmt. Die Standardabweichung des Mittelwertes Se ist ein MaB fiir die 

b 
Ubereinstimmung von p% und př. Für eine große Versuchsanzahl Ngxp, gilt 


lim s—=0, (4.10) 


NExp. > Py 


sodass die gemittelte zufallige Abweichung Ep gegenüber den systematischen 


Abweichungen bpx vernachlassigbar klein wird (Coleman und Steele, 2009). 


Die kombinierte statistische Abweichung aus den systematischen und zufäl- 
ligen Messabweichungen wurde ohne diese Approximation zu 


2 2 
Ups = (b) + (555) (4.11) 
bestimmt. Mit Hilfe der kombinierten statistischen Abweichung und der 
Student-t-Verteilung wurde die Messunsicherheit 


Ups 95% = t95 %(Vp*) ` Ups (4.12) 
berechnet, aus der sich das Konfidenzintervall |p% + Up» 95 go | um px ergab, 


welches (pk ) real 10 95 % der nach diesem Vorgehen durchgeführten Messungen 
enthält (JCGM 100:2008, 2010). Die effektive Anzahl der Freiheitsgrade Ypž 
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zur Bestimmung des Koeffizienten to5 4, wurde mit der Formel von Welch und 
Satterthwaite (JCGM 100:2008, 2010) zu 


2 
(er) 
PSSST (4.13) 
(sie) [Ne= 1) + X (bye a). 
angenähert. Die hierfür notwendige Anzahl der Freiheitsgrade vp; der syste- 
matischen Abweichungen wurde durch 


-2 
Ab,» ; 
er pa (4.14) 
i 2 box i 

*, 


approximiert (JCGM 100:2008, 2010). Für eine hohe Anzahl effektiver 
Freiheitsgrade, das heißt kleine Unsicherheiten Abp> i in den systematischen 
Abweichungen Dox i> gilt 


lim Upos% = 1,96ups. (4.15) 


Für den längsangeströmten Zylinder zeigen die Mittelwerte der relativen 
Basisdrücke p* als Schätzwerte der wahren Messgröße (px) „ den in Ab- 


bildung 4.7 gezeigten radialen Basisdruckverlauf. Die Fehlerbalken zeigen 
fiir jeden Messwert die Messunsicherheiten +U p 95 % aus Tabelle 4.6. Die 
Messergebnisse haben eine mittlere relative Abweichung von 3,5 % gegenüber 
den Messergebnissen von Leopold (1993), die mit dem gleichen Windkanal- 
aufbau für einen Zylinder mit einem Durchmesser von 38,66 mm bestimmt 
worden waren. Die Abweichungen zwischen den Messergebnissen von Leopold 
(1993) und den Messungen für den längsangeströmten Zylinder liegen somit 
im Bereich der Messunsicherheiten. Außerdem ergaben Messungen mit 
einem zusätzlichen, außerhalb des Windkanals montierten Messgeber (GE 
Sensing, 2007), der mit der zentralen Öffnung der Basis verbunden war, einen 
gemittelten relativen Druck von (0,610 + 0,002)im Zentrum der Basis. Die 
Messergebnisse des externen Messgebers weichen damit lediglich um 1,6 % 
von dem Messergebnis des zentral in das Modell integrierten Messgebers ab. 
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Abbildung 4.7: Gemessener, zeitlich gemittelter, radialer Druckverlauf am Heck eines längsan- 
geströmten Zylinders. 


Tabelle 4.6: Messabweichungen und resultierende Messunsicherheit für den Basisdruck des längs- 
angeströmten Zylinders. 


r/R 0,00 0,25 0,40 0,55 0,70 0,825 


Nexp. 25 18 22 3 3 22 
pe 0,62 0,64 0,65 0,61 0,66 0,64 


S— <0,01 <0,01 <0,01 <0,01 <0,01 <0,01 


bps 0,01 0,01 0,01 0,01 0,01 0,01 
Ups 0,01 0,01 0,01 0,01 0,01 0,01 
Ups os% 0,02 0,02 0,02 0,02 0,02 0,02 
U D. * 95% 
= 4% 3% 3% 4% 3% 4% 
Ph 
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5.1 Orientierung des wandnahen 
Strömungsfeldes 


Abbildung 5.1 zeigt die Orientierung der wandnahen Strömungsfelder für 
Strömungskörper mit unterschiedlichen Leitflächenanstellwinkeln A. Für A = 
0° ergab sich das in Abbildung 5.la gezeigte, nahezu rotationssymmetrische 
Strömungsbild. Im zeitlichen Mittel bildete sich im Zentrum der Basis ein 
Stagnationspunkt aus und das auf das Modell aufgebrachte Öl strömte zur 
äußeren Kante der Basis. Dieses Verhalten ist charakteristisch für eine turbu- 
lente achsensymmetrische Überschallnachlaufströmung (Herrin und Dutton, 
1994). Das periodische Verhalten der wandnahen Strömung an der Basis mit 
zwei Stagnationspunkten aufgrund der asymmetrischen Düsenströmung, wie 
es in Abbildung 4.2 für den längsangeströmten Zylinder ohne Leitflächen 
gezeigt wurde, konnte für das Modell mit Leitflächen nicht festgestellt werden. 


In Abbildung 5.2a wurde der vor den stumpfen Leitflächenvorderkanten entste- 
hende abgelöste Verdichtungsstoß skizziert, der in der äußeren Strömung in 
einen schrägen Verdichtungsstoß übergeht. Der abgelöste Stoß führt in der 
Grenzschicht zu einer Strömungsablösung, die sich aufgrund der wandnahen 
Unterschallströmung stromauf des abgelösten Stoßes befindet. Aufgrund der 
Strömungsablösung in der Grenzschicht bildet sich ein schräger Ablösungsstoß 
aus. Unmittelbar vor der Leitflächenvorderkante wird die abgelöste Strömung 
der Grenzschicht durch einen weiteren Verdichtungsstoß komprimiert und in 
eine Unterschallströmung überführt (Dolling und Bogdonoff, 1980). Durch 
die Strömungsablösung vor den Leitflächenvorderkanten entstehen somit die 
in Abbildung 5.2b skizzierten Hufeisenwirbel. Aufgrund der Hufeisenwirbel 
löst die Strömung zwischen den Leitflächen von der zylindrischen Modell- 
oberfläche ab, sodass sich die in Abbildung 5.1a erkennbaren Bereiche ergaben, 
in welche das stromauf aufgetragene Öl nicht eindringen konnte. 
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Ablösung spiralförmiges 
durch Ver- 


ie + d 7 


Ablösungs- * 


Muster 


| Befestigungs- Befestigungs- 
schrauben (x) schrauben (x) 


uinsförmige 
Ablösung 


Ablösungs- 


stoß u x 


Befestigungs- 
sChrauben (x) 


(c) A = 32° 


Abbildung 5.1: Olvisualisierung der Orientierung des wandnahen Strömungsfeldes fiir Strö- 
mungskörper mit unterschiedlichen Leitflächenanstellwinkeln A (Weidner et al., 2016, 2017). 
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Interaktion 
benachbarter 
Stöße 


abgelöster Stoß 


EIER 


Ablösungsstoß 


abgelöster_ 
Stoß 7 
Grenzschicht 


(a) Skizze der Verdichtungsstöße an den (b) Dreidimensionale Darstellung der Hufeisen- 
Leitflächenvorderkanten wirbel um die Leitflächen 


Abbildung 5.2: Skizzen zur Entstehung der Hufeisenwirbel an den Leitflächenvorderkanten. 


Abbildung 5.2b stellt außerdem dar, wie die von zwei benachbarten Leitflächen 
hervorgerufenen Verdichtungsstöße zwischen den beiden Leitflächen mitein- 
ander interagieren. Anschließend treffen die Verdichtungsstöße auf der 
Oberfläche der benachbarten Leitfläche auf. Dort erzeugen sie eine weitere 
Strömungsablösung, die ebenfalls in Abbildung 5.1a zu erkennen ist. 


In Abbildung 5.1b ist ein spiralförmiges Muster zu erkennen, welches sich 
an der Basis des Strömungskörpers mit A = 16° ausbildete. Die Azimutal- 
komponente der Geschwindigkeit ist eine Folge des durch die Leitflächen 
eingebrachten Dralls, der in das Strömungsgebiet an der Basis transportiert 
wurde. Außerdem zeigt die Ölansammlung im Zentrum der Basis, dass das 
Öl in Richtung der Modellachse transportiert wurde. Folglich wurde die 
radiale Strömungsrichtung in der Nähe der Basis gegenüber dem Modell 
mit A = 0° umgekehrt. Dies wurde durch die transienten Aufnahmen der 
Ölströmung bestätigt und deutet auf die in den Simulationen von Hruschka 
und Leopold (2015) gezeigte Veränderung der Strömungsstruktur aufgrund 
des eingebrachten Dralls hin. 


Für das Modell mit A = 32° wurde das Öl, wie in Abbildung 5.1c gezeigt, 
sowohl von der äußeren Kante der Basis als auch von deren Zentrum zu 
einer ringförmigen Olansammlung bei r/R ~ 0,8 transportiert. Dort löste die 
Strömung von der Modelloberfläche an der Basis ab. Im Gegensatz zu dem 
Modell mit A = 16° ist an der Basis keine Bewegung des Öls in azimutaler 
Richtung erkennbar. Die Orientierung des wandnahen Stromungsfeldes für das 
Modell mit A = 32° deutet somit auf eine weitere strukturelle Veränderung der 
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Strömung in der Nähe der Basis hin. Außerdem befindet sich die Strömungs- 
ablösung in Abbildung 5.1c aufgrund der Verdichtungsstöße an den Leitflä- 
chenvorderkanten gegenüber den Abbildungen 5.1a und 5.1b weiter stromauf. 
Bei den beobachteten Veränderungen des Strömungsfeldes für das Modell mit 
A = 32° handelt es sich vermutlich nicht um die Folge einer vollständigen 
Blockierung der Strömung zwischen den Leitflächen, wie die im Folgenden 
beschriebenen Messungen der Stromungsgeschwindigkeiten gezeigt haben. 


5.2 Geschwindigkeitsfeld im Nachlauf 


Die Abbildungen 5.3 und 5.4 zeigen die mittels PTV-Methode bestimmten, 
zeitlich gemittelten relativen Axial- und Radialgeschwindigkeiten in der 
zentralen Ebene hinter den Strömungskörpern mit den unterschiedlichen 
Leitflächenanstellwinkeln A. Die azimutale Lage Ayırpiv der PIV-Ebene 
relativ zur Position der Leitflächen an der Modellhinterkante wurde für die 
Modelle mit Leitflächenanstellwinkeln von 0°, 16° und 32° zu (18 +1)°, 
(27 + 2)° beziehungsweise (7 +3)° bestimmt. Die Axialgeschwindigkeiten 
unmittelbar hinter den Leitflachen (LF) war fiir alle Modelle in der gleichen 
Größenordnung wie für den längsangeströmten Zylinder in Abbildung 4.3. 
Eine vollständige Blockierung der Strömung zwischen den Leitflächen lag 
somit vermutlich nicht vor. Der Betrag der Axialgeschwindigkeiten unmit- 
telbar hinter den Leitflächen des Modells mit A = 32° wurde gegenüber dem 
Modell mit 2 = 0° jedoch aufgrund des höheren Anstellwinkels der Leitflächen 
um etwa 10 % verringert, wie in Abbildung 5.3 zu erkennen ist. 


Die in den Abbildungen 5.3a und 5.4a gezeigten Geschwindigkeitsfelder für 
A = 0° zeigen das typische Rezirkulationsgebiet einer turbulenten, achsen- 
symmetrischen Überschall-Nachlaufströmung. Im Vergleich zu den Ergeb- 
nissen für den längsangeströmten Zylinder in Abbildung 4.3 sind in der äußeren 
Strömung schräge Verdichtungsstöße zu erkennen, die von den Leitflächen- 
vorderkanten hervorgerufen wurden. Die axiale Position des hinteren Stagna- 
tionspunkts wurde für das Modell mit A = 0° zu x/R = 2,90 + 0,05 bestimmt. 
Im Vergleich zum längsangeströmten Zylinder ergab sich somit eine tenden- 
zielle Verlängerung des Rezirkulationsgebietes. Dies ist neben den durch 
die Leitflächen erzeugten in Strömungsrichtung orientieren Wirbelstruk- 
turen hauptsächlich eine Folge der geringeren Machzahl unmittelbar vor der 
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Ablösung aufgrund der Verdichtungsstöße an den Leitflächenvorderkanten. 
Die geringere Machzahl vor der Ablösung resultierte in einem erhöhten Basis- 
druck relativ zum Druck unmittelbar vor der Ablösung (Lamb und Oberkampf, 
1995). Deshalb wurde die Expansion an der Modellhinterkante im Vergleich 
zum längsangeströmten Zylinder abgeschwächt, sodass auch der Grad der 
Strömungsumlenkung in Richtung der Modellachse abnahm. 


Die axiale Position der maximalen Rückströmgeschwindigkeit wurde für das 
Modell mit 2 = 0° zu x/R = 1,41 + 0,05 bestimmt und lag damit im Mittel 
um 16% näher an der Basis als für den längsangeströmten Zylinder. Der 
Betrag der maximalen relativen Rückströmgeschwindigkeit ist mit 0,27 + 0,02 
im Vergleich zum längsangeströmten Zylinder um 4 % höher. Dies ist eine Folge 
der Hufeisenwirbel, die an den Leitflächenvorderkanten erzeugt wurden und 
in der Scherschicht zu einer erhöhten Durchmischung zwischen der äußeren 
Strömung und dem Rezirkulationsgebiet führten (Bourdon und Dutton, 2001, 
Janssen und Dutton, 2005). Im Vergleich zum längsangeströmten Zylinder 
entnahm die Scherschicht einen größeren Anteil des Fluides aus dem Rezir- 
kulationsgebiet (Bourdon und Dutton, 2001, Janssen und Dutton, 2005). Um 
die Massenerhaltung im Rezirkulationsgebiet zu gewährleisten, wurde die dem 
Rezirkulationsgebiet entnommene Masse im Bereich des hinteren Stagnations- 
punktes in das Rezirkulationsgebiet zurückgeführt. Der stromaufgerichtete 
Strömungsimpuls an der Achse wurde folglich durch die stärkere Durchmi- 
schung in der Scherschicht erhöht. 


Insgesamt ergaben sich im Geschwindigkeitsfeld der Nachlaufströmung hinter 
dem Modell mit A = 0° im Vergleich zum längsangeströmten Zylinder Verän- 
derungen in der Größenordnung von bis zu 16 %. Die Struktur der Nachlauf- 
strömung wurde durch die Präsenz der Leitflächen jedoch nicht entscheidend 
verändert. 


Abbildung 5.3b zeigt, dass die Strömung, die an der Hinterkante des Modells 
mit A = 16° ablöste, die zentrale Achse des Strömungsfeldes nicht erreichte. 
Die äußere Strömung wurde ab x/R ~ 15 und r/R ~ 0,6 parallel zur 
Modellachse ausgerichtet. Hinter dem Rekompressionsstoß hatte sich der 
Abstand der Scherschicht zur Rotationsachse gegenüber der Modellhinter- 
kante nahezu halbiert, sodass sich die Azimutalgeschwindigkeit aufgrund der 
Drehimpulserhaltung in etwa verdoppelte. Aufgrund der hohen Azimutal- 
geschwindigkeit wurde die Anzahl der für die PIV-Messungen notwendigen 
Partikel im Rotationszentrum verringert. In der Folge waren ab x/R > 3 
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Abbildung 5.3: Zeitlich gemittelte, relative Axialgeschwindigkeiten (PIV-Messung) im Nachlauf 
der Modelle mit Leitflächen (LF), eingezeichnet sind die zeitlich gemittelten Stromlinien, die 
Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten (——) sowie das Ende der 
Lavaldüse (+++ ), der Expansionsfächer an der Modellhinterkante (---), die Rekompression 
(— - —) und die Isotache U, = 0 (Weidner et al., 2017). 


im Bereich nahe der Achse nicht ausreichend Partikel vorhanden, um eine 
verlässliche Geschwindigkeitsmessung vornehmen zu können, was zu den 
Abbildungen 5.3 und 5.4 sichtbaren Bereichen ohne Messergebnis im fernen 
Nachlauf führte. 


In Abbildung 5.3b ist außerdem zu erkennen, dass ein Teil der Scherschicht- 
strömung durch die Rekompression in Richtung der Basis umgelenkt wurde. 
In der Folge bildete sich für das Modell mit A = 16° die maximale relative 
Rückströmung von —0,15 + 0,02 bei x/R = 0,30 + 0,05 aus. Im Gegensatz 
zum längsangeströmten Zylinder und dem Modell mit A = 0° befand sich 
die maximale Rückströmung nicht mehr an der Achse, sondern bei r/R = 
0,67 + 0,05. Zusammen mit der Scherschicht bildete diese Rückströmung unter 
Vernachlässigung der zusätzlich vorhandenen Azimutalgeschwindigkeiten den 
in Abbildung 5.3b gezeigten äußeren Wirbelring 413. Außerdem ist ein gegen- 
läufig rotierender, innerer Wirbelring {œ direkt an der Achse zu erkennen. Der 
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x/R 


Abbildung 5.4: Zeitlich gemittelte, relative Radialgeschwindigkeiten (PIV-Messung) im Nachlauf 
der Modelle mit Leitflächen (LF), eingezeichnet sind die zeitlich gemittelten Stromlinien, die 
Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten (——) sowie das Ende der 
Lavaldüse (+: ), der Expansionsfächer an der Modellhinterkante (---), die Rekompression 
(— - —) und die Isotache U, = 0 (Weidner et al., 2017). 


innere Wirbelring {Ð resultierte, wie bereits in Abbildung 5.1b gezeigt wurde, 
in einer Umkehr der radialen Strömungsrichtung in der Nähe der Basis. Die 
in den Abbildungen 5.3b und 5.4b gezeigten Axial- und Radialgeschwindig- 
keiten wurden dabei von der azimutalen Stromungsgeschwindigkeit überlagert, 
deren Effekte bereits durch die Ölanstrichmethode gezeigt wurden. Da die 
Strömung der Scherschicht die Strömungsachse nicht erreichte, entstand durch 
den inneren Wirbelring {2 eine bis in den fernen Nachlauf reichende, strom- 
abgerichtete Strömung G), die in dem von Hruschka und Leopold (2015) 
simulierten Wirbelschlauch entlang der Achse resultierte. 


Im Gegensatz zum Modell mit A = 16° zeigt Abbildung 5.3c für A = 32° an der 
Achse eine stromaufgerichtete Strömung vom fernen bis in den nahen Nachlauf. 
Die maximale relative Rückströmung von —0,37 + 0,02 bei x/R = 1,09 + 0,05 
befand sich im Gegensatz zu dem Modell mit A = 16° in der Nähe der 
Strömungsachse. Lediglich in der Nähe der Basis existierte an der Zentralachse 
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ein lokales Strömungsgebiet mit U, > 0, welches durch den für die PIV- 
Methode notwendigen Partikelmassenstrom verursacht wurde (Anhang B), der 
durch das Zentrum der Modellbasis in das Nachlaufgebiet eingebracht worden 
ist. 


Die Scherschicht, die sich stromab der Hinterkante des Modells mit A = 32° 
ausbildete, war kaum in Richtung der Zentralachse geneigt und wurde bei 
x/R x 1 und r/R = 0,8 wieder nahezu parallel zur Modellachse ausge- 
richtet. Aufgrund der geringen Neigung der Scherschicht war die Machzahl 
senkrecht zum Rekompressionsstoß gering, sodass die Strömung nur gering- 
fügig komprimiert wurde. Deshalb konnte fast die gesamte Strömung den 
entstehenden Druckgradienten mühelos überwinden. Das Fluid strömte als 
Wirbelschlauch stromab und umgab ein zentrales Totwassergebiet. Wie im 
Falle des klassischen Nachlaufes ergab sich als Folge der Massenerhaltung 
im Totwassergebiet eine Rückströmung an der Zentralachse. Die durch die 
Scherschicht mitgenommene Fluidmasse wurde an der Achse in Richtung 
der Basis zurückgeführt. Die Strömungsstruktur des Nachlaufes hinter dem 
Modell mit 4 = 32° unterscheidet sich sowohl von dem klassischen Rezir- 
kulationsgebiet als auch von der für A = 16° gezeigten Nachlaufstruktur mit 
stromabgerichtetem Wirbelschlauch entlang der Achse. Das Strömungsfeld 
weist jedoch Ähnlichkeiten mit der Struktur von Wirbelschläuchen auf, wie 
sie in Überschallströmungen durch die Interaktion mit Verdichtungsstößen 
entsteht. Für das Modell mit A = 32° könnte dies durch die an den Windkanal- 
wänden reflektierten Verdichtungsstöße der Leitflächen ausgelöst worden sein, 
die jedoch erst außerhalb des PIV-Messfeldes mit der rotierenden Strömung 
an der Modellachse interagierten. Dies wurde in Kapitel 7.1.2 mittels numeri- 
scher Simulationen genauer untersucht. Auch ohne Interaktion des rotierenden 
Nachlaufes mit Verdichtungsstößen muss die Massenerhaltung im Totwasser- 
gebiet gewährleistet sein, sodass ein Teil der Strömung im fernen Nachlauf in 
Richtung der Basis umgelenkt werden muss. 


Des Weiteren zeigt Abbildung 5.4, dass die von den Leitflächenvorderkanten 
erzeugten Verdichtungsstöße mit ansteigendem Leitflächenanstellwinkel A zu 
einer höheren Radialgeschwindigkeit in der äußeren Strömung führten. Für 
A = 32° istder Verdichtungsstoß selbst nicht erkennbar, da dieser, wie bereits in 
Abbildung 5.1c gezeigt wurde, aufgrund des hohen Leitflächenanstellwinkels 
A vollständig außerhalb des PIV-Messfeldes lag. 
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5.3 Druckfeld 


5.3.1 Statischer Druck vor der Modellhinterkante 


Die Abbildungen 5.5a bis 5.5c zeigen die mittels PSP-Methode ermittelten 
und mit dem Druck pœ vor den Leitflächen normierten Oberflächendruck- 
verteilungen im Bereich der Leitflächen mit unterschiedlichen Anstellwinkeln 
A. Der Druckanstieg vor den Leitflächen wurde durch Leitflächen mit höherem 
Anstellwinkel A verstärkt. Wie bereits in den Abbildungen 5.1 und 5.4 zu 
erkennen war, befand sich der durch die abgelösten Stöße an den Leitflächen- 
vorderkanten hervorgerufene Ablösungsstoß außerdem weiter stromauf. 


Unmittelbar vor den Leitflächen ist in Abbildung 5.5a ein zweiter Druckanstieg 
zu erkennen, welcher aus der in Abbildung 5.2a skizzierten Verzögerung der 
abgelösten Strömung auf Unterschall (M < 1) resultierte. Außerdem ist die 
Interaktion der abgelösten Verdichtungsstöße mit der Grenzschicht an den 
Oberflächen der benachbarten Leitflächen zu erkennen. Der Druckanstieg 
aufgrund des Verdichtungsstoßes resultierte in der lokalen Strömungsablösung, 
die bereits in Abbildung 5.1a gezeigt wurde. Stromab lag die Strömung wieder 
an der Leitflächenoberfläche an, sodass der dort gemessene Druck nahezu dem 
statischen Druck vor den Leitflächen entspricht. Da die Leitflächen in Abbil- 
dung 5.5a keinen Anstellwinkel haben, ist die Umströmung der Leitflächen 
symmetrisch. 


Im Gegensatz zu dem Modell mit 4 = 0° ist für die angestellten Leitflächen in 
den Abbildungen 5.5b und 5.5c ein gegenüber der freien Anströmung erhöhter 
Druck an den der Strömung zugewandten Seiten und ein verminderter Druck 
an den der Strömung abgewandten Seiten der Leitflächen zu erkennen. Die 
Umströmung der Leitflächen ist vergleichbar mit einer schrägangeströmten 
ebenen Platte. Auf der der Strömung zugewandten Seite der ebenen Platte 
bildet sich an der Vorderkante ein Verdichtungsstoß aus, der die Strömungs- 
geschwindigkeit verringert. Auf der der Strömung abgewandten Seite bildet 
sich an der Vorderkante der ebenen Platte eine Expansion aus, sodass die 
Strömung beschleunigt wird. Die unterschiedlichen Drücke auf den beiden 
Seiten einer einzelnen ebenen Platte werden erst stromab der Hinterkante 
ausgeglichen, indem die Strömung auf der Überdruckseite expandiert und auf 
der Unterdruckseite komprimiert wird. 
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Abbildung 5.5: Zeitlich gemittelte, mit Pœ normierte, Druckverteilung an der zylindrischen 
Modelloberfläche sowie den Leitflächen, gemessen mittels PSP-Methode (mit und ohne zylin- 
drische Modellverlängerung). 
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Für das Modell mit A = 16° interagierten die an zwei benachbarten Leitflä- 
chenvorderkanten entstehenden Expansionen und Verdichtungsstöße jedoch 
miteinander. Der Verdichtungsstoß komprimierte dabei die Strömung auf der 
Unterdruckseite und der Expansionsfächer verminderte den Druck auf der 
Überdruckseite. Auf diese Weise wurden die Drücke auf den jeweiligen Seiten 
einer Leitfläche nahezu angeglichen, bevor die Strömung an der Leitflächen- 
hinterkante ablöste. 


Da der Anstellwinkel der Leitflächen von A = 32° für das in Abbildung 5.5c 
gezeigte Modell größer war als der Prandtl-Meyer-Winkel bei M = 2 von 
etwa 26° (Oswatitsch, 1952), löste die Strömung auf der Unterdruckseite der 
Leitflächen ab. Dadurch verringerte sich gegenüber dem Modell mit A = 16° 
der Strömungsquerschnitt zwischen den Leitflächen, sodass der abgelöste Stoß 
an der Leitflächenvorderkante verstärkt wurde. In der Folge befand sich auch 
der Ablösungsstoß vor den Leitflächenvorderkanten weiter stromauf. 


Um den Einfluss der Strömung zwischen den Leitflächen auf den statischen 
Druck in der Ebene der Strömungsablösung genauer untersuchen zu können, 
wurden die Modelle hinter den Leitflächen mit einem 30 mm langen Zylinder 
verlängert. Der Durchmesser der Modellverlängerung entsprach dem Basis- 
durchmesser D der Modelle. Der zusätzlich angebrachte zylindrische Teil 
wirkte sich, wie in den Abbildungen 5.5d bis 5.5f gezeigt, kaum auf die 
Strömung zwischen den Leitflächen aus. 


In Abbildung 5.5d ist das Druckfeld des Modells mit A = 0° zusammen mit 
der Modellverlängerung zu sehen. Hinter den Leitflächen ist auf der Modell- 
verlängerung eine lokale Expansion aufgrund der finiten Dicke der Leitflächen 
zu erkennen. Der gemessene Druck unmittelbar hinter den Leitflächen ist für 
0 < x/R < 0,5 im Mittel etwa 10 % geringer als vor den Leitflächen. 


Abbildung 5.5e zeigt, dass der mittlere Druck an der Oberfläche der Modellver- 
längerung für A = 16° um fast 30 % geringer ist als vor den Leitflächen. Dies 
ist sowohl durch die stärkeren Verdichtungsstöße an den Leitflächenvorder- 
kanten als auch die eingebrachten azimutalen Stromungsgeschwindigkeiten zu 
erklären, die eine radiale Geschwindigkeitskomponente erzeugten und damit in 
einem gegenüber dem Freistrom verringerten Druck an der Modelloberfläche 
resultierten. Die in Abbildung 5.5e erkennbaren Gebieten unterschiedlichen 
Druckes an der Oberfläche der Modellverlängerung wurden durch die an den 
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Oberflächen der Leitflächen reflektierten Expansionen und Verdichtungsstöße 
hervorgerufen. 


Für A = 32° zeigt Abbildung 5.5f eine Verminderung des mittleren statischen 
Druckes hinter den Leitflächen von etwa 50% gegenüber der Anströmung. 
Wie an dem gegenüber dem Modell mit A = 16° stärkeren Druckanstieg vor 
den Leitflächen zu erkennen ist, stellten die Leitflächen mit 32° Anstellwinkel 
auch aufgrund der Strömungsablösungen an deren Vorderkanten ein größeres 
Hindernis für die Strömung dar. Deshalb wurde die Strömung stromab der 
Leitflächen auch stärker expandiert. Außerdem nahmen wegen des höheren 
Anstellwinkels der Leitflächen die Azimutalgeschwindigkeiten zwischen den 
Leitflächen zu. Wie in Abbildung 5.4 zu sehen ist, wurde die Strömung deshalb 
im Vergleich zu dem Modell mit A = 16° radial stärker abgelenkt, sodass der 
statische Druck an der Modelloberfläche zusätzlich vermindert wurde. 


Mit Hilfe der PSP-Methode wurde gezeigt, dass die Leitflächen den statischen 
Druck im Bereich der Strömungsablösung an der Modellhinterkante abhängig 
vom Anstellwinkel A verminderten. Die hierzu verwendete zylindrische 
Modellverlängerung hatte stromauf der Modellhinterkante keinen messbaren 
Einfluss auf die Strömung. In Kapitel 7.1 wurden die PSP-Ergebnisse außerdem 
zur Validierung der numerischen Simulationen verwendet. 


5.3.2 Basisdruck 


Abbildung 5.6 zeigt den gemessenen radialen Druckverlauf am Heck der 
Strömungskörper im Vergleich zu den Messergebnissen des längsangeströmten 
Zylinders aus Kapitel 4.2.3. Die Abweichungen und resultierenden Unsicher- 
heiten der Messergebnisse sind in den Tabellen 5.1 bis 5.3 angegeben und 
weisen ähnliche absolute Messunsicherheiten auf, wie sie auch für den längs- 
angeströmten Zylinder in Tabelle 4.6 bestimmt wurden. 


Dem in Abbildung 5.5d gezeigten Druckabfall von etwa 10% hinter den 
Leitflachen des Modells mit A = 0° wirkte die verminderte Machzahl 
aufgrund der durch die Leitflächen erzeugten Verdichtungsstöße entgegen. 
Die im Vergleich zum längsangeströmten Zylinder verminderte Machzahl 
resultierte in einem erhöhten Verhältnis zwischen Basisdruck und statischem 
Druck an der Modellhinterkante (Lamb und Oberkampf, 1995). Insgesamt 
jedoch verringerten die Leitflächen aufgrund des vor der Ablösung vermin- 
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Abbildung 5.6: Gemessener, zeitlich gemittelter radialer Druckverlauf am Heck der Strömungs- 
körper mit Leitflächen im Vergleich zum längsangeströmten Zylinder (Weidner et al., 2017). 


derten statischen Druckes den Basisdruck gegenüber dem längsangeströmten 
Zylinder um 6 %, wie in Abbildung 5.6 zu erkennen ist. Der gemessene Druck- 
abfall liegt somit in dem von Moore et al. (1992) gezeigten Bereich von 5 % 
bis 10%, der sich für Leitflächen mit einem Verhältnis von Dicke zu Länge 
von 0,05 bei einer Anströmung mit M = 2 ergibt. 


Für A = 16° ist in Abbildung 5.6 für r/R > 0,5 ein Druckabfall von etwa 
50 % gegenüber den Leitflächen ohne Anstellwinkel zu erkennen. Dies ist zum 
einen durch den in Abbildung 5.5e gezeigten geringeren Druck unmittelbar 
hinter den Leitflächen zu erklären. Zum anderen wurde der Massenabfluss 
aus dem Gebiet unmittelbar hinter dem Strömungskörper durch die stukturelle 
Veränderung des nahen Nachlaufgebietes von einem Totwassergebiet hin zu 
einem Strömungsgebiet ohne klassischen Stagnationspunkt begünstigt, was zu 
einer weiteren Verminderung des Basisdruckes führte. Außerdem ergab sich 
für r/R < 0,5 ein radialer Druckgradient, der in einem Druckminimum im 
Zentrum der Basis resultierte. Dies ist auf die in Abbildung 5.1b visualisierten 
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Tabelle 5.1: Messabweichungen und resultierende Messunsicherheit für den Basisdruck des 
Modells mit A = 0°. 


r/R 0,00 0,25 0,40 0,55 0,70 0,825 
Nexp. 26 20 25 4 5 25 
pe 0,58 0,58 0,60 0,57 0,59 0,58 


S— <0,01 <0,01 <0,01 <0,01 <0,01 <0,01 


bps 0,01 0,01 0,01 0,01 0,01 0,01 

Ups 0,01 0,01 0,01 0,01 0,01 0,01 

Up* o5% 002 0,02 0,02 0,02 0,02 0,02 
= 

U DX ‚95% 

I 4% 4% 4% 4% 4% 4% 
Po 


Tabelle 5.2: Messabweichungen und resultierende Messunsicherheit fiir den Basisdruck des 
Modells mit A = 16°. 


r/R 0,00 0,25 0,40 0,55 0,70 0,825 
Nexp. 13 8 12 4 5 12 
pe 0,14 0,19 0,23 0,26 0,27 0,26 
Sox 0,01 <0,01 0,01 <001  <0,01 <001 
b 
bps 0,01 0,01 0,01 0,01 0,01 0,01 
Ups 0,02 0,01 0,01 0,01 0,01 0,01 
Ups 95% 0,03 0,02 0,02 0,02 0,02 0,02 
Pe 9 
U,x 95% 
> 24% 11% 11% 8% 8% 8% 
Po 
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Tabelle 5.3: Messabweichungen und resultierende Messunsicherheit für den Basisdruck des 
Modells mit A = 32°. 


r/R 0,00 0,25 0,40 0,55 0,70 0,825 


Nexp. 32 33 34 7 1 34 
pe 0,22 025 0,25 0,20 0,22 0,20 
Soe <0,01 <0,01 <0,01 <0,01 - <0,01 
bys 001 001 0.01 0,01 0,01 0,01 
Ups 001 001 0,01 0,01 0,01 0,01 
Upos% 002 0,02 0,02 0,02 0,02 0,02 
ane 10% 8% 9% 11% 9% 11% 
i 


Azimutalgeschwindigkeiten U, in der Nähe der Basis zurückzuführen, die 
einen radialen Druckgradienten 


ð U2 
= =p— (5.1) 
r r 


erzeugten, wie er auch in Gaszentrifugen existiert (Kemp, 2009). 


Bei einem Anstellwinkel der Leitflächen von A = 32° zeigt Abbildung 5.6 
ein gegenüber dem Modell mit A = 16° kaum verändertes mittleres Druck- 
niveau an der Basis. Die Verringerung des Basisdruckes gegenüber dem 
längsangeströmten Zylinder entspricht etwa dem Druckunterschied von 50 %, 
der mittels PSP-Methode bereits im Bereich der Strömungsablösung festge- 
stellt wurde. Im Vergleich zu A = 16° ergab sich eine geringere Diskrepanz 
zwischen dem in Abbildung 5.5f gezeigten Druckabfall unmittelbar hinter 
den Leitflächen und dem Druckabfall an der Basis. Dies könnte auf die 
Interaktion des Nachlaufes mit den an den Windkanalwänden reflektierten 
Verdichtungsstößen zurückgeführt werden, die an den Leitflächenvorder- 
kanten erzeugt wurden. Der resultierende Druckanstieg im Nachlauf würde 
einem gegenüber dem Druckabfall an den Leitflächen erhöhten Druckabfall an 
der Basis entgegenwirken. Außerdem könnte der Druckanstieg in Verbindung 
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mit der erneuten Änderung der Strömungsstruktur zur Wiederausbildung 
eines hinteren Stagnationspunktes geführt haben, der in den PIV-Ergebnissen 
(Abbildungen 5.3c und 5.4c) jedoch nicht zu erkennen ist. Im Vergleich zu den 
Druckmessungen für A = 16° ist der für A = 32° gemessene Druck außerdem 
nahezu unabhängig von der radialen Koordinate r. Dies ist eine Folge der 
gegenüber A = 16° vernachlässigbaren Azimutalgeschwindigkeiten in der 
Nähe der Basis, wie sie aus Abbildung 5.1c hervorgehen. 


Die zusätzlich durchgeführten Basisdruckmessungen mit drucksensitiver Farbe 
(PSP) in Kapitel 5.5 bestätigen das mit ansteigendem Anstellwinkel vermin- 
derte Druckniveau sowie den radialen Druckgradienten für A = 16°. 


5.4 Ergebnisse für weitere 
Leitflächenanstellwinkel 


Zusätzlich zu den zuvor gezeigten experimentellen Ergebnissen wurden weitere 
Experimente mit Leitflächenanstellwinkeln A von 8° und 24° durchgeführt, 
um die Transitionspunkte zwischen den unterschiedlichen Nachlaufstrukturen 
genauer bestimmen zu können. Die Ergebnisse für A = 8° zeigen ähnliche 
Strömungsstrukturen wie sie für das Modell mit A = 16° beobachtet wurden 
und bestätigen somit die für A = 16° beschriebenen strukturellen Verände- 
rungen des Strömungsfeldes. Die Transition von der klassischen Nachlauf- 
struktur zu der veränderten Struktur, wie sie für A = 8° und A = 16° beobachtet 
wurde, findet somit für Leitflächenanstellwinkel von weniger als 8° statt. Die 
Ergebnisse für das Modell mit einem Leitflächenanstellwinkel von 24° liegen 
hingegen in der Nähe des Transitionspunktes zwischen den für A = 16° und 
A = 32° beobachteten Nachlaufstrukturen. 


5.4.1 Orientierung des wandnahen Strömungsfeldes 


Abbildung 5.7 zeigt die Orientierung der wandnahen Strömungsfelder für die 
Modelle mit Leitflächenanstellwinkeln von 8° und 24°. In Abbildung 5.7a ist 
an der Basis des Modells mit A = 8° ein spiralförmiges Stromungsmuster zu 
erkennen, wie es zuvor bereits für einen Leitflächenanstellwinkel A von 16° 
(Abbildung 5.1b) beobachtet wurde. Da das Ölgemisch, das auf die Basis des 


80 


5.4 Ergebnisse für weitere Leitflächenanstellwinkel 


(a) A= 8° (b) A = 24° 


Abbildung 5.7: Orientierung der wandnahen Strömungsfelder für Strömungskörper mit Leit- 
flächenanstellwinkeln A von 8° und 24° (Weidner et al., 2016). 


Modells mit A = 8° aufgetragen wurde, eine höhere Viskosität hatte als das für 
das Modell mit A = 16° verwendete Ölgemisch, ist das Muster jedoch weniger 
deutlich zu erkennen. 


Für A = 24° zeigt Abbildung 5.7b einen charakteristischen Ablösungsring 
an der Basis, wie er auch in Abbildung 5.1c für A = 32° zu sehen ist. Der 
Radius des Ablösungsringes ist mit r/R ~ 0,6 geringer als für das Modell mit 
A = 32°. Da der Leitflächenanstellwinkel von 24° geringfügig kleiner als der 
Prandtl-Meyer-Winkel einer Strömung mit M = 2 von etwa 26° (Oswatitsch, 
1952) war, kam es im Vergleich zu dem Modell mit A = 32° nur zu kleineren 
Ablösungen an der Leitflächenvorderkante. Die Leitflächen stellten somit ein 
geringeres Hindernis für die Strömung dar, sodass der Verdichtungsstoß vor 
den Leitflächen schwächer war als bei dem Modell mit A = 32°. Deshalb 
ist auch der in Abbildung 5.7b erkennbare Abstand des Ablösungsstoßes von 
der Leitflächenvorderkante im Vergleich zu dem in Abbildung 5.1c gezeigten 
Abstand für das Modell mit A = 32° deutlich geringer. 
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5.4.2 Geschwindigkeitsfeld im Nachlauf 


Die Abbildungen 5.8 und 5.9 zeigen die PIV-Messergebnisse der axialen und 
radialen Geschwindigkeitsfelder für die Modelle mit Leitflächenanstellwinkeln 
A von 8° und 24°. 


Für A = 8° zeigen die Abbildungen 5.8a und 5.9a ein ähnliches Stromungsbild 
wie die Abbildungen 5.3b und 5.4b für A = 16°. Die geringeren Azimutal- 
geschwindigkeiten, im Vergleich zu dem Modell mit A = 16°, führten in 
einem mit dem Fluidelement dV mitrotierenden Koordinatensystem zu gerin- 
geren Zentrifugalkräften auf das Volumenelement dV. Deshalb wurde die 
Scherschicht im Vergleich zu dem Modell mit A = 16° erst bei geringerem 
Abstand zur Zentralachse und somit weiter stromab bei x/R ~ 2,5 parallel 
zur Modellachse ausgerichtet. Der Bereich an der Achse, in dem aufgrund der 
verminderten Partikelanzahl keine Geschwindigkeitsmessungen durchgeführt 
werden konnten, beginnt aufgrund der geringeren Azimutalgeschwindigkeiten 
erstbeix/R > 4,5. 


Die Abbildungen 5.8b und 5.9b zeigen für A = 24° eine deutliche Ablenkung 
der Scherschicht in Richtung der Modellachse und unterscheiden sich somit von 
den Ergebnissen für A = 32°. Der auf diese Weise entstandene, in Strömungs- 
richtung ansteigende, Druckgradient in der Scherschicht war deshalb in der 
Lage, einen Teil der Strömung in Richtung der Basis umzulenken. Somit 
ist das gemessene Geschwindigkeitsfeld in diesem Bereich den Ergebnissen 
für die Modelle mit Leitflächenanstellwinkeln von 8° und 16° ähnlicher als 
dem in den Abbildungen 5.3c und 5.4c gezeigten Strömungsfeld für das 
Modell mit A = 32°. Andererseits erstreckte sich das stromaufgerichtete 
Strömungsgebiet direkt an der Basis bis in die Nähe der Achse, sodass sich 
der Nachlauf des Modells mit 2 = 24° diesbezüglich von den Ergebnissen 
für geringere Leitflächenanstellwinkel A von 8° und 16° unterscheidet. Die 
gemessenen Geschwindigkeiten in der Nähe der Basis zeigen zudem einen im 
Vergleich zu Abbildung 5.7b größeren Ablösungsring bei r/R ~ 0,75. Mögli- 
cherweise wurde das Messergebnis für das Modell mit A = 24° durch den für 
die PIV-Messungen notwendigen Partikelmassenstrom im Zentrum der Basis 
beeinflusst. Der zentrale Massenstrom könnte, obwohl für die in Anhang B.2 
untersuchten Konfigurationen kein Einfluss des Massenstroms auf die globale 
Strömungsstruktur festgestellt wurde, für das Modell mit A = 24° die zentrale 
stromaufgerichtete Strömung im fernen Nachlauf daran gehindert haben, bis 
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Abbildung 5.8: Zeitlich gemittelte, relative Axialgeschwindigkeiten im Nachlauf der Modelle mit 
A = 8° und A = 24° (PIV-Messung), eingezeichnet sind die zeitlich gemittelten Stromlinien, die 
Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten (——) sowie durch das Ende 
der Lavaldüse (------- ), der Expansionsfächer an der Modellhinterkante (— — —), die Rekompression 
(— - —), die Isotache Ux = 0, das gegenläufig rotierende torusförmige Wirbelpaar Ay Q> und 
der zentrale Wirbelschlauch ©) (Weidner et al., 2017). 
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x/R 


x/R 


Abbildung 5.9: Zeitlich gemittelte, relative Radialgeschwindigkeiten im Nachlauf der Modelle mit 
A = 8° und A = 24° (PIV-Messung), eingezeichnet sind die zeitlich gemittelten Stromlinien, die 
Verdichtungsstöße hervorgerufen durch die Leitflächenvorderkanten (——) sowie durch das Ende 
der Lavaldüse (+--+: ), der Expansionsfächer an der Modellhinterkante (- — —), die Rekompression 
(— - —), die Isotache U, = 0, das gegenläufig rotierende torusförmige Wirbelpaar 11 2 und 
der zentrale Wirbelschlauch ©) (Weidner et al., 2017). 
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zur Basis vorzudringen, und somit die Struktur des Nachlaufs maßgeblich 
verändert haben. 


5.4.3 Druckfeld 


Statischer Druck vor der Modellhinterkante 


Abbildung 5.10 zeigt die mit dem statischen Druck vor den Leitflächen 
normierte Druckverteilung an der Oberfläche der Modelle mit einem Anstell- 


fe 
a 
m 
-3 -2 -1 0 1 
x/R 
(c) A = 8° mit Verlängerung 
x 
Fre, 
Sn 
-3 -2 -1 0 -3 -2 -1 0 1 
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(b) A = 24° (d) A = 24° mit Verlängerung 
P/P% 
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Abbildung 5.10: Zeitlich gemittelte, normierte Druckverteilung vor der Modellhinterkante 
gemessen mittels PSP-Methode für Leitflächenanstellwinkel A von 8° und 24°. 
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winkel von 8° und 24°. Für A = 8° ist in Abbildung 5.10a das durch den 
abgelösten Stoß an den stumpfen Vorderkanten der Leitflächen entstehende 
und in Abbildung 5.2a skizzierte Stoßsystem, bestehend aus Ablösungsstoß 
und nahezu senkrechtem Verdichtungsstoß, zu erkennen. Auf der zylindri- 
schen Oberfläche zwischen den Leitflächen sind auch die an den Vorderkanten 
entstehende Expansion auf der Unterdruckseite und die Kompression auf 
der Uberdruckseite der Leitflächen sichtbar. Bei x/R ~ -0,5 überlagerten 
sich die Expansionen und die Kompressionen, sodass der Druckunterschied 
zwischen der Über- und Unterdruckseite der Leitflächen ausgeglichen wurde. 
Unmittelbar hinter den Leitflächen ist in Abbildung 5.10c eine Verminderung 
des statischen Druckes, im Vergleich zur freien Anströmung, von etwa 20 % 
aufgrund der azimutalen Strömungsablenkung zu erkennen. 


Aufgrund des hohen Anstellwinkels der Leitflächen von A = 24° ist in Abbil- 
dung 5.10b im Vergleich zu geringeren Anstellwinkeln A ein stärkerer Druckan- 
stieg vor den Leitflächen zu erkennen. Unmittelbar hinter den Leitflächen wurde 
der statische Druck hingegen um etwa 40 % gegenüber der freien Anströmung 
vermindert und liegt somit zwischen den gemessenen statischen Drücken für 
die Modelle mit Anstellwinkeln von 16° und 32°. 


Basisdruckmessungen 


In Abbildung 5.11 wird der radiale Verlauf des relativen Basisdruckes für 
den längsangeströmten Zylinder und die Modelle mit den Leitflächenanstell- 
winkeln von 0°, 8°, 16°, 24° und 32° verglichen. Die Messabweichungen und 
resultierenden Unsicherheiten der Druckmessungen für die Modelle mit A = 8° 
und A = 24° sind in den Tabellen 5.4 und 5.5 angegeben. 


Der relative Basisdruck, der für das Modell mit A = 8° gemessen wurde, 
liegt zwischen den zuvor gemessenen Werten für die Modelle mit A = 0° und 
A = 16°. Der radiale Gradient des Basisdruckes ist aufgrund der geringeren 
Azimutalgeschwindigkeiten im Nachlauf weniger ausgeprägt als für 2 = 16°. 


Für das Modell mit A = 24° wurde, abgesehen von der Druckmessung 
im Zentrum der Basis, der geringste Druck aller Modelle gemessen. Der 
gemessene Druckabfall gegenüber dem längsangeströmten Zylinder beträgt 
bis zu 80 %. Der stärkere Verdichtungsstoß vor den Leitflächen aufgrund der 
ausgeprägteren Strömungsablösungen an den Leitflächenvorderkanten führte 
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Abbildung 5.11: Gemessener, zeitlich gemittelter Druckverlauf am Heck der Strömungskörper 
(Weidner et al., 2017). 


für das Modell mit A = 32° im Vergleich zu dem Modell mit A = 24° zu einer 
geringeren Machzahl unmittelbar vor der Modellhinterkante. Auf diese Weise 
erhöhte sich der Basisdruck relativ zum statischen Druck unmittelbar vor 
der Modellhinterkante (Lamb und Oberkampf, 1995). Außerdem führten die 
stärkeren Verdichtungsstöße durch die Reflexion an den Windkanalwänden zu 
einer stärkeren Kompression der Strömung im Nachlauf. Der höhere Druck im 
fernen Nachlauf resultierte durch die zentrale Rückströmung in einer weiteren 
Erhöhung des Basisdruckes. Der Druckanstieg an der Basis für das Modell 
mit A = 32° im Vergleich zu dem Modell mit A = 24° ist somit eine Folge der 
Strömungsablösungen an den Leitflächenvorderkanten und nicht unmittelbar 
auf höhere Drallraten im Nachlauf zurückzuführen. 
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Tabelle 5.4: Messabweichungen und resultierende Messunsicherheit für den Basisdruck des 


Modells mit A = 8°. 


r/R 0,00 0,25 0,40 0,55 0,70 0,825 
Nexp. 21 15 20 3 3 20 
pe 041 0,43 0,45 0,45 0,47 0,45 
SE 001 <001 <O001 <O01 <001 <001 
b 
bps 0,01 0,01 0,01 0,01 0,01 0,01 
Ups 0,01 0,01 0,01 0,01 0,01 0,01 
Up* o5% 002 0,02 0,02 0,02 0,02 0,02 
Py, 9 
U D. x 95% 
= 6% 5% 5% 5% 5% 5% 
Po 


Tabelle 5.5: Messabweichungen und resultierende Messunsicherheit fiir den Basisdruck des 


Modells mit A = 24°. 


r/R 0,00 0,25 0,40 0,55 0,70 0,825 
Nexp. 16 16 16 0 0 16 
pe 0,20 0,12 0,13 - - 0,08 
are 001 <001 <00l - - <0,01 
b 
bpt 0,01 0,01 0,01 - - 0,01 
Ups 0,01 0,01 0,01 - - 0,01 
Ups os% 0,03 0,02 0,02 - - 0,02 
U tet 95% 
> 14% 18 % 16 % - - 25% 
Po 
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5.5 Zusätzliche Basisdruckmessungen mittels 
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Abbildung 5.12: PSP-Messergebnisse für den Basisdruck (Gauthier et al., 2017) im Vergleich zu 


den Messungen mittels Druckgebern. 


Entsprechend den in Kapitel 4.2.3 beschriebenen Druckmessungen im Bereich 
der Leitflächen wurden auch an der Modellbasis Messungen mit druck- 
sensitiver Farbe (PSP) durchgeführt (Gauthier et al., 2017). Abbildung 5.12 
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vergleicht die radiale Basisdruckverteilung der Messungen mittels PSP- 
Methode und Druckgebern. Die Ergebnisse der PSP-Messungen (Gauthier 
et al., 2017) wurden zeitlich und azimutal gemittelt sowie mit dem in 
Kapitel 4.2.3, ebenfalls mittels PSP-Methode gemessenen, statischen Druck 
im Freistrom normiert. Für die PSP-Messungen wurde in Abbildung 5.12 die 
Standardabweichung opsp der azimutalen Druckverteilung vom Mittelwert 
der Messungen angegeben. Die PSP-Messungen ergaben für r/R < 0.9 keine 
azimutalen Druckverteilungen, die die Messergebnisse der zufällig an der 
Basis verteilten Druckbohrungen (Kapitel 5.3.2 und 5.4.3) beinflusst haben 
könnten. 


Die PSP-Ergebnisse (Gauthier et al., 2017) entsprechen qualitativ den 
Messergebnissen aus Kapitel 5.3.2, die mittels in das Modell integrierter 
Druckgeber gemessen wurden. Der gemessene Basisdruck des Modells mit 
geraden Leitflächen (A = 0°) ist geringer als für den längsangeströmten 
Zylinder. Im Allgemeinen verringerte sich der mittlere Basisdruck mit anstei- 
gendem Leitflächenanstellwinkel A. Eine Ausnahme stellt der Anstieg des 
Basisdruckes des Modells mit A = 32° gegenüber dem Modell mit 2 = 24° 
dar, der sich sowohl in den Messungen mit Druckgebern als auch mittels 
PSP-Methode ergab und in Kapitel 5.4 genauer betrachtet wurde. Außerdem 
zeigen beide Messmethoden für das Modell mit A = 16° einen radialen Druck- 
gradienten. Für das Modell mit A = 8° ergab sich für die Messungen mittels 
Druckgebern ein geringfügiger Druckanstieg vom Zentrum zur äußeren Kante 
der Basis, der in den PSP-Ergebnissen jedoch kaum erkennbar ist. 


5.6 Fazit zu den experimentellen Ergebnissen 


Die experimentellen Ergebnisse für die untersuchten Modelle mit unterschied- 
lichen Leitflächenanstellwinkeln A zeigen die drei in Abbildung 5.13 schema- 
tisch dargestellten unterschiedlichen Strukturen der Nachlaufströmung. 


Für 4 = 0° bildete sich im zeitlichen Mittel die in Abbildung 5.13a skizzierte 
Nachlaufströmung mit einem klassischen Rezirkulationsgebiet aus. Im Ver- 
gleich zu den Messergebnissen für den längsangeströmten Zylinder in Kapitel 4 
hatten die Leitflächen einen stabilisierenden und homogenisierenden Ein- 
fluss auf das Rezirkulationsgebiet gezeigt, sodass die, durch die Düse verur- 
sachte, Strömungsasymmetrie im Nachlauf vermindert wurde. Nichtrotierende 


90 


5.6 Fazit zu den experimentellen Ergebnissen 


Rezirkulationsgebiet 


(a) A = 0° 


Umkehr der axialen 
Strömungsrichtung 


zentraler, 
stromabgerichteter 
EU Wirbelschlauch 


(b) A = 8° und A = 16° 


zentrale 


. Rückströmung 
Beschleunigung der 


inneren Strömung 


(c) A = 32° 
Abbildung 5.13: Skizzen der zeitlich gemittelten Nachlaufstrukturen für Modelle mit unterschied- 
lichen Leitflächenanstellwinkeln A. 
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Leitflächen ohne Anstellwinkel resultierten in einem um etwa 6% vermin- 
derten Basisdruck im Vergleich zum längsangeströmten Zylinder, veränderten 
die Struktur des Nachlaufs jedoch nicht. 


In Abbildung 5.13b ist die Struktur des Nachlaufs für die Modelle mit Leitflä- 
chenanstellwinkeln A von 8° und 16° dargestellt. Die rotierende Strömung löst 
an der Modellhinterkante ab und wird aufgrund des gegenüber der Anströmung 
geringeren Basisdruckes in Richtung der Achse beschleunigt. Die auf diese 
Weise konvergierende Strömung erzeugt einen in Strömungsrichtung anstei- 
genden Druckgradienten, der die axiale und radiale Strömungsgeschwindigkeit 
verringert. Außerdem nimmt die notwendige Zentripetalbeschleunigung für 
eine kreisförmige Bewegung mit konstantem Radius mit abnehmendem 
Abstand der Strömung zur Achse zu, sodass die Strömung nicht bis zur 
zentralen Achse vordringen kann. Ein Teil der Strömung wird deshalb bereits 
vor dem Erreichen der Achse in Richtung der Modellhinterkante umgelenkt. 
Auf diese Weise entsteht ein äußerer Wirbelring <1), wie er auch in den 
Messergebnissen für das Geschwindigkeitsfeld (Abbildungen 5.3b und 5.8a) 
unter Vernachlässigung der azimutalen Geschwindigkeitskomponenten zu 
erkennen ist. Der gegenläufig rotierende, innere Wirbelring {2> lässt das 
Fluid aufgrund der verminderten Zentrifugalkräfte in der Nähe der Basis- 
oberfläche von der Modellhinterkante zum Zentrum der Basis strömen und 
verursacht zusammen mit der Rotation der Strömung um die Modellachse die 
spiralförmige Orientierung des wandnahen Strömungsfeldes. Die Azimutalge- 
schwindigkeiten im nahen Nachlaufgebiet erzeugen außerdem einen radialen 
Druckgradienten an der Basis. An der Achse bildet sich ein stromabgerichteter 
Wirbelschlauch mit geringem axialem Strömungsimpuls aus. 


Für A = 32° bildete sich die in Abbildung 5.13c dargestellte Strömungsstruktur 
aus. Die an der Modellhinterkante ablösende, rotierende Strömung wird nur 
geringfügig in Richtung der Rotationsachse abgelenkt. Die äußere Strömung 
erfährt im Gegensatz zu den Nachlaufströmungen mit kleineren Leitflächen- 
anstellwinkeln A nur geringe Druckgradienten in Strömungsrichtung und wird 
deshalb kaum in Richtung der Basis umgelenkt. Durch die sich ausbildende 
Scherschicht wird jedoch ein Teil der inneren Strömung von der äußeren 
Strömung mitgerissen. Die Massenerhaltung für die innere Strömung resul- 
tiert in einer Rückströmung an der Achse. An der Modellbasis wird die 
zentrale Rückströmung in Richtung der Modellhinterkante umgelenkt. Durch 
den stark ansteigenden Abstand zur Rotationsachse reduziert sich aufgrund der 
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Drehimpulserhaltung eine eventuell vorhandene Winkelgeschwindigkeit der 
zentralen Rückströmung, sodass durch die Visualisierung der Orientierung des 
wandnahen Strömungsfeldes keine Azimutalkomponente festgestellt werden 
konnte. Deshalb weist der Basisdruck im Gegensatz zu dem Modell mit A = 16° 
auch keinen radialen Gradienten auf. Die Rückströmung löst von der Basis 
ab, bevor sie die Außenkante der Basis erreicht, sodass in den Visualisie- 
rungen der Orientierung des wandnahen Strömungsfeldes ein charakteristi- 
scher Ablösungsring zu erkennen ist. 


Die zusätzlich durchgeführten Messungen für A = 24° zeigen den Übergang 
zwischen der in Abbildung 5.13b gezeigten Nachlaufstruktur mit stromab- 
gerichtetem, zentralem Wirbelschlauch und der in Abbildung 5.13c skizzierten 
Nachlaufstruktur mit einer Rückströmung an der Achse. 


Insgesamt erbringen die experimentellen Ergebnisse den Nachweis, dass sich 
die Struktur der Nachlaufströmung aufgrund des eingebrachten Dralls verän- 
derte. Die für A = 8° und A = 16° gemessenen Veränderungen konnten bisher 
nur durch numerische Simulationen von Hruschka und Leopold (2015) gezeigt 
werden. Außerdem konnte mit der in Abbildung 5.13c skizzierten Nachlauf- 
struktur, wie sie experimentell für das Modell mit 2 = 32° nachgewiesen 
wurde, eine weitere Veränderung der Nachlaufströmung gezeigt werden. 
Die stabilitätstheoretischen Untersuchungen von Jimenez-Gonzälez et al. 
(2014) sagen für sehr hohe Drallraten strukturelle Veränderungen einer 
laminaren, inkompressiblen Nachlaufströmung voraus, verweisen jedoch 
auf die Ausbildung von zwei getrennten Rezirkulationsblasen entlang der 
Zentralachse. Somit unterscheiden sich die vorliegenden experimentellen 
Ergebnisse von diesen stabilitätstheoretischen Untersuchungen nicht nur durch 
die turbulenten und kompressiblen Strömungsbedingungen, sondern auch 
durch die sich ausbildenden, grundsätzlich unterschiedlichen Strömungsstruk- 
turen. Im Folgenden wurden deshalb numerische Simulationen durchgeführt, 
mittels derer unter anderem die Ursachen für diese in der Literatur bisher nicht 
beschriebene Strukturveränderung der Nachlaufströmung untersucht wurden. 
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Um einen detaillierten Einblick in das Strömungsfeld zu ermöglichen, wur- 
den die experimentellen Untersuchungen mit numerischen Simulationen 
nachgebildet. Zusätzlich zum Nachweis der strukturellen Veränderungen der 
Nachlaufströmung in Kapitel 5 wurden die experimentellen Ergebnisse für 
die Modelle mit Leitflächenanstellwinkeln von 0°, 16° und 32° in Kapitel 7 
als Grundlage zur Validierung der benutzten Simulationsmethoden und ihrer 
Ergebnisse verwendet. Die validierten Simulationsergebnisse ermöglichten 
anschließend eine tiefergehende Analyse der beobachteten Strömungsstruk- 
turen. Strömungsgrößen, die experimentell nicht oder nur mit hohem Auf- 
wand messbar waren, konnten somit den Simulationsergebnissen entnommen 
werden. Außerdem wurde in Anhang B.2 durch numerische Simulationen 
gezeigt, dass der für die Vermessung des Geschwindigkeitsfeldes mittels PIV- 
Methode notwendige Massenstrom in den nahen Nachlauf die Struktur der 
Nachlaufströmung für die Leitflächenanstellwinkel von 0°, 16° und 32° nicht 
grundlegend veränderte. Im Folgenden werden zunächst die numerischen 
Methoden anhand von Strömungssimulationen für den längsangeströmten 
Zylinder beschrieben. 


Die kompressiblen Nachlaufströmungen wurden kontinuumsmechanisch mit- 
tels Erhaltungsgleichungen für Masse, Impuls und Energie beschrieben. Zur 
Beschreibung des Zusammenhangs zwischen den thermischen Zustandsgrößen 
wurde, basierend auf der Annahme eines thermisch idealen Gases, die allge- 
meine Gasgleichung verwendet (Oswatitsch, 1952). Die Abhängigkeit der 
spezifischen Wärmekapazität des Fluids von der Temperatur wurde durch ein 
Polynom achten Grades im Temperaturbereich zwischen 100K und 1000K 
für trockene Luft beschrieben (ANSYS, 2013). Die thermische Leitfähigkeit 
des Fluids wurde mittels kinetischer Gastheorie und die Viskosität mit der 
Drei-Koeffizienten-Methode für trockene Luft nach Sutherland aus der lokalen 
Temperatur bestimmt (Anderson, 2006). 
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Das resultierende System gekoppelter, nichtlinearer, partieller Differential- 
gleichungen zweiter Ordnung, deren Impulserhaltungsgleichungen als Navier- 
Stokes-Gleichungen bezeichnet werden (Oertel et al., 2009), hat keine allge- 
meingültige analytische Lösung. Deshalb wurden die Erhaltungsgleichungen 
und das Rechengebiet diskretisiert, sodass das Gleichungssystem mittels 
numerischer Verfahren gelöst werden konnte (Ferziger und Peri¢, 2008). Für 
die numerischen Simulationen dieser Arbeit wurde eine Finite-Volumen- 
Methode mit einem kommerziellen, druckbasierten, gekoppelten, impliziten 
Gleichungslöser (ANSYS, 2013) verwendet. 


6.1 Turbulenzmodellierung 


Als Direkte Numerische Simulation (DNS) werden Simulationen mit einer 
räumlichen und zeitlichen Diskretisierung beschrieben, die es ermöglicht 
sämtliche Wirbelskalen aufzulösen (Ferziger und Perić, 2008). Um die klein- 
sten Wirbel räumlich aufzulösen, ist eine hohe räumliche Diskretisierungs- 
dichte nötig, sodass ein hoher Rechenaufwand entsteht. Wegen der beschränk- 
ten Rechenkapazitäten sind DNS-Simulationen heutzutage nur für Strömungen 
mit geringer geometrischer Komplexität und niedriger Reynoldszahl möglich 
(Spalart, 2000, Laurien und Oertel, 2011, Spalart und Venkatakrishnan, 2016). 
Die Abbildungen 2.1 und 2.2 zeigen, dass DNS-Simulationen für Überschall- 
nachlaufströmungen mit einer gegenüber den experimentellen Untersuchungen 
reduzierten Reynoldszahl zur Reduktion der Rechenzeiten qualitativ gute 
Ergebnisse erzielen können. Quantitativ ergeben sich dadurch jedoch deutliche 
Unterschiede zwischen den Ergebnissen der Simulationen und der Experi- 
mente. 


Für Strömungen mit hoher geometrischer Komplexität oder hoher Reynolds- 
zahl kann der Rechenaufwand numerischer Simulationen verringert werden, 
indem das Gleichungssystem um eine Turbulenzmodellierung erweitert wird 
(Laurien und Oertel, 2011). Der Einfluss turbulenter Schwankungen wird 
hierzu von einem Turbulenzmodell beschrieben, sodass im Allgemeinen ein 
gegenüber einer DNS-Simulation gröberes Rechengitter verwendet werden 
kann. Abhängig davon welche Teile der Turbulenz modelliert werden, können 
die Simulationsmethoden, wie in Abbildung 2.3 dargestellt, in unterschiedliche 
Gruppen eingeteilt werden. 
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Eine Möglichkeit zur Turbulenzmodellierung ist die Verwendung Reynolds- 
gemittelter Navier-Stokes-Gleichungen (RANS). Die Strömungsgrößen wer- 
den hierzu in einen zeitlich konstanten Mittelwert und einen im zeitlichen 
Mittel verschwindenden Schwankungswert aufgeteilt. Der Schwankungswert 
beinhaltet sämtliche turbulenten Fluktuationen deren Einfluss auf die Strömung 
mit Hilfe eines Turbulenzmodells beschrieben wird. Für die Turbulenzmodel- 
lierung wird dazu oftmals die Annahme isotroper Turbulenz getroffen (Laurien 
und Oertel, 2011). Im Vergleich zu einer DNS-Simulation ist der Rechen- 
aufwand einer RANS-Simulation um mehrere Größenordnungen geringer, 
sodass auch Strömungen mit hoher Reynoldszahl und geometrischer Komple- 
xität berechnet werden können. Insbesondere Strömungen, die an Körpern 
anliegen, können mittels RANS-Methoden mit einem geringen Rechen- 
aufwand gut beschrieben werden. In Ablösungsgebieten ist die Annahme 
isotroper Turbulenz häufig nicht gültig, weshalb RANS-Simulationen abgelöste 
Strömungen oftmals nicht korrekt beschreiben können. RANS-Methoden 
eignen sich deshalb nur in geringem Maße zur Turbulenzmodellierung in 
Überschallnachlaufströmungen, wie auch in den Abbildungen 2.1 und 2.2 zu 
erkennen ist (Laurien und Oertel, 2011). 


Bei Large-Eddy-Simulationsmethoden (LES) werden die räumliche Diskreti- 
sierung des Rechengebietes und die zeitliche Diskretisierung der Erhaltungs- 
gleichungen so gewählt, dass die großskaligen, nichtisotropen Wirbel aufgelöst 
werden können (Sagaut, 2002). Ausschließlich die kleinsten Wirbel werden 
mit einem Subgrid-Scale-Turbulenzmodell (SGS), basierend auf der Annahme 
isotroper Turbulenz, beschrieben (Kolmogorov, 1991). Die Abbildungen 2.1 
und 2.2 zeigen, dass LES-Methoden sehr gut geeignet sind, um Überschall- 
nachlaufströmung zu beschreiben. Die kleinskalige Turbulenz in der Grenz- 
schicht, die für eine gute quantitative Übereinstimmung der LES-Ergebnisse 
mit den Messergebnissen ebenfalls durch das Rechengitter aufgelöst werden 
muss, resultiert im Vergleich zu RANS-Methoden jedoch in höheren Rechen- 
zeiten. Deshalb sind LES-Methoden für parametrische Studien von Überschall- 
nachlaufströmungen nicht geeignet. 


Um Strömungen, die sowohl abgelöste als auch anliegende Strömungsbereiche 
enthalten, effizient simulieren zu können, wurden hybride Methoden ent- 
wickelt, die eine Mischung aus RANS- und LES-Methoden darstellen. Eine 
dieser hybriden Methoden ist die Detached-Eddy-Simulationsmethode (DES), 
die die großskaligen Wirbel in den abgelösten Strömungsbereichen mittels 


97 


6 Numerische Simulationsmethoden 


LES-Methode auflöst und die kleinskaligen Wirbel sowie die Wirbel in der 
Nähe fester Oberflächen mittels Turbulenzmodellierung beschreibt (Spalart 
et al., 1997). Diese Methode wurde ursprünglich für Ablösungen hinter 
Tragflügeln, basierend auf dem Wirbelviskositätsmodell von Spalart und 
Allmaras (1992), entwickelt. Strelets (2001) erweiterte diese Methode um das 
k-w-SST-Wirbelviskositätsmodell von Menter (1994), um damit Strömungs- 
ablösungen in der Grenzschicht besser berechnen zu können. 


Um den Einfluss der großskaligen, nichtisotropen Wirbel im Nachlauf abbilden 
zu können ohne die kleinskalige Turbulenz auflösen zu müssen, wurde für 
die folgenden numerischen Simulationen eine DES-Methode (Spalart, 2001) 
verwendet. Wie in Abbildung 2.1 gezeigt, sind numerische Simulationen mit 
DES-Methoden in der Lage, Überschallnachlaufströmungen qualitativ gut zu 
beschreiben. Zur Modellierung der kleinskaligen Wirbel und der Turbulenz 
in der Grenzschicht wurde das k-w-SST-Wirbelviskositätsmodell nach Menter 
(1994) verwendet. Das k-w-SST-Modell ist in der Lage, Strömungsablösungen 
abzubilden, die im Bereich der Leitflächen und damit im Gültigkeitsbereich 
der RANS-Formulierung auftreten (Menter, 1994). Die SST-Formulierung 
kombiniert hierfür die Anwendbarkeit der k-w-Formulierung im Bereich der 
reinviskosen Unterschicht innerhalb der Grenzschicht mit der geringen Abhän- 
gigkeit der k-e-Formulierung von den Turbulenzeigenschaften im Freistrom 
(Menter, 1994). 


6.2 Diskretisierung des Rechengebietes 


Das Rechengebiet umfasste den gesamten in Abbildung 3.1 gezeigten Ver- 
suchsaufbau und wurde, wie in Abbildung 6.1 gezeigt, durch eine zylindrische 
Gitterschnittstelle in zwei Teilgebiete unterteilt. Das äußere Teilgebiet beinhal- 
tete die Wände des Windkanals sowie die Modellhalterung und umschloss das 
innere Teilgebiet mit den unterschiedlichen Modellkonfigurationen. Die axiale 
Ausdehnung der Ruhekammer wurde durch die Modellhalterung bestimmt, 
während die axiale Ausdehnung der Messkammer so gewählt wurde, dass 
die Strömung das Rechengebiet mit M > 1 verließ. Auf diese Weise wurde 
sichergestellt, dass die Randbedingungen am Auslass des Rechengebietes keine 
Auswirkungen auf das Nachlaufgebiet hatten. Am Einlass wurde entsprechend 
den mittleren experimentellen Bedingungen ein Ruhedruck von 4,8 - 10° Pa 
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adiabate Wand Auslass 


zylindrische Gitterschnittstelle 


Einlass 


Abbildung 6.1: Aufteilung des verwendeten Rechengitters in ein äußeres und ein inneres Rechen- 
gebiet, die über eine zylindrische Gitterschnittstelle miteinander verbunden sind. 


und eine Ruhetemperatur von 295 K vorgegeben. Die Wände des Windkanals 
und die Modelloberflächen wurden als reibungsbehaftete Wände modelliert. 
Da sich bei den Windkanalversuchen die Modelloberfläche innerhalb von weni- 
gen Sekunden asymptotisch auf anschließend nahezu konstante Temperaturen 
abkühlte, sodass zwischen Strömung und Modell kein Wärmeaustausch mehr 
stattfand, wurden die Wände in den Simulationen außerdem als adiabat appro- 
ximiert. 


Für beide Teilgebiete wurden strukturierte Netze mit hexaederförmigen Zellen 
verwendet. Das äußere Teilgebiet umfasste etwa 2 - 10° Gitterzellen. Für das 
innere Teilgebiet wurden die drei in Tabelle 6.1 aufgeführten Rechengitter 
mit unterschiedlichen räumlichen Auflösungen erzeugt. Um die großska- 
ligen Wirbel mittels LES-Methode auflösen zu können, befand sich die 
höchste Gitterzellendichte, wie in Abbildung 6.2 für die mittlere Gitterauf- 
lösung gezeigt, im Nachlauf der unterschiedlichen Strömungskörper. Die 
Zellen wurden in diesem Bereich möglichst kubisch und gleichmäßig gewählt 
(Spalart, 2001). Der in Tabelle 6.1 angegebene Zeitschritt Ar der DES- 
Simulationen wurde für jedes der Rechengitter so gewählt, dass sich die 
Strömungsinformation im Nachlauf innerhalb eines Zeitschrittes konvektiv um 
maximal eine Zelle ausbreitete und sich die CFL-Zahl zu Cor = UAt/Ax < 1 
ergab. 


Wie in Abbildung 3.6 gezeigt, wurde die Grenzschicht an der zylindrischen 
Modelloberfläche mit einem dimensionslosen Wandabstand des wandnächsten 
Zellzentrums von y* ~ 1 bis in die reinviskose Unterschicht (Schlichting und 
Gersten, 2006) aufgelöst. In Abbildung 3.6 wurde außerdem bereits gezeigt, 
dass die simulierte Grenzschichtdicke ô von etwa 5 mm den experimentellen 
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Abbildung 6.2: Inneres Rechengitter mit mittlerer Gitterauflösung (m = 2) für den längsange- 
strömten Zylinder (a) und das Modell mit A = 0° (b). 
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Tabelle 6.1: Übersicht über die räumliche Diskretisierung des inneren Rechengebietes 


Gitterbezeichnung grob mittel fein 
Gitterindex m 3 2 1 
gesamt Ng 3,5 - 10° 7.106 14 : 106 
2 Anströmung 29 36 44 
=  Leitflächen 22 28 34 
N ne 
3 # naher Nachlauf 87 109 134 
3 ferner Nachlauf 18 23 29 
© 
3 äußere Strömung 11 18 23 
a 
= 3 Leitflächen 26 28 34 
< a 
Nachlauf 74 91,5 116,5 
azimutal 324 396 516 
numerischer Zeitschritt At[s] 3,0 - 1077 2,5 -1077 2,0 -1077 


Werten entspricht. Unter der Annahme, dass die Grenzschicht der Leitflächen 
im Vergleich zur Grenzschicht der zylindrischen Oberfläche einen geringen 
Einfluss auf die Strömungsablösung an der Modellhinterkante und somit auf 
das Nachlaufgebiet hat, wurde die wandnahe Strömung an der Leitflächen- 
oberfläche mittels Wandfunktionen (Launder und Spalding, 1974) model- 
liert, sodass die Rechenzeit minimiert werden konnte (ANSYS, 2013). An 
der Leitflächenoberfläche wurde der wandnormale Abstand des wandnächsten 
Zellzentrums deshalb zu y* ~ 30 gewählt. Für die Modellierung wurde eine 
Kombination der Ansätze von White und Christoph (1971) sowie Huang et al. 
(1993) verwendet. Auf diese Weise konnten sowohl die Kompressibilität der 
Strömung als auch der Einfluss von Druckgradienten berücksichtigt werden 
(ANSYS, 2013), welche beispielsweise durch die an benachbarten Leitflä- 
chenvorderkanten verursachten Verdichtungsstöße und Expansionen erzeugt 
wurden. 
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6.3 Diskretisierung der Erhaltungsgleichungen 


Die Flüsse der Erhaltungsgrößen wurden als lokales Riemann-Problem an 
den Zellflächen beschrieben. Der Impuls der Strömung im Zellzentrum wurde 
mit dem Bounded-Central-Differencing-Verfahren nach Leonard (1991) auf 
die Zellflächen projiziert. Dieses Verfahren basiert auf der von van Leer 
(1977, 1979) beschriebenen Grundstruktur für Diskretisierungsverfahren 
höherer Ordnung nach Godunov (1959). Zur Wahrung der Stabilität und zur 
Vermeidung unphysikalischer, numerischer Oszillationen bei größtmöglicher 
Genauigkeit verwendet das Schema für die Projektionen entweder ein reines 
zentrales Differenzenverfahren zweiter Ordnung, eine Mischung aus dem 
zentralen Differenzenverfahren zweiter Ordnung mit einem Upwind-Verfahren 
zweiter Ordnung (Barth und Jespersen, 1989) oder ein Upwind-Verfahren 
erster Ordnung (ANSYS, 2013). Das Upwind-Verfahren erster Ordnung, 
mit seiner im Vergleich zu den Verfahren zweiter Ordnung hohen numeri- 
schen Diffusion, wird ausschließlich zur Stabilisierung des numerischen 
Verfahrens benötigt, wenn die Verwendung der Diskretisierungsverfahren 
zweiter Ordnung zu unphysikalischen Oszillationen im Strömungsfeld führen 
(Leonard, 1988, Gaskell und Lau, 1988). Das Bounded-Central-Differencing- 
Verfahren (ANSYS, 2013) nach Leonard (1991) zeichnet sich aufgrund des 
verwendeten zentralen Differenzenverfahrens zweiter Ordnung gegenüber 
reinen Upwind-Verfahren, wie sie häufig für Überschallströmungen verwendet 
werden, durch eine geringe numerische Diffusion aus und ist somit für die 
Verwendung innerhalb einer DES-Methode geeignet (ANSYS, 2013). Für 
die Projektion der restlichen Erhaltungsgrößen aus dem Zellzentrum auf 
die Zellflächen wurde ein Diskretisierungsverfahren nach van Leer (1979) 
verwendet. Das Verfahren interpoliert ebenfalls zwischen dem numerisch 
stabileren Upwind-Verfahren zweiter Ordnung (Barth und Jespersen, 1989) und 
dem aufgrund seiner geringeren numerischen Diffusion genaueren zentralen 
Differenzenverfahren zweiter Ordnung (ANSYS, 2013). Für die Projektion des 
Druckes an die Zellwände wurde ein zentrales Differenzenverfahren zweiter 
Ordnung verwendet (ANSYS, 2013). Die zeitliche Diskretisierung wurde 
mittels eines impliziten Verfahrens zweiter Ordnung vorgenommen, welches 
zur Vermeidung unphysikalischer, numerischer Oszillationen zwischen einem 
Upwind-Verfahren erster und zweiter Ordnung interpoliert (ANSYS, 2013). 
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Die Gradienten des Strömungsfeldes wurden mittels der Methode der kleinsten 
Fehlerquadrate, basierend auf den zellzentrierten Strömungsgrößen, berechnet. 
Um die in Verbindung mit den Upwind-Diskretisierungsverfahren zweiter 
Ordnung entstehenden Oszillationen der Gradienten zu minimieren, wurden 
die Gradienten mit einem Minimum-Modulus-Verfahren nach Barth und 
Jespersen (1989) beschränkt. 


6.4 Einfluss der räumlichen 
Diskretisierungsdichte auf 
das Simulationsergebnis 


Zur Bestimmung des Einflusses der Diskretisierungsdichte auf das Simulati- 
onsergebnis wurde der über die Zeitdauer Ary und die Fläche der Basis Ap 
gemittelte Basisdruckkoeffizient 


1 Po 
= —dA-1 
sat y N, (6.1) 
Atm 


C ane 
PET Ning yM2, 


fiir die unterschiedlichen Rechengitter m betrachtet. Das zeitliche Mittelungs- 
intervall der transienten Simulationsergebnisse Arm wurde für alle Rechen- 
gitter zur Erreichung statistischer Konvergenz der Schwankungsgrößen zu 
7,5 - 107° s gewählt. Dies entspricht etwa 10 vollständigen Durchströmungen 
des inneren Rechengitters mit einer Strömungsgeschwindigkeit von U«. 
Bezogen auf den Basisdurchmesser ergab sich eine entdimensionalisierte 
Simulationszeit von (U../D)Aty = 97,5. Die Größenordnung der verwendeten 
entdimensionalisierten Simulationszeit entspricht der entdimensionlisierten 
Simulationszeit, die bei den DNS-Ergebnissen von Sandberg (2012) zur Errei- 
chung statistischer Konvergenz der Schwankungsgrößen notwendig gewesen 
ist. 
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Für eine ausreichende relative Gitterauflösung h*, = Axm/Axm=ı und zeitliche 
Mittelung nähert sich das Ergebnis für ein Gitter m dem Wert Cp, m=o für ein 
theoretisches Gitter m = 0 mit unendlich feiner relativer Auflösung 


lim hy, ~ lim ~—— = 0 (6.2) 


an. Der Einfluss der Diskretisierungsdichte auf das Simulationsergebnis wurde 
durch die von der relativen Gitterauflösungen 4% , dem konstanten Koeffizienten 
a und dem ebenfalls konstanten Exponenten g abhängige Abweichung 


a(hy,)4 = Cop.m = Cpy,m=0 (6.3) 


des Basisdruckkoeffizientens Cam von dem theoretischen Wert Cp, m=0 
beschrieben. Zur Bestimmung der Konstanten a und g wurden die Ergeb- 
nisse Cp,,m für die drei Gitterauflösungen m = 1,2,3 verwendet. Da sich 
die Gitterauflösungen jeweils um den Faktor 7 unterschieden, ergibt sich bei 
asymptotischer Konvergenz die beobachtete Ordnung der Konvergenz zu 


Cpy,m=2-Cpy,m=1 


In (7) 


q = (6.4) 


Außerdem gilt 


C m=2 — C m= 
q = rZ Por (6.5) 
ee 
sodass sich der theoretische Grenzwert für m — 0 durch eine nach Richardson 
(1910) durchgeführte Extrapolation zu 


Cp,,m=2 Fa Cp,,m=1 
eae 


Cpy,.m=0 = Cp,,m=1 — (6.6) 


ergibt (Roache, 1994). 


Für die in Tabelle 6.1 angegebenen Rechengitter wurde 7 zu 1,26 bestimmt 
und liegt damit in dem von Roache (1994) empfohlenen Bereich zur Durch- 
führung einer Richardson-Extrapolation numerischer Ergebnisse. Für den 
längsangeströmten Zylinder ergab sich mit den in Tabelle 6.2 angegebenen 
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Tabelle 6.2: Zeitlich und räumlich gemittelter Basisdruckkoeffizient für verschiedene räumliche 
Diskretisierungsdichten und Modelle, berechnet mittels DES-Methode. 


m h* Cova 
m Zylinder =0@ A= 16° A = 32° 
1 1 -0,124 —0,137 —0,256 -0,301 
2 VE -0,127 -0,139 -0,253 -0,303 
3 alee -0,131 -0,140 -0,256 -0,305 
Gitterindex m 
0 1 2 3 
I I T 
0, 118 - 
£ -0,122 + “a 4 
S -0,124- E 4 
3 : 
2 . 
E ", 
3 -0127-) o Cpm a 4 
a Extrapolation nach = 
nn Richardson (1910) 
extrapolierter Wert 
E fürh* 30 
extrapolierter Wert © 
0,131 | DO rn! o 4 
l | | l 
0 g! 1 7 n? 


Relative Gitterauflösung h* 


Abbildung 6.3: Richardson-Extrapolation des mittleren Basisdruckkoeffizienten Cp, für den 


längsangeströmten Zylinder. 
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Werten die beobachtete Konvergenzordnung zu q = 1,88. Die beobachtete 
Konvergenzordnung ist naturgemäß geringfügig kleiner als die formale 
Ordnung qf = 2 des verwendeten numerischen Schemas (Celik et al., 2005). 
Abbildung 6.3 zeigt die Extrapolation der numerischen Simulationsergeb- 
nisse für den längsangeströmten Zylinder auf ein Gitter mit unendlich feiner 
Auflösung h* — 0. Die Abweichung des auf dem feinsten Gitter m = 1 
berechneten Basisdruckkoeffizienten von dem extrapolierten Wert beträgt 5 %. 
Die sich daraus ergebende Unsicherheit der numerischen Ergebnisse wurde als 
ausreichend betrachtet, um die experimentell beobachteten Druckänderungen 
an der Modellbasis von bis zu 50 % wiedergeben zu können. Durch eine weitere 
Verdoppelung der Anzahl der numerischen Gitterzellen gegenüber dem Gitter 
m = 1 und die daraus resultierende relative Gitterauflösung h* = 77! ist eine 
Verringerung der Abweichung zwischen dem simulierten und dem extrapo- 
lierten Wert auf 3 % zu erwarten. Aufgrund der damit einhergehenden höheren 
Rechenzeiten und der beschränkten vorhandenen Rechenkapazitäten wurde 
in der vorliegenden Arbeit jedoch auf Simulationen mit höheren Gitterauflö- 
sungen verzichtet. 


Für die Modelle mit Leitflächen wurde keine asymptotische Konvergenz mit 
einer Konvergenzordnung in der Größenordnung der formalen Ordnung gr 
festgestellt. Vreman et al. (1996) zeigten, dass für LES-Simulationen eine 
höhere Gitterauflösung nicht zwangsläufig zu kleineren numerischen Fehlern 
führt. Eine solche nichtasymptotische Konvergenz (Celik et al., 2005) zeigen 
die Simulationsergebnisse von Kawai und Fujii (2005) in Abbildung 2.1. 
Entgegen der allgemeinen Erwartung, dass eine LES-Simulation nur einen 
Bruchteil der turbulenten kinetischen Energie einer DNS-Lösung auflösen 
kann, kann ein zu grobes Rechengitter dazu führen, dass die aufgelöste turbu- 
lente kinetische Energie gegenüber einer DNS-Lösung überschätzt wird (Celik 
et al., 2005). Dies passiert, wenn die aufgelösten Schubspannungen und somit 
die Dissipationsrate der aufgelösten turbulenten kinetischen Energie von der 
LES-Simulation unterschätzt werden (Celik et al., 2005). 


Aufgrund der nichtasymptotischen Konvergenz der Simulationen wurde die 
Richardson-Extrapolation für die Modelle mit Leitflächen mit der formalen 
Konvergenzordnung des numerischen Schemas qf = 2 sowie ausschließlich 
basierend auf dem feinen (m = 1) und dem mittleren Gitter (m = 2) durchge- 
führt. 
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Für die Simulationsergebnisse konnte anschließend ein Unsicherheitsbereich 
Cp,,m(1 + GCI =) für den numerisch berechneten Basisdruckkoeffizienten 


angegeben werden. Der Unsicherheitsbereich ist nach Roache (1994) unter 
Berücksichtigung eines Sicherheitsfaktors Fs als Gitterkonvergenzindex 


Bere u _ 
B Fs o für m = 1, 
GClp,.m = BESTEN... (6.7) 
|Coy.m—Cpy.m—i]74 fiir We 2 3 
Cpy.m-1(79-1) ir 


definiert. Da die Konvergenzordnung für den längsangeströmten Zylinder 
direkt aus den Simulationsergebnissen bestimmt werden konnte, wurde der 
Sicherheitsfaktor Fs für dieses Modell zu 1,25 gewählt (Roache, 1997). Für 
die Modelle mit Leitflächen wurde aufgrund der angenommenen Konvergenz- 
ordnung ein höherer Sicherheitsfaktor von Fs = 3 gewählt, um den größeren 
Unsicherheiten in der Bestimmung von Cp,,m=0 gegenüber einer berechneten 
Konvergenzordnung Rechnung zu tragen (Roache, 2003). 


Abbildung 6.4 zeigt die berechneten Basisdruckkoeffizienten der unterschied- 
lichen Modellkonfigurationen in Abhängigkeit von der Anzahl der Gitterzellen. 
Die Fehlerbalken geben für jedes Modell den berechneten GCI), ‚m des jewei- 
ligen Rechengitters m an. Da der GCI,,,m für die Modelle mit Leitflächen 
ausschließlich auf den Gittern m = 1 und m = 2 basiert, wurde für das Gitter 
m = 3 kein Fehlerbalken eingezeichnet. In Abbildung 6.4 ist für die Modelle 
mit Leitflächen zu erkennen, dass die Unterschiede im Basisdruckkoeffizienten 
zwischen den einzelnen Modellen größer sind als der Unsicherheitsbereich der 
numerischen Simulationen. Die Simulationsergebnisse sind somit hinreichend 
genau, um die grundlegenden Veränderungen des Basisdruckes für die unter- 
schiedlichen Modellkonfigurationen mit Leitflächen abbilden zu können. 


Abbildung 6.5 zeigt das axiale Geschwindigkeitsfeld im Nachlauf des längsan- 
geströmten Zylinders für die unterschiedlichen Gitterauflösungen m. Die entdi- 
mensionalisierte Axialgeschwindigkeit wurde für die numerischen Simula- 
tionen durch eine dichtegewichtete azimutale Mittelung der zeitlich gemittelten 
Axialgeschwindigkeit U, über Nọ = 361 äquidistante, azimutale Gitterpunkte 
zu 
No = 
T M (6.8) 
Uc% Lia) Pi 
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0,0 


l 
> 
a 
| 
| 


Basisdruckkoeffizient Cp, 


Anzahl der Gitterzellen Ng 


—8— Zylinder --- A=16° --- Cp, (Pv = 0) = -0,36 
A 1=0° -©- Az=32° 


Abbildung 6.4: Abhängigkeit des Basisdruckkoeffizientens vom Rechengitter, eingezeichnet ist 
auch der theoretisch minimale Basisdruckkoeflizient von —0,36 für pp = 0. 


DES-Simulation, m = 3, zeitlich und azimutal gemittelt 
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4 T T T_T T T T =a Par a) T T T T T Er T T T 
> > >| DES-Simulation, m = 2, zeitlich und azimutal gemittelt | 4 
> s aa 4 

‘ a = > > 
> > 1 
> > 
> > 4 
> > 4 
me => ~ 
8 10 
(b)m=2 
4 T T al T T = ae T T T T T BT T T T 
> So >= DES-Simulation, m = 1, zeitlich und azimutal gemittelt | 4 
> > 2 A 
7 = = 
Y oie = R = = 
> a = 4 
A Tig SA = == 2 => = 7 > > == 
= > 


2 nn. é 
| = SS SOS —— on PIV-Messung, zeitlich gemittelt| | 
| 1 1 1 l ` = È Saa Sen 1 1 l 1 N 1 l L L l 1 1 
0 2 4 6 8 10 


(c) m = | und PIV-Messergebnis 


Ux 


U 


es EZ | i 
xsim x,PIV 


-0,50 0,25 0,00 0,25 0,50 0,75 1,00 1,25 


Abbildung 6.5: Vergleich des auf unterschiedlichen Rechengittern m mittels DES-Methode 
simulierten, zeitlich und azimutal gemittelten, relativen Axialgeschwindigkeitsfeldes für einen 
längsangeströmten Zylinder, eingezeichnet sind die zeitlich gemittelten Stromlinien, der Verdich- 
tungsstoß, hervorgerufen durch das Ende der Lavaldüse (+++ ), der Expansionsfächer an der 
Modellhinterkante (---), die Rekompression (— - —) und die Isotachen UX = 0; zusätzlich 
werden die numerischen Ergebnisse für m = 1 mit den PIV-Ergebnissen verglichen (c). 
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Tabelle 6.3: Quantitativer Vergleich der numerischen Ergebnisse für den längsangeströmten 
Zylinder auf den unterschiedlichen Rechengittern m. 


xsp/R min Ur X min Ux /R X min Ux /xsp 
m=3 2,69 +0,05 -0,25 + 0,01 1,51 + 0,05 (56 + 2) % 
m=2 2,84+0,05 -025 + 0,01 1,76 + 0,05 (62 + 2) fo 
m=1 2,95 +0,05 -0,25 + 0,01 1,77 + 0,05 (60 + 2) % 

SCH, GCs. GO a 3 
m=3 5 fo — — — 
m=2 4% 6% 6% — 
m=1 3% 4% 4% _ 


bestimmt. Qualitativ unterscheiden sich die mit unterschiedlichen Gitterauflö- 
sungen erzeugten numerischen Ergebnisse in Abbildung 6.5 nur geringfügig. 
Quantitativ ergaben sich die in Tabelle 6.3 angegebenen Unterschiede in der 
Position des hinteren Stagnationspunktes xsp sowie der maximalen relativen 
Rückströmgeschwindigkeit min UX und deren axialer Position x|min us. Der 
GCI, basierend auf diesen Größen, liegt für das feine Gitter (m = 1) im 
Bereich des GCIp,,m=1, der auf dem Basisdruckkoeffizienten basiert. Die 
Unsicherheiten der numerischen Ergebnisse in der Größenordnung von +5 % 
entsprechen den in Abbildung 2.1 gezeigten Abweichungen ähnlicher Simula- 
tionen von den experimentellen Ergebnissen von Herrin und Dutton (1994). 
Die Simulation auf dem Gitter m = 1 überschätzte die Länge des Rezirku- 
lationsgebietes gegenüber den PIV-Ergebnissen um 4 % und weist somit eine 
geringere relative Abweichung von den experimentellen Ergebnissen auf als 
die in Abbildung 2.2 gezeigten Simulationsergebnisse vergleichbarer hybrider 
RANS/LES-Methoden. 


Die qualitativ und im Vergleich zu anderen DES-Simulationen in Abbil- 
dung 2.2 auch quantitativ gute Ubereinstimmung der Simulationsergebnisse 
mit den Messergebnissen fiir das relative Axialgeschwindigkeitsfeld in Abbil- 
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dung 6.5c sowie für das relative Radialgeschwindigkeitsfeld in Abbildung 6.6 
mit FE 

* _ dint PiUr,i 

Sim. = — N. (6.9) 
Us D 1 Pi 
zeigt die ausreichende Auflösung des feinen Rechengitters (m = 1) zur Be- 
schreibung der grundlegenden strömungsmechanischen Prozesse des Nach- 
laufs. Im Folgenden wurden die Ergebnisse der DES-Simulationen zur Analyse 
der Strömungsfelder verwendet, die mit dem feinsten Gitter (m = 1) erreicht 
wurden, da dies auch für die detaillierte Visualisierung der Ergebnisse von 
Vorteil ist. 


DES-Simulation, m = 1, zeitlich und azimutal gemittelt 


PIV-Messung, zeitlich gemittelt | | 


U U, piv 


r Sim? 


Abbildung 6.6: Vergleich des auf dem Rechengittern m = 1 mittels DES-Methode simulierten, 
zeitlich und azimutal gemittelten, relativen Radialgeschwindigkeitsfeldes für einen längsange- 
strömten Zylinder mit PIV-Messergebnissen, eingezeichnet sind die zeitlich gemittelten Strom- 
linien, der Verdichtungsstoß, hervorgerufen durch das Ende der Lavaldiise (------- ), der Expansions- 
fächer an der Modellhinterkante (---), die Rekompression (— - —) und die Isotachen UX = 0. 
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6.5 Einfluss der Methode zur 
Turbulenzmodellierung auf 
das Simulationsergebnis 


Relativer Basisdruck p* 


| 
0,2 0,4 0,6 0,8 
Relativer Basisradius r/R 


-:-:- RANS-Simulation, m = 1 © Messung (Kapitel 5.3.2) 
—— DES-Simulation, m = 1 x Leopold (1993) 


p 
oO 


+ JOC, jai + eem) 


Pho m= 


(DES-Simulation, m = 1) 


Abbildung 6.7: Vergleich von RANS- und DES-Methode anhand des zeitlich gemittelten, radialen 
Basisdruckverlaufes eines längsangeströmten Zylinders. 


Abbildung 6.7 vergleicht den in Kapitel 5.3.2 gemessenen radialen Basis- 
druckverlauf des längsangeströmten Zylinders mit Simulationsergebnissen, die 
mit unterschiedlichen Turbulenzmodellierungsverfahren erzielt wurden. Die 
gezeigten DES- und RANS-Simulationen (Anhang A) wurden jeweils auf dem 
feinen Rechengitter (m = 1) durchgeführt. Die RANS-Ergebnisse zeigen an der 
äußeren Kante der Basis einen Druckunterschied von >10 % gegenüber dem 
gemessenen Basisdruck. Dies wurde auch bei ähnlichen RANS-Simulationen 
von Sandberg und Fasel (2006b) sowie Fasel und Sandberg (2006) beobachtet. 
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Aufgrund der zeitlichen Mittelung der RANS-Methode ergab sich im Zentrum 
der Basis ein Stagnationspunkt mit einem lokalen Druckmaximum. Von dort 
wurde die Strömung zur äußeren Kante der Basis beschleunigt, was in dem 
gezeigten radialen Druckgradienten resultierte. In der Nähe der äußeren Kante 
wurde die Strömung aufgrund der Strömungsumlenkung abgebremst, sodass 
der Druck wieder anstieg. Insgesamt sind RANS-Methoden aufgrund der 
zeitliche Mittelung der Strömungsprozesse nicht in der Lage, die experi- 
mentellen Ergebnisse für Überschallnachlaufströmungen hinreichend genau 
wiederzugeben (Baurle et al., 2001). 


Die Ergebnisse der DES-Simulation zeigen im Gegensatz zu den RANS- 
Ergebnissen eine im zeitlichen Mittel nahezu homogene Druckverteilung an 
der Basis und eine gute Übereinstimmung mit den experimentellen Ergeb- 
nissen. Die Unsicherheit des mittels DES-Methode simulierten mittleren Basis- 
druckes wurde aus dem quadratischen Mittel des Gitterkonvergenzindexes 
GCI,,,‚m=ı und der doppelten Standardabweichung 20, „, des in Abbildung 6.7 
gezeigten azimutalen Mittelwertes des Basisdruckes bestimmt. Die Standard- 
abweichung des Mittelwertes ist für die DES-Simulationen des längsange- 
strömten Zylinders gegenüber den numerischen Unsicherheiten vernachläs- 
sigbar klein. Der Vergleich von RANS- und DES-Methode für den längsange- 
strömten Zylinder bestätigt die Ergebnisse aus den Abbildungen 2.1 und 2.2, 
dass turbulenzauflösende Simulationsverfahren besser geeignet sind, um die 
strömungsmechanischen Prozesse einer Nachlaufströmung zu beschreiben. 


Um das Auflösungsvermögen der durchgeführten DES-Simulationen beurtei- 
len zu können, wurde das auf den Geschwindigkeitsfluktuationen U’ basierende 
aufgelöste turbulente Energiespektrum 


(U*(Ki)U(Ki)) 


ae (6.10) 


3 
Eyry(x) =) 2n(KiKi) 
i=l 


in Abhängigkeit von den Wellenzahlen x; der drei Raumrichtungen untersucht, 
die sich aus der Fourier-Transformierten des Geschwindigkeitsfeld U(x;) und 
deren komplexer Konjugation U*(x;) ergaben. Die Auswertung wurde in dem 
in Abbildung 6.8 gezeigten kubischen Bereich hinter den Strömungskörpern 
(0 < x/R < 10,-5 < y/R < 5 und -5 < z/R < 5) vorgenommen. 
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kubischer Auswertungsbereich 
des numerisch aufgelösten 
turbulenten Energiespektrums 


Modellhalterung 


Strömungskörper 


EIO 
Ruhekammer Lavaldüse Messkammer 


Abbildung 6.8: Skizze des Auswertungsbereiches des aufgelösten turbulenten Energiespektrums. 


Abbildung 6.9 zeigt das innerhalb dieses Bereiches aufgelöste und mit dem 
Basisradius R skalierte turbulente Energiespektrum El u lle)R = Evu (K) SR 
(Davidson, 2018), bezogen auf die turbulente Längenskala |, = 27/k, 
für die unterschiedlichen Modellkonfigurationen. Wirbel mit einer großen 
dimensionslosen Längenskala (lk/R = 1) zerfallen in kleinere Wirbel und 
transportieren so die turbulente kinetische Energie zu den kleineren turbu- 
lenten Längenskalen. Für isotrope Turbulenz gilt nach Kolmogorov (1991) 
Ply x ake (Davidson, 2018). Abbildung 6.9 zeigt, dass die turbulenten 
Skalen fiir alle Modellkonfigurationen bis in den Bereich nahezu isotroper 
Turbulenz (lk/R < 1) aufgelöst wurden. Für Wirbelgrößen lk/R < 0,6 
weicht das aufgelöste Energiespektrum von der Gesetzmäßigkeit der isotropen 
Turbulenz nach Kolmogorov (1991) ab, da ein Teil der turbulenten Strukturen 
abhängig von der lokalen Gitterauflösung durch das numerische Verfahren 
nicht mehr aufgelöst werden konnte und deren turbulente kinetische Energie 
somit an das SGS-Turbulenzmodell übertragen wurde. Als maximales lokales 
Auflösungsvermögen des numerischen Gitters wurde die doppelte Kanten- 
länge der Gitterzellen 2Ax = 2 (AV)! > basierend auf dem Zellvolumen 
angenommen. Da das Rechengitter im nahen Nachlauf turbulente Strukturen 
mit einer relativen Wirbelgröße /,/R bis 0,1 und in der äußeren Strömung 
lediglich bis 0,6 auflösen konnte, erfolgte ein gradueller Übergang zwischen 
aufgelöster und modellierter turbulenter kinetischer Energie. Das SGS- 
Turbulenzmodell wurde in der abgelösten Strömung ausschließlich für Wirbel- 
größen /./R < 0,6 verwendet, die Teil der isotropen Feinstruktur-Turbulenz 
waren und sich deshalb mittels Wirbelviskositätsmodellen gut beschreiben 
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Abbildung 6.9: Aufgelöstes turbulentes Energiespektrum in Abhängigkeit von der Wirbelgröße 


für die DES-Simulationen des längsangeströmten Zylinders und der Modelle mit Leitflächen, 


eingezeichnet ist die Abhängigkeit El eh oc g" > für isotrope Turbulenz nach Kolmogorov 


(1991) sowie das typische Auflösungsvermögen des Rechengitters für turbulente Strukturen im 
Nachlauf (- — —) und in der äußeren Strömung (— —). 


ließen (Laurien und Oertel, 2011). Außerdem zeigt Abbildung 6.9, dass für 
höhere Leitflächenanstellwinkel insgesamt mehr turbulente kinetische Energie 
durch das numerische Schema aufgelöst wurde. Die aufgelöste turbulente 
kinetische Energie des Nachlaufs mit A = 32° ist mehr als doppelt so hoch, 
wie die für den längsangeströmten Zylinder. 


115 


7 Numerische Ergebnisse und 
Interpretation der Messergebnisse 


Bevor die Messergebnisse zusammen mit den Simulationsergebnissen in 
Kapitel 7.2 genauer analysiert werden konnten, wurden die experimentellen 
und numerischen Ergebnisse in Kapitel 7.1 zunächst auf Übereinstimmung 
bezüglich der grundlegenden Strömungsprozesse untersucht. 


7.1 Vergleich experimenteller und 
numerischer Ergebnisse 


Für den Vergleich mit den numerischen Ergebnissen wurden im Folgenden 
exemplarisch für die drei, in Abbildung 5.13 dargestellten, unterschied- 
lichen Nachlaufstrukturen die experimentellen Ergebnisse für die Modelle 
mit Leitflächenanstellwinkeln A von 0°, 16° und 32° verwendet. Ziel dieses 
Vergleiches war es den Gültigkeitsbereich der numerischen Simulations- 
ergebnisse für die unterschiedlichen Strömungsstrukturen zu bestimmen. 


7.1.1 Statischer Druck an der Modelloberfläche 
im Bereich der Leitflächen 


Zunächst wurde die Qualität der numerischen Simulationen stromauf der 
Strömungsablösung an der Modellhinterkante untersucht. Nur wenn die 
Simulationsergebnisse nicht bereits vor der Strömungsablösung an der Modell- 
hinterkante maßgeblich von den experimentellen Ergebnissen abweichen, ist 
ein anschließender Vergleich der experimentellen und numerischen Ergebnisse 
im Nachlauf sinnvoll. 
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x/R x/R 
(a) A =0° (b) A = 16° 


P/P» 


Abbildung 7.1: Simulierter, zeitlich gemittelter, mit p.. normierter Druck an der zylindrischen 
Modelloberfläche und an den Leitflächen. 


118 


7.1 Vergleich experimenteller und numerischer Ergebnisse 


Abbildung 7.1 zeigt die simulierten, zeitlich gemittelten, mit pœ normierten 
statischen Druckverteilungen an den Modelloberflächen. Die Simulationser- 
gebnisse zeigen insbesondere bezüglich der Positionen der Grenzschichtab- 
lösungen vor den Leitflächen eine gute Übereinstimmung mit den in Abbil- 
dung 5.5 gezeigten PSP-Messergebnissen. Quantitative Vergleiche zwischen 
den PSP-Ergebnissen und dem simulierten Druck an der Modelloberfläche im 
Bereich der Leitflächen waren aufgrund der in Kapitel 4.2.3 beschriebenen 
Temperaturabhängigkeit der PSP-Ergebnisse nur begrenzt möglich, sodass im 
Rahmen dieser Arbeit darauf verzichtet wurde. 


Von dem in Kapitel 3.1.1 mittels Druckgeber bei x/R = -1,925 gemessenen 
Druck weicht der simulierte, gemittelte statische Druck vor den Leitflächen 
mit 0,59 - 10° Pa um etwa 3 % ab. Die Abweichung liegt somit im Bereich der 
Unsicherheiten von Simulation und Messung. 


Da sich vor der Strömungsablösung an der Modellhinterkante keine maßgeb- 
lichen Abweichungen zwischen dem gemessenen und simulierten Strömungs- 
feld ergeben haben, wurde im Folgenden die Gültigkeit der numerischen Ergeb- 
nisse im Nachlauf überprüft. 


7.1.2 Axiales und radiales Geschwindigkeitsfeld 


Die simulierten, zeitlich und azimutal gemittelten Felder der axialen und 
radialen Geschwindigkeiten wurden in den Abbildungen 7.2 und 7.3 mit 
den experimentellen Ergebnissen in der zentralen Ebene hinter den Modellen 
verglichen. 


Die Abbildungen 7.2a und 7.3a zeigen im Nachlauf des Modells mit 2 = 0° 
eine gute Übereinstimmung der simulierten axialen und radialen Geschwindig- 
keiten mit den PIV-Messergebnissen. Die Abweichung zwischen den experi- 
mentellen und numerischen Ergebnissen von 3 % in der Länge des Rezirkula- 
tionsgebietes liegt innerhalb der experimentellen und numerischen Unsicher- 
heiten. 


Für das Modell mit A = 16° bestätigt das in Abbildung 7.2b gezeigte 
numerische Ergebnis die Ausbildung des experimentell beobachteten, gegen- 
läufig rotierenden, torusförmigen Wirbelpaares sowie der zentralen, stromab- 
gerichteten Strömung in Form eines Wirbelschlauches. Abbildung 7.3b zeigt, 
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DES-Simulation, m = 1, zeitlich und azimutal gemittelt 


PIV-Messung, zeitlich gemittelt, |] 
Agir—opiv = (7 + 3)° 


0 2 4 6 8 10 


U* 


hi SS lll a lCOT—a 
x ,Sim.’?~ x ,PIV 


U 
—0,50 —0,25 0,00 0,25 0,50 0,75 1,00 1,25 


Abbildung 7.2: Vergleich der gemessenen (PIV) und simulierten (DES), relativen Axialgeschwin- 
digkeiten hinter Modellen mit Leitflächen, eingezeichnet sind die zeitlich gemittelten Stromlinien, 
die Verdichtungsstöße, hervorgerufen durch die Leitflachenvorderkanten (——) sowie deren Refle- 
xionen von den Windkanalwänden (— —), der Verdichtungsstoß am Ende der Lavaldiise (------: ), 
der Expansionsfächer an der Modellhinterkante (-—--), die Rekompression (— - —) und die 
Isotachen UX = 0. 


dass die PIV-Messung die radiale Geschwindigkeit unmittelbar hinter dem 
Modell (0 < r/R < 1) gegenüber dem Simulationsergebnis um bis zu 
10% der Freistromgeschwindigkeit U~ überschätzte. Die Partikel der PIV- 
Messungen wurden aus dem Wirbelzentrum an der Modellachse verdrängt, da 
sie aufgrund ihrer im Vergleich zur Strömung höheren Dichte und der damit 
für eine kreisförmige Bewegung notwendigen höheren Zentripetalkräfte eine 
positive radiale Relativgeschwindigkeit gegenüber der Strömung aufwiesen. 


Das in Abbildung 7.2c gezeigte numerische Ergebnis für das Modell mit A = 
32° bestätigt die in den Experimenten gezeigte Umkehr der axialen Strömungs- 
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DES-Simulation, m = T, zeitlich und azimutal gemittelt 


PIV-Messung, zeitlich gemittelt, |] 
Agir—opiv = (7 + 3)° 

l fi fi fi l fi 
0 2 4 6 8 10 


UF sim? UF prv 


Abbildung 7.3: Vergleich der gemessenen (PIV) und simulierten (DES), relativen Radial- 


E O a A ; ; we ; ; ene : 
geschwindigkeiten U* pīy beziehungsweise U, hinter Modellen mit Leitflächen, einge- 


zeichnet sind die zeitlich gemittelten Stromlinien, die Verdichtungsstöße, hervorgerufen durch 
die Leitflächenvorderkanten (——) sowie deren Reflexionen von den Windkanalwänden (— —), 
der Verdichtungsstoß am Ende der Lavaldüse (+ ), der Expansionsfächer an der Modellhinter- 
kante (- ——), die Rekompression (— - —) und die Isotache Už = 0 


richtung an der Achse im Vergleich zum Nachlauf des Modells mit A = 16°. 
In Abbildung 7.3c ist zu erkennen, dass die PIV-Messung auch für das Modell 
mit A = 32° die radialen Geschwindigkeiten gegenüber den numerischen 
Ergebnissen in der Größenordnung von 5 % der Freistromgeschwindigkeit Ux 
überschätzt. Allerdings liegt der Bereich in dem die radialen Geschwindig- 
keiten überschätzt wurden im Gegensatz zu dem Modell mit A = 16° zwischen 
r/R = 1 und r/R = 1,6. Die geringen Abweichungen zwischen den numeri- 
schen und experimentellen Ergebnissen im Bereich 0 < r/R < 1 hinter dem 
Modell mit A = 32° deuten darauf hin, dass die Azimutalgeschwindigkeiten 
dort gegenüber dem Modell mit A = 16° vernachlässigbar sind. 
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Übergang zwischen Düse und Messkammer 
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Abbildung 7.4: Skizze der an den Leitflachenvorderkanten entstehenden und an den Windkanal- 
wänden reflektierten Verdichtungsstöße sowie des VerdichtungsstoBes am Ende der Lavaldüse. 


Abbildung 7.4 skizziert die Verdichtungsstöße, die an den Übergängen 
zwischen Düse und Messkammer entstanden. Diese Verdichtungsstöße sind 
auch in den in Abbildung 7.3 gezeigten experimentellen Ergebnissen zu 
erkennen. Weniger ausgeprägt sind diese Verdichtungsstöße in den numeri- 
schen Ergebnissen, da der Übergang zwischen Düse und Messkammer bei 
dem verwendeten Rechengitter nahezu glatt war. Zusätzlich zu dem äußerst 
schwachen Düsenstoß zeigt Abbildung 7.4 die an den Leitflächenvorderkanten 
entstehenden Verdichtungsstöße. Diese im Vergleich zum Düsenstoß stärkeren 
Verdichtungsstöße wurden an den Wänden des Windkanals reflektiert. Die 
reflektierten Verdichtungsstöße lagen außerhalb des PIV-Messbereiches und 
sind deshalb in den experimentellen Ergebnissen nicht zu erkennen. Die 
numerischen Ergebnisse in Abbildung 7.3 zeigen, dass sich die durch 
die Leitflächen erzeugten Verdichtungsstöße sowie deren Reflexionen von 
den Windkanalwänden mit zunehmendem Leitflächenanstellwinkel weiter 
stromauf befanden. In Abbildung 7.3c ist des Weiteren zu erkennen, dass die 
Nachlaufströmung bei x/R ~ 9 mit den reflektierten Verdichtungsstößen inter- 
agierte. Die Strömung wird durch die reflektierten Verdichtungsstöße zunächst 
radial in Richtung der Achse umgelenkt. Ein Teil des Fluids kann den daraus 
resultierenden ansteigenden Druckgradienten nicht überwinden und strömt 
deshalb anschließend stromauf in Richtung des Modells. Ohne die reflek- 
tierten Verdichtungsstöße würde es im fernen Nachlauf aufgrund der viskosen 


124 


7.1 Vergleich experimenteller und numerischer Ergebnisse 


Reibung zu einer Abschwächung des Wirbels kommen. Somit würden die 
zentrifugalen Volumenkräfte abnehmen und das Fluid in Richtung der Achse 
vordringen können, sodass es vermutlich auch ohne reflektierte Verdichtungs- 
stöße im Nachlauf zu einer Umkehr der axialen Strömungsrichtung kommen 
würde. 


7.1.3 Turbulente kinetische Energie 


Ein kritischer Punkt turbulenzmodellierender numerischer Verfahren ist die 
korrekte Beschreibung der turbulenten Prozesse. Um die Gültigkeit der 
gewählten Modellierungen zu überprüfen, wurden in Abbildung 7.5 die 
gemessene (PIV) und die mittels DES-Methode simulierte, entdimensionali- 
sierte aufgelöste turbulente kinetische Energie 


Kan = e (7.1) 


basierend auf den axialen und radialen Geschwindigkeitsfluktuationen UL 
und U/ im Nachlauf verglichen. Sowohl die PIV- als auch die DES-Methode 
waren in ihren räumlichen Auflösungen durch die verwendeten diskreten Gitter 
begrenzt, sodass weder die numerischen noch die experimentellen Ergebnisse 
die kleinsten turbulenten Schwankungen enthalten. Beide Methoden waren 
jedoch in der Lage, die turbulenten Wirbelstrukturen in der Größenordnung 
l«/R =, 0,1 im Nachlauf aufzulösen, die den größten Anteil an der gesamten 
turbulenten kinetischen Energie beinhalteten. Deshalb sind die experimentellen 
und numerischen Ergebnisse bezüglich der turbulenten kinetischen Energie 
trotz lokal unterschiedlicher Auflösungen näherungsweise vergleichbar. 


Die Abbildungen 7.5a und 7.5b zeigen, dass die turbulenten Wirbelstrukturen 
unmittelbar nach der Ablösung der Grenzschicht an der Modellhinterkante 
von der DES-Simulation nicht aufgelöst wurden. In den numerischen Simula- 
tionen fand an dieser Stelle ein unphysikalischer Übergang zwischen der 
modellierten Turbulenz der Grenzschicht und der aufgelösten Turbulenz im 
Nachlauf statt (Forsythe et al., 2002a). Im Bereich der Rekompression hingegen 
überschätzten die numerischen Simulationen die turbulente kinetische Energie 
im Vergleich zu den experimentellen Ergebnissen. Dies ist zum einen auf 
das verwendete k-w-SST-Wirbelviskositätsmodell zurückzuführen, welches in 
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Abbildung 7.5: Vergleich des gemessenen (PIV) und simulierten (DES), zeitlich gemittelten 
Feldes der relativen turbulenten kinetischen Energie k*, u, eingezeichnet sind die zeitlich 
gemittelten Stromlinien, die Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten 
(——) sowie deren Reflexionen von den Windkanalwänden (— —), der Verdichtungsstoß am Ende 


der Lavaldüse (+++: ), der Expansionsfächer an der Modellhinterkante (— -= -), die Rekompression 
(— - —) und die Isotachen UX = 0 
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Bereichen starker Druckgradienten die turbulente kinetische Energie tenden- 
ziell überschätzt (Oliver et al., 2007, Tan und Jin, 2011). Zum anderen können 
LES-Simulationsverfahren im Vergleich zu DNS-Lösungen bei zu geringer 
Gitterauflösung die Schubspannungen und damit die Dissipationsrate der 
turbulenten kinetischen Energie unterschätzen (Celik et al., 2005). Somit kann 
die turbulente kinetische Energie gegenüber DNS-Lösungen und experimen- 
tellen Ergebnissen überschätzt werden (Celik et al., 2005). 


Der Vergleich von Abbildung 7.5a und Abbildung 7.5b zeigt, dass die 
Leitflachen des Modells mit 4 = 0° die turbulente kinetische Energie in 
der Scherschicht gegenüber dem längsangeströmten Zylinder geringfügig 
erhöhten. Dies führte zu einer veränderten Durchmischung zwischen der 
äußeren Strömung und dem Rezirkulationsgebiet und beeinflusste somit den 
Basisdruck (Sandberg und Fasel, 2006b). Bourdon und Dutton (2002) hatten 
gezeigt, dass azimutale Wirbelstrukturen den Basisdruck erhöhen, während 
in Strömungsrichtung ausgerichtete turbulente Strukturen den Basisdruck 
vermindern (Bourdon und Dutton, 2001, Janssen und Dutton, 2005). 


In Abbildung 7.5c ist auch für das numerische Simulationsergebnis des Modells 
mit A = 16° die im Vergleich zum Experiment geringere turbulente kinetische 
Energie unmittelbar nach der Strömungsablösung zu erkennen. Ansonsten 
stimmen die experimentellen und numerischen Ergebnisse für die Modelle 
mit A = 0° und A = 16° trotz der begrenzten räumlichen Auflösungen von 
PIV- und DES-Methode gut überein. 


Für das Modell mit A = 32° zeigt sich in Abbildung 7.5d sowohl in den experi- 
mentellen als auch den numerischen Ergebnissen ein gegenüber den anderen 
Modellen erhöhtes Niveau der turbulenten kinetischen Energie. Obwohl die 
Simulation für das Modell mit A = 32° die turbulente kinetische Energie 
unmittelbar nach der Ablösung der Strömung überschätzte, zeigen die numeri- 
schen und experimentellen Ergebnisse ansonsten eine gute Übereinstimmung. 


7.1.4 Wandschubspannungen 


Die in Abbildung 7.6 gezeigten Richtungen der simulierten, zeitlich gemit- 
telten Wandschubspannungen rw an der Basis stimmen qualitativ mit den 
Strömungssichtbarmachungen mittels Ölanstrichverfahrens aus Abbildung 5.1 
überein. Quantitativ zeigt Abbildung 7.6b für das Modell mit A = 16°, dass 
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Abbildung 7.6: Zeitlich gemittelte, simulierte Wandschubspannungen an der Basis der Modelle 
mit Leitflächen. 
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die Wandschubspannungen im Vergleich zu dem Modell mit A = 0° in Ab- 
bildung 7.6a größer sind. Dies ist die Folge des höheren Gesamtimpulses 
der Strömung im Nachlaufgebiet für das Modell mit A = 16° aufgrund der 
zusätzlich vorhandenen azimutalen Strömungsgeschwindigkeiten. 


Abbildung 7.6c zeigt die hauptsächlich in positiver radialer Richtung orien- 
tierten Wandschubspannungen für das Modell mit A = 32°. Eine Rückströmung 
in Richtung des Basiszentrums ist nur lokal zu erkennen. Erst bei einem 
zehnmal höheren zeitlichen Mittelungsintervall wäre ein gemitteltes Ergebnis 
zu erwarten gewesen, bei dem die azimutalen Strömungsasymmetrien vernach- 
lässigbar geworden wären (Forsythe et al., 2002b). Um den Rechenaufwand 
gering zu halten, wurde stattdessen eine zusätzliche azimutale Mittelung durch- 
geführt, wie es von Forsythe et al. (2002b) für einen ählichen Fall vorge- 
schlagen wurde. Daraus resultierte das bei r/R=0,87 eingezeichnete lokale 
Minimum der zeitlich und azimutal gemittelten Wandschubspannungen, dessen 
radiale Position um nur 8,8 % von der in Abbildung 5.1c erkennbaren kreis- 
förmigen Strömungsablösung abweicht. Abbildung 7.6c zeigt außerdem in 
Übereinstimmung mit den experimentellen Ergebnissen, dass die Azimutal- 
komponente der wandnahen Strömung im Nachlauf des Modells mit A = 32° 
gegenüber dem Modell mit A = 16° vernachlässigbar klein ist. Die Wandschub- 
spannungen für A = 32° haben somit eine mit dem Modell mit A = 0° 
vergleichbare Magnitude. Das Strömungsgebiet unmittelbar hinter dem Modell 
mit A = 32° lässt sich deshalb analog zu dem Rezirkulationsgebiet des längs- 
angeströmten Zylinders und des Modells mit A = 0° als Totwassergebiet mit 
geringem Strömungsimpuls charakterisieren. 


7.1.5 Basisdruck 


Abbildung 7.7 vergleicht den gemessenen Basisdruckverlauf mit den Ergeb- 
nissen der numerischen Simulation. Die Simulation für das Modell mit A = 0° 
war in der Lage, den niedrigeren Basisdruck im Vergleich zum längsange- 
strömten Zylinder in der gemessenen Größenordnung von 6 % wiederzugeben. 
Für die Modelle mit Leitflächenanstellwinkeln von 16° und 32° reduzierte sich 
der Basisdruck im Mittel auf weniger als die Hälfte des Druckes des Modells 
mit A = 0° und entspricht damit qualitativ den Messergebnissen. 
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Abbildung 7.7: Vergleich der gemessenen Drücke an der Basis mit den mittels DES-Methode 
simulierten, zeitlich und azimutal gemittelten, radialen Basisdruckverläufen des längsangeströmten 
Zylinders und der Modelle mit Leitflächen. 
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Für das Modell mit 4 = 16° bildete sich ein radialer Druckgradient aus, der 
in den numerischen Ergebnissen ausgeprägter ist als in den Messergebnissen. 
Der Druck am äußeren Rand der Basis ist in den numerischen Ergebnissen 
somit höher als in den Messergebnissen, während der Druck im Zentrum der 
Basis geringer ist als in den Experimenten. Der Anstieg der Unsicherheiten des 
Simulationsergebnisses im Bereich des äußeren Randes (r/R > 0,85) resul- 
tiert aus der höheren Standardabweichung des Mittelwertes o,,,, der azimu- 
talen Druckverteilung für dieses Modell, die auch in Abbildung 5.5e auf der 
Modellverlängerung erkennbar ist. Im Umkehrschluss zeigen die vergleichs- 
weise kleinen Unsicherheitsbereiche für r/R < 0,85, deren Magnitude in 
diesem Bereich außerdem hauptsächlich auf die Unsicherheiten der numeri- 
schen Berechnungen (GCI,,) zurückzuführen ist, dass die zufällige azimutale 
Verteilung der Druckbohrungen an der Basis selbst für das Modell mit A = 16° 
keinen entscheidenden Einfluss auf die experimentellen Ergebnisse hatte. 


Für 4 = 32° stimmen Simulation und Messungen bezüglich des gegenüber dem 
Modell mit A = 16° vernachlässigbaren radialen Gradienten gut überein. Beide 
Ergebnisse zeigen ähnlich den Ergebnissen für den längsangeströmten Zylinder 
und das Modell mit A = 0° nur eine geringe Abhängigkeit des Basisdruckes von 
der radialen Position. Die Simulation überschätzte jedoch im Vergleich zu den 
Messungen die Veränderung des Basisdruckes relativ zum längsangeströmten 
Zylinder, sodass der für A = 32° simulierte Basisdruck geringer ist als das 
Messergebnis. 


7.1.6 Fazit zum Vergleich zwischen den Simulationen 
und Messungen 


Die durchgeführten DES-Simulationen zeigen eine gute qualitative und 
teilweise auch quantitative Übereinstimmung mit den experimentellen Ergeb- 
nissen aus Kapitel 5. Die DES-Simulationen waren in der Lage, die in 
Kapitel 5 gemessenen strukturellen Veränderungen der Nachlaufströmung 
abzubilden. Die quantitativen Abweichungen zwischen den experimentellen 
und numerischen Ergebnissen sind im Allgemeinen geringer als die Unter- 
schiede, die sich aufgrund der unterschiedlichen Anstellwinkel der Leitflächen 
ergaben. Die DES-Ergebnisse eignen sich aufgrund ihrer guten qualitativen 
Übereinstimmung mit den Messergebnissen für eine tiefergehende Analyse 
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der strömungsmechanischen Zusammenhänge zur Erklärung der in Kapitel 5 
experimentell nachgewiesenen Veränderungen des Stromungsfeldes. 


Außerdem konnten die Ergebnisse der numerischen Simulationen bereits 
Aufschluss über die Ursache der axialen Strömungsumlenkung im fernen 
Nachlauf des Modells mit A = 32° geben, die mit experimentellen Methoden 
nicht untersucht werden konnte. 


7.2 Interpretation der Mess- und 
Simulationsergebnisse 


Im Folgenden wurden die numerischen Ergebnisse dazu verwendet, die 
Azimutalgeschwindigkeiten und das Druckfeld im Nachlauf zu untersuchen, 
die experimentell nicht vermessen werden konnten. 


7.2.1 Azimutalgeschwindigkeit 


Abbildung 7.8 zeigt die numerischen Ergebnisse der azimutal und zeitlich 
gemittelten relativen Azimutalgeschwindigkeiten 


N, m 
Api pil, yl 

U% = Fr (7.2) 
Us DE Pi 


für die Modelle mit den Leitflächenanstellwinkeln A von 16° und 32°. 


Für das Modell mit A = 16° istin Abbildung 7.8a zu erkennen, dass die relative 
Azimutalgeschwindigkeit unmittelbar nach der Ablösung der Strömung an der 
Modellhinterkante etwa 0,5 beträgt. Anschließend konvergierte die Strömung 
in Richtung der Modellachse, sodass die relative Azimutalgeschwindigkeit 
aufgrund der Drehimpulserhaltung aufbis zu 0,70 anstieg. Ein Teil des Drehim- 
pulses wurde in dem Gebiet mit stromaufgerichteter Strömung in Richtung der 
Basis transportiert und versetzte somit die Strömung unmittelbar hinter dem 
Modell in eine rotierende Bewegung um die Zentralachse. 
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(a) A = 16° 


(b) A = 32° 
U; ME = | TE 
<0,0 0,2 0,4 0,6 0,8 1,0 


Abbildung 7.8: Zeitlich und azimutal gemittelte, relative Azimutalgeschwindigkeit im Nachlauf 
der Modelle mit angestellten Leitflächen, eingezeichnet sind die zeitlich gemittelten Stromlinien, 
die Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten (——) sowie deren Refle- 
xionen von der Windkanalwand (— —), der Verdichtungsstoß am Ende der Lavaldiise (------: ), 
der Expansionsfächer an der Modellhinterkante (- --), die Rekompression (— - —), die Isotache 
U% = 0 und die maximale relative Azimutalgeschwindigkeit U, Fas a 
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(a) A = 16° 


(b) A = 32° 
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Abbildung 7.9: Zeitlich und azimutal gemittelte Drallrate im Nachlauf der Modelle mit an- 
gestellten Leitflächen, eingezeichnet sind die zeitlich gemittelten Stromlinien, die Verdichtungs- 
stöße, hervorgerufen durch die Leitflächenvorderkanten (——) sowie deren Reflexionen von der 
Windkanalwand (— —), der Verdichtungsstoß am Ende der Lavaldüse (------- ), der Expansions- 
fächer an der Modellhinterkante (- - -), die Rekompression (— - —) und die Isotache UX = 0 
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Abbildung 7.8b zeigt die relative Azimutalgeschwindigkeit im Nachlauf des 
Modells mit A = 32°. Bereits unmittelbar nach der Ablösung der Strömung von 
der Modellhinterkante erreichte die Strömung die maximale relative Azimutal- 
geschwindigkeit von 0,76. Da aufgrund der in Strömungsrichtung nur gering- 
fügig ansteigenden Druckgradienten kaum konvektiver Transport des Dreh- 
impulses in das Strömungsgebiet unmittelbar hinter dem Modell mit 2 = 
32° stattfand, war die relative Azimutalgeschwindigkeit von <0,3 in diesem 
Bereich, im Vergleich zu dem Modell mit A = 16°, vernachlässigbar gering. 
Die geringen Azimutalgeschwindigkeiten in der Nähe der Basis des Modells 
mit A = 32° wurden bereits aus den Ölvisualisierungen in Abbildung 5.1 
erkenntlich. 


Der Vergleich der Drallraten der Strömung fro. = UZUoo/(2ar) in Abbil- 
dung 7.9 zeigt, dass die maximal erreichte Rotationsrate von etwa 8000 Hz 
für das Modell mit A = 16° trotz des geringeren Anstellwinkels 2 nahezu 
doppelt so hoch ist wie für das Modell mit A = 32°. Der Drehimpuls konnte 
für A = 16° unmittelbar nach der Modellhinterkante bis in die Nähe der 
Modellachse vordringen, sodass die höchsten Drallraten für dieses Modell im 
Bereich x/R < 2 entstanden. Für das Modell mit A = 32° wurde die maximale 
Drallrate erst bei x/R ~ 10 erreicht, als die Strömung durch die reflektierten 
Verdichtungsstöße der Leitflächenvorderkanten in Richtung der Modellachse 
umgelenkt wurde. 


7.2.2 Radiale Druckgradienten und Zentrifugalkrafte 


Der Vergleich von Abbildung 7.3 mit den Abbildungen 7.2 und 7.8 zeigt, 
dass die Radialgeschwindigkeiten fiir die Modelle mit angestellten Leitflachen 
um nahezu eine Größenordnung geringer waren als die Axial- und Azimutal- 
geschwindigkeiten. Diese kreisförmige Bewegung der Strömung mit einem 
nahezu konstanten Radius um die Zentralachse ist nur möglich, wenn 
innerhalb des Strömungsfeldes eine in Richtung der zentralen Achse wirkende 
Zentripetalkraft pus /rdV auf jedes Volumenelement dV wirkte. In einem 
mit dem Fluidelement dV mitrotierenden Koordinatensystem muss aufgrund 
der dort verschwindenden Azimutalgeschwindigkeit eine in positiver radialer 
Richtung wirkende zentrifugale Volumenkraftdichte pU, a /r eingeführt werden, 
um die Bewegung des Fluidelements dV korrekt beschreiben zu können. 
In Gaszentrifugen ergibt sich in diesem mitrotierenden Koordinatensystem 
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beispielsweise ein Gleichgewicht zwischen der zentrifugalen Volumenkraft- 
dichte und dem radialen Druckgradienten (Kemp, 2009). Der Zusammenhang 
zwischen Azimutalgeschwindigkeiten und radialem Druckgradienten im 
Nachlauf wurde außerdem bereits in den Messungen des Basisdruckes für 
das Modell mit A = 16° in Abbildung 5.6 sichtbar. In Abbildung 7.10 wurden 
deshalb die zentrifugale Volumenkraftdichte pu; / r in einem mitrotierenden 
Koordinatensystem mit dem radialen Druckgradienten öp/ör im Nachlauf der 
Modelle mit Leitflachen miteinander verglichen. 


Für das Modell mit A = 0° zeigt Abbildung 7.10a die radialen Druckgra- 
dienten, die sowohl aufgrund der Verdichtungsstöße und Expansionen im 
Bereich der Leitflächen entstanden als auch aufgrund der Expansion an der 
Modellhinterkante sowie der Rekompression im Nachlauf. Die zentrifugalen 
Volumenkraftdichte war aufgrund der vernachlässigbaren Azimutalgeschwin- 
digkeit der Strömung im Nachlauf des Modells mit A = 0° erkennbar gering. 


Für diein den Abbildungen 7.10b und 7.10c gezeigten Modelle mit angestellten 
Leitflächen waren die zentrifugalen Volumenkraftdichten aufgrund der in der 
Strömung erzeugten Rotationsraten sichtbar höher als für das Modell mit 
A = O°. Für das Modell mit 2 = 32° zeigt Abbildung 7.10c unmittelbar 
hinter den Leitflächen ein Gebiet hoher Zentrifugalkräfte. Den zentrifugalen 
Volumenkräften wirkte ein positiver radialer Druckgradient von annähernd 
gleicher Größe entgegen. Das lokale Gleichgewicht zwischen radialem 
Druckgradienten und zentrifugaler Volumenkraftdichte führte dazu, dass 
die Strömung nach der Ablösung an der Modellhinterkante nur geringfügig 
in Richtung der Zentralachse umgelenkt wurde. Der radiale Druckgradient 
im Bereich 1 < r/R < 1,6, der sich stromab des Modells ergab, resultierte 
in einem im Vergleich zu dem Modell mit 4 = 0° verminderten Druck im 
Rezirkulationsgebiet. 


Für das Modell mit A = 16° waren die in Abbildung 7.10b gezeigten zentri- 
fugalen Volumenkräfte unmittelbar hinter den Leitflachen (x/R > 0 und 1 < 
r/R < 1,625) um eine Größenordnung geringer als für das Modell mit A = 32°. 
Deshalb überwog der radiale Druckgradient unmittelbar hinter der Modellhin- 
terkante die zentrifugale Volumenkraftdichte, sodass die Strömung zunächst in 
Richtung der Zentralachse umgelenkt wurde. Durch die Drehimpulserhaltung 
nahmen die konvektiv in Richtung der Zentralachse transportierten Azimu- 
talgeschwindigkeiten mit abnehmendem Abstand zur Rotationsachse zu, wie 
in Abbildung 7.8a zu erkennen ist. Insgesamt wuchsen die zentrifugalen 
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x/R 
(b) A = 16° 
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Abbildung 7.10: Vergleich des simulierten, zeitlich und azimutal gemittelten, radialen Druckgra- 
dienten und der zentrifugalen Volumenkraftdichte in einem mit dem Fluidelement dV mitrotie- 
renden Koordinatensystem für Modelle mit Leitflächen, eingezeichnet sind die zeitlich gemittelten 


Stromlinien, die durch die Leitflächenvorderkanten (——) hervorgerufenen Verdichtungsstöße 
sowie deren Reflexionen von der Windkanalwand (— —), der Verdichtungsstoß am Ende der 
Lavaldüse (+: ), der Expansionsfächer an der Modellhinterkante (---), die Rekompression 


(— - —) und die Isotache UX = 0. 


Volumenkräfte somit durch den konvektiven Transport der Azimutalgeschwin- 
digkeiten in Richtung der Zentralachse stark an. Bevor die Strömung die 
Zentralachse erreichte, wurde die Strömung deshalb bei x/R ~ 2 parallel 
zur Modellachse ausgerichtet. Der in Strömungsrichtung ansteigende Druck- 
gradient im Bereich der Rekompression war im Gegensatz zu dem Modell 
mit A = 32° stark genug, dass Teile der Strömung zurück zur Basis strömten. 
Auf diese Weise wurde ein Teil des Drehimpulses in das Strömungsgebiet 
unmittelbar hinter dem Modell transportiert, sodass die dort eingebrachten 
zentrifugalen Volumenkräfte in dem in Abbildung 7.10b gezeigten radialen 
Druckgradienten resultierten. 
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7.2.3 Ähnlichkeit zum Wirbelaufplatzen 
in Überschallströmungen 


Die Geschwindigkeitsfelder für das Modell mit A = 32° in den Abbil- 
dungen 7.2c, 7.3c und 7.8b zeigen eine Strömung, die um ein nahezu 
stagnierendes Totwassergebiet mit stromaufgerichteter Axialgeschwindigkeit 
rotiert. Stromab wird dieses Gebiet begrenzt durch die an den Windkanal- 
wänden reflektierten Verdichtungsstöße, die durch die Leitflächenvorderkanten 
hervorgerufen wurden. Settles und Cattafesta (1993) hatten ein ähnliches 
Strömungsfeld, das, wie in Abbildung 7.11 skizziert, durch einen konischen 
Verdichtungsstoß ausgelöst worden war, als Wirbelaufplatzen in Überschall- 
strömungen (Delery et al., 1984) beschrieben. Die Rezirkulationsblase hatte 
sich, ähnlich zu dem Strömungsfeld im Nachlauf des Modells mit A = 32°, 
hauptsächlich stromauf des konischen Verdichtungsstoßes ausgebildet. 


Im Nachlauf des Modells mit A = 32° erstreckte sich die Rezirkulationsblase 
stromauf bis zur Modellbasis. Deshalb ist im Gegensatz zu den experimentellen 
Untersuchungen von Settles und Cattafesta (1993) kein Übergang in Form eines 
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Abbildung 7.11: Skizze zum Wirbelaufplatzen in Überschallströmungen, ausgelöst durch einen 
konischen Verdichtungsstoß (Settles und Cattafesta, 1993). 
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abgelösten Verdichtungsstoßes zwischen dem ungestörten Wirbelschlauch 
und der Rezirkulationsblase erkennbar. Trotz der durchgeführten numerischen 
Simulationen war es somit nicht möglich die Strömungsbedingungen der 
ungestörten Wirbelströmung zu bestimmen. Deshalb konnte auch kein exakter 
Vergleich mit den von Delery et al. (1984) sowie Settles und Cattafesta (1993) 
bestimmten Grenzbedingungen für Wirbelaufplatzen in Überschallströmungen 
durchgeführt werden. 


Abbildung 7.12 zeigt die Grenzbedingungen des Wirbelaufplatzens in Über- 
schallströmungen (Delery et al., 1984, Settles und Cattafesta, 1993) in Abhän- 
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Abbildung 7.12: Vergleich der Strömungsbedingungen im Nachlauf der Modelle mit Leitflächen 
mit den Grenzbedingungen für Wirbelaufplatzen in Unter- und Überschallströmungen (Hall, 1972, 
Delery et al., 1984, Settles und Cattafesta, 1993). 
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gigkeit von der Machzahl M, normal zum auslösenden Verdichtungssto und 
der Wirbelstärke 
U y,max 


Ux ‚ext 


mit der maximalen Umfangsgeschwindigkeit U, max innerhalb des Schlauch- 
wirbels und der Axialgeschwindigkeit U, ext der äußeren Strömung (Délery, 
1994). Zusätzlich ist die Grenzbedingung für Wirbelaufplatzen in inkompres- 
siblen Strömungen nach Hall (1972) angegeben. 


S= 


(7.3) 


Die Machzahl Mn wurde fiir die Modelle mit Leitflächen aus den der Simulation 
entnommenen statischen Drücke vor (pı) und hinter (p2) den reflektierten 
Leitflächenstößen in der äußeren Strömung zu 


(7.4) 


bestimmt (Anderson, 2003). Die nominale Wirbelstärke wurde zu S = tan 4 
in Abhängigkeit vom Leitflächenanstellwinkel A bestimmt und in Abbil- 
dung 7.12 dargestellt. Die maximal erreichte Wirbelstärke im Nachlauf 
entspricht jedoch der, in Abbildung 7.8 eingezeichneten, maximalen relativen 
Azimutalgeschwindigkeit U% ma, und wurde in Abbildung 7.12 durch die 
Fehlerbalken gekennzeichnet. Sowohl die nominale als auch die maximale 
Wirbelstärke des Nachlaufs hinter dem Modell mit A = 32° waren höher als 
die von Settles und Cattafesta (1993) bestimmte notwendige Wirbelstärke, 
um bei einer Stoß-Machzahl Mp = 1 zum Aufplatzen eines Schlauchwirbels 
zu führen. Für den Nachlauf des Modells mit A = 16° konnte aufgrund der 
Abweichungen zwischen der nominalen und maximalen Wirbelstärke keine 
eindeutige Zuordnung zu dem von Settles und Cattafesta (1993) bestimmten 
über- oder unterkritischen Bereich vorgenommen werden. 


7.3 Fazit zu den numerischen Ergebnissen 


Die in Kapitel 5 experimentell beobachteten und beschriebenen Veränderungen 
der Nachlaufstruktur konnten mittels DES-Simulationen nachgebildet werden. 
Die Simulationsergebnisse zeigen, dass die Veränderungen im Nachlauf des 
Modells mit A = 16° hauptsächlich durch die mit der Rotation der Strömung 
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um die Zentralachse eingebrachten Zentrifugalkräfte hervorgerufen wurden, 
die sich in einem mit dem Fluidelement dV mitrotierenden Koordinaten- 
system ergeben. Die Zentrifugalkräfte hinderten die Scherschicht daran bis 
zur zentralen Achse des Strömungsfeldes vorzudringen. Auf diese Weise 
entstanden das gegenläufig rotierende torusförmige Wirbelpaar und der strom- 
abgerichtete Wirbelschlauch entlang der Zentralachse. Für einen Leitflächen- 
anstellwinkel von 32° waren die zentrifugalen Volumenkräfte so hoch, dass 
die Strömung im Bereich der Modellhinterkante nur geringfügig in Richtung 
der Strömungsachse abgelenkt wurde. In der Folge wurde die Rekompression 
im Nachlauf abgeschwächt, sodass nur ein kleiner Anteil der Strömung in 
Richtung der Basis umgelenkt wurde. Erst die an den Windkanalwänden 
reflektierten Verdichtungsstöße der Leitflächen erzeugten im fernen Nachlauf 
einen Druckgradienten, der einen Teil der Strömung in Richtung der Basis 
umlenkte. Die auf diese Weise erzeugte rotierende Strömung um ein Totwas- 
sergebiet weist Ähnlichkeiten mit dem Aufplatzen von Wirbelschläuchen in 
Überschallströmungen auf. Ob es sich dabei um das gleiche Phänomen handelt, 
konnte unter anderem aufgrund der nicht messbaren Strömungsbedingungen 
des ungestörten Schlauchwirbels jedoch nicht abschließend geklärt werden. 
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Im Rahmen der vorliegenden Arbeit wurde der Einfluss des Dralls auf 
den Überschallnachlauf längsangeströmter zylindrischer Körper untersucht. 
Zunächst wurden Windkanalexperimente mit rotierenden Strömungskörpern 
durchgeführt. Obwohl ein Druckabfall an der Basis von bis zu 10 % in Abhän- 
gigkeit von Modelldrehrate und -geometrie gemessen wurde, waren die im 
Nachlauf erzeugten Drallraten nicht hoch genug, um zu einer Änderung der 
Strömungsstruktur zu führen. Höhere Drallraten konnten mit festmontierten 
Strömungskörpern erzeugt werden, die den Drall mittels vor der Modell- 
hinterkante angebrachter, gegenüber der Anströmung angestellter Leitflächen 
in die Strömung einbrachten. Die Modelle hatten unterschiedliche Leitflächen- 
anstellwinkel A. Auf diese Weise konnte der für diese Untersuchung notwendige 
Drall mit akzeptablem technischem Aufwand erzeugt werden. Zusätzlich 
durchgeführte numerische Simulationen wurden anhand der experimentellen 
Ergebnisse validiert. Es wurde gezeigt, dass die verwendete turbulenzauf- 
lösende Grobstruktursimulationsmethode (DES) im Gegensatz zu RANS- 
Methoden geeignet ist rotierende Überschallnachlaufströmungen qualitativ 
abzubilden. 


Für einen Strömungskörper mit einem Leitflächenanstellwinkel A von 0° 
zeigen die Mess- und Simulationsergebnisse im Vergleich zu einem längs- 
angeströmten Zylinder ohne Leitflächen ausschließlich quantitative Verän- 
derungen. Aufgrund der veränderten turbulenten Durchmischung und des 
geringeren statischen Druckes vor der Modellhinterkante verringerten sich 
beispielsweise sowohl der gemessene als auch der simulierte Basisdruck um 
6%. 


Im Gegensatz zu den Modellen, die keine Rotation im Nachlauf erzeugten, 
konnten für die Modelle mit angestellten Leitflächen experimentell und 
numerisch zwei grundsätzlich unterschiedliche, von der Drallrate abhängige 
Nachlaufstrukturen gezeigt werden. Die strukturellen Änderungen konnten 
experimentell sowohl im axialen und radialen Geschwindigkeitsfeld des 
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Nachlaufes als auch der Orientierung des wandnahen Strömungsfeldes und dem 
Druck an der Modellbasis gezeigt werden. Die Änderungen der Strömungs- 
struktur gingen dabei mit einer Verringerung des Basisdruckes von bis zu 80 % 
gegenüber einer nichtrotierenden Nachlaufströmung einher. Die Simulationen 
gaben zusätzlichen Aufschluss über die azimutale Geschwindigkeitsverteilung 
und den Druck im Strömungsfeld. Außerdem konnten die Nachlaufströmungen 
mit Hilfe der numerischen Simulationen auch außerhalb des experimentell 
zugänglichen Messbereiches untersucht werden. 


Für Leitflächenanstellwinkel A von 8° und 16° wurde durch die PIV-Messungen 
in Verbindung mit den numerischen Ergebnissen für A = 16° gezeigt, dass die 
an der Modellhinterkante ablösende Strömung durch die eingebrachten Zentri- 
fugalkräfte parallel zur Modellachse ausgerichtet wurde, bevor die Strömung 
die Achse erreichte. Ein Teil der Strömung konnte, ähnlich zu der klassi- 
schen Nachlaufströmung, den im Bereich der Strömungsumlenkung entste- 
henden Druckgradienten nicht überwinden und wurde in Richtung der Basis 
umgelenkt. Auf diese Weise wurde ein Teil des azimutalen Strömungsimpulses 
konvektiv in Richtung der Basis transportiert. Im Gegensatz zu einer klassi- 
schen Überschallnachlaufströmung wurde die Strömung aufgrund des vorhan- 
denen Dralls zusätzlich radial in Richtung der Modellhinterkante umgelenkt 
und bildete somit einen äußeren Wirbelring aus, der einen gegenläufig rotie- 
renden, inneren Wirbelring umgab. Aufgrund des inneren Wirbelrings änderte 
sich die radiale Richtung der wandnahen Strömung an der Modellbasis im 
Vergleich zur klassischen Nachlaufstruktur. Durch die Überlagerung mit der in 
das Rezirkulationsgebiet eingebrachten azimutalen Geschwindigkeitskompo- 
nente strömte das Fluid in der Nähe der Modellbasis spiralförmig in Richtung 
der Modellachse. Außerdem entstand durch den in das Rezirkulationsgebiet 
transportierten Drehimpuls und die daraus resultierenden Zentrifugalkräfte ein 
radialer Basisdruckgradient. An der zentralen Achse des Strömungsfeldes löste 
die Strömung von der Basis ab und bildete einen zentralen, stromabgerichteten 
Wirbelschlauch aus. 


Für einen Leitflächenanstellwinkel von A = 32° konnten sowohl die experi- 
mentellen als auch die numerischen Ergebnisse eine weitere, bisher nicht 
beschriebene, Veränderung der Nachlaufstruktur nachweisen. Durch den hohen 
Anstellwinkel der Leitflächen waren die Zentrifugalkräfte in der ablösenden 
Strömung unmittelbar hinter der Modellhinterkante so hoch, dass die Strömung 
nur geringfügig in Richtung der Modellachse abgelenkt wurde. In der Folge 
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war die Rekompression aufgrund der Umlenkung der Strömung parallel zur 
Zentralachse so gering, dass die Strömung den entstehenden Druckgradienten 
mühelos überwinden konnte. Im Vergleich zur klassischen Überschallnach- 
laufströmung und der für A = 8° und A = 16° beschriebenen Nachlaufstruktur 
war die unmittelbar hinter der Strömungsablösung erzeugte Rückströmung 
in Richtung der Basis vernachlässigbar. Allerdings beschleunigte die Scher- 
schicht aufgrund viskoser und turbulenter Effekte einen Teil der inneren 
Strömung in axialer Richtung. Aufgrund der Massenerhaltung kam es deshalb 
im Bereich der Strömungsachse zu einer stromaufgerichteten Strömung. Die 
numerischen Simulationen konnten zeigen, dass dieses stromaufgerichtete 
Strömungsfeld im fernen Nachlauf durch die an den Windkanalwänden 
reflektierten Verdichtungsstöße der Leitflächenvorderkanten begrenzt wurde. 
Die reflektierten Verdichtungsstöße erzeugten in der Strömung des fernen 
Nachlaufs einen Druckgradienten in Strömungsrichtung, den ein Teil der 
Strömung nicht überwinden konnte und somit an der Zentralachse in Richtung 
der Basis umgelenkt wurde. 


Weitere experimentelle Untersuchungen für einen Strömungskörper mit einem 
Leitflächenanstellwinkel A von 24°, zeigten in den Visualisierungen des 
wandnahen Strömungsfeldes und den Basisdruckmessungen Ähnlichkeiten zu 
dem Strömungsfeld, das für A = 32° bestimmt wurde. Die Messergebnisse 
für das axiale und radiale Geschwindigkeitsfeld des Nachlaufes hingegen sind 
den Ergebnissen für A = 8° und A = 16° ähnlich. Der zentrale Massenstrom, 
der für die Vermessung des Geschwindigkeitsfeldes mittels Particle-Image- 
Velocimetry-Methode notwendig war, könnte das Messergebnis beeinflusst 
haben. Vermutlich findet die Transition zwischen den beiden in dieser Arbeit 
erstmals experimentell nachgewiesenen Überschallnachlaufstrukturen bei 
einem Leitfächenanstellwinkel A von etwas weniger als 24° statt. 


Die Struktur der Nachlaufströmung hinter dem Modell mit einem Leitflächen- 
anstellwinkel A von 32° weist des Weiteren Ähnlichkeiten mit dem Phänomen 
des Wirbelaufplatzens in Überschallströmungen auf, das bei der Interaktion 
von Schlauchwirbeln mit Verdichtungsstößen entstehen kann. Allerdings 
konnte der ferne Nachlauf, in welchem der durch die Leitflächen erzeugte 
Schlauchwirbel mit den an den Leitflächenvorderkanten entstehenden und 
den Windkanalwänden reflektierten Verdichtungsstößen interagiert, mit dem 
verwendeten Versuchsaufbau experimentell nicht vermessen werden. 
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In nachfolgenden Arbeiten sollte deshalb experimentell untersucht werden, ob 
die für A = 32° beobachteten Veränderungen der Nachlaufstruktur tatsächlich 
aus einer Art Wirbelaufplatzen im fernen Nachlauf resultieren oder ob die 
Ähnlichkeiten der Strömungsstrukturen zufälliger Natur sind. Verbunden mit 
dieser Fragestellung ist auch die Untersuchung, ob sich die Nachlaufströmung 
für A = 32° verändert, wenn es zu keiner Interaktion der zentralen Strömung mit 
den an den Windkanalwänden reflektierten Verdichtungsstößen der Leitflächen 
kommt. 


Außerdem sollten mit einer parametrischen Studie die Grenzbedingungen 
für die beobachteten Strukturänderungen des Nachlaufs genauer untersucht 
werden. Für Leitflächenanstellwinkeln von weniger als 8° und im Bereich von 
etwa 24° könnten zur Bestimmung der Transitionspunkte zwischen den drei 
unterschiedlichen Nachlaufstrukturen beispielsweise numerische Simulationen 
verwendet werden, wie sie in der vorliegenden Arbeit durchgeführt wurden. 
Interessant wäre auch eine Untersuchung der Abhängigkeit der strukturellen 
Veränderungen des Nachlaufs von der Machzahl der Anströmung und einem 
Modellanstellwinkel. 


Aufgrund der veränderten Strömungsstruktur und der damit verbundenen 
veränderten turbulenten Mischungseigenschaften sind des Weiteren Untersu- 
chungen zur Interaktion mit Gasausblasungen und Verbrennungsprozessen im 
Rezirkulationsgebiet interessant. 
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A RANS-Simulationen 


A.1 Numerische Methoden 


Für die RANS-Simulationen wurde ein druckbasierter, gekoppelter Glei- 
chungslöser (ANSYS, 2013) basierend auf einer Finite-Volumen-Methode 
verwendet. Als Turbulenzmodell wurde das k-w-SST-Wirbelviskositätsmodell 
nach Menter (1994) verwendet. 


Die verwendeten Diskretisierungsmethoden sind identisch mit den, für die 
DES-Simulationen in den Kapiteln 6 und 7 eingesetzten, Methoden, die in 
Kapitel 6 beschrieben wurden. Lediglich für die Projektion des Strömungs- 
impulses aus dem Zellzentrum auf die Zellflächen zur Lösung des Riemann- 
Problems wurde ein anderes Diskretisierungsverfahren angewendet. Das für die 
RANS-Simulationen verwendete Diskretisierungsverfahren dritter Ordnung 
nach van Leer (1979) interpoliert aus Stabilitätsgründen zwischen einem 
Upwind-Verfahren zweiter Ordnung (Barth und Jespersen, 1989) und einem 
zentralen Differenzenverfahren zweiter Ordnung (ANSYS, 2013). 


Tabelle A.1: Zeitlich und räumlich gemittelter Basisdruckkoeffizient für verschiedene räumliche 
Diskretisierungsdichten und Modelle, berechnet mittels RANS-Methode. 


m hx Com. 

= Zylinder A=0° A= 16° A = 32° 
1 1 —0,164 -0,166 -0,275 0,340 
2 u -0,166 -0,165 -0,275 -0,341 
3 3 35 0,169 -0,163 0,274 -0,341 
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Abbildung A.1: Abhängigkeit der mittels RANS-Methode berechneten Basisdruckkoeffizienten 
von der räumlichen Auflösung des Rechengitters, eingezeichnet ist auch der theoretisch minimale 
Basisdruckkoeffizient von —0,36 für pp = 0. 


Tabelle A.1 zeigt die über die Basisfläche Ap gemittelten Basisdruckkoeffizi- 
enten Cy fiir die Rechengitter, die auch fiir die DES-Simulationen verwendet 
wurden. Aufgrund der fiir RANS-Simulationen relativ hohen Gitterauflösung 
ist die in Abbildung A.2 gezeigte Abhängigkeit des Basisdruckkoeffizientens 
Cp, von der Anzahl der Gitterzellen äußerst gering. Der Gitterkonvergenzindex 
ist im Vergleich zu den DES-Simulationen um nahezu eine Größenordnung 
geringer und damit im Verhältnis zu den Unterschieden des Basisdruckkoef- 
fizienten zwischen den verschiedenen Modellen so klein, dass dieser in Ab- 
bildung A.1 die gleiche Größenordnung wie die eingezeichneten Symbole 
aufweist. Trotzdem waren die RANS-Simulationen nicht in der Lage, den 
gemessenen und mittels DES-Methode simulierten Unterschied des Basis- 
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druckes in der Größenordnung von 6% zwischen dem längsangeströmten 
Zylinder und dem Modell mit A = 0° wiederzugeben. Außerdem zeigt Ab- 
bildung A.2, dass sich in den RANS-Simulationen für den Nachlauf hinter 
dem Modell mit A = 32° ein deutlich niedrigerer Basisdruck als in den Experi- 
menten ergab. 


A.2 Ergebnisse 


Abbildung A.2 vergleicht die mittels RANS-Methode berechneten Axial- 
geschwindigkeiten im Nachlauf der Modelle mit Leitflächen mit den Ergeb- 
nissen der PIV-Messungen aus Kapitel 5.2. Qualitativ entsprechen die mittels 
RANS-Methode simulierten Geschwindigkeitsfelder für A = 0° und A = 16° 
den Messergebnissen. Die quantitative Abweichung von über 10% der 
mittels RANS-Methode simulierten Länge des Rezirkulationsgebietes für 
das Modell mit A = 0° von den Messergebnissen ist jedoch höher als für 
die DES-Ergebnisse. Außerdem war die RANS-Simulation für A = 32° 
im Gegensatz zur DES-Simulation nicht in der Lage, die experimentell 
beobachtete Strömungsstruktur qualitativ wiederzugeben. 


Abbildung A.3 vergleicht die mittels RANS-Methode simulierte turbulente 
kinetische Energie k mit den PIV-Messergebnissen. Die turbulente kinetische 
Energie der RANS-Ergebnisse wurde hierfür zu k* = k/U2, normiert. Im 
Allgemeinen zeigen die Ergebnisse der RANS-Simulation im Nachlauf eine 
höhere turbulente kinetische Energie als die PIV-Ergebnisse. Dies ist teilweise 
durch die in den Simulationsergebnissen enthaltenen azimutalen sowie 
durch kleinskalige Turbulenz verursachten Geschwindigkeitsfluktuationen 
zu erklären, die experimentell nicht vermessen werden konnten. Außerdem 
überschätzt das verwendete k-w-SST-Wirbelviskositätsmodell (Menter, 1994) 
die turbulente kinetische Energie aufgrund der Druckgradienten im Bereich der 
Rekompression tendenziell (Oliver et al., 2007, Tan und Jin, 2011). Zusätzlich 
kann es in den Bereichen nichtisotroper Turbulenz, wie beispielsweise in der 
Scherschicht, zu systematischen Abweichungen gegenüber den experimen- 
tellen Ergebnissen gekommen sein, da das verwendete Turbulenzmodell auf 
der Annahme isotroper Turbulenz beruht (Menter, 1994). 
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RANS-Simulation, azimutal gemittelt 


PIV-Messung, zeitlich gemittelt, |] 
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U* 
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Abbildung A.2: Vergleich des gemessenen (PIV) und simulierten (RANS), azimutal gemit- 
telten, axialen Geschwindigkeitsfeldes, eingezeichnet sind die zeitlich gemittelten Stromlinien, 


die Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten (——) sowie deren Refle- 
xionen von der Windkanalwand (— —), der Verdichtungsstoß am Ende der Lavaldüse (------- ), der 
Expansionsfächer an der Modellhinterkante (-—--), die Rekompression (— - —) und die Isotache 
Ux = 0. 


A.3 Fazit zu den RANS-Simulationen 


Das verwendete Wirbelviskositätsmodell war nicht in der Lage, die Auswir- 
kungen nichtisotroper Turbulenz korrekt zu beschreiben, sodass der Einfluss 
der Leitflachen auf die Stromungsbedingungen in der Scherschicht nicht 
korrekt wiedergegeben wurde. Deshalb konnte der gemessene und mittels 
DES-Methode simulierte Unterschied des Basisdruckes zwischen dem längs- 
angeströmten Zylinder und dem Modell mit A = 0° in der Größenordnung 
von 6% durch die RANS-Simulationen nicht gezeigt werden. Durch die 
Verwendung eines Reynolds-Stress-Modells, mit welchem auch anisotrope 
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RANS-Simulation, azimutal gemittelt 
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Abbildung A.3: Vergleich der gemessenen (PIV) und simulierten (RANS), azimutal gemittelten, 


turbulenten kinetischen Energie kh yp beziehungsweise k*, eingezeichnet sind die zeitlich 
eI: 


gemittelten Stromlinien, die Verdichtungsstöße, hervorgerufen durch die Leitflächenvorderkanten 
(——) sowie deren Reflexionen von der Windkanalwand (— —), der Verdichtungsstoß am Ende 
der Lavaldüse (+ ), der Expansionsfächer an der Modellhinterkante (— — —), die Rekompression 
(— - —) und die Isotache Ux = 0. 


Turbulenz modelliert werden kann, könnte die Übereinstimmung der Simula- 
tionsergebnisse mit den Messungen mutmaßlich verbessert werden. 


Außerdem war die RANS-Simulation für das Modell mit A = 32° nicht in 
der Lage, die Interaktion des rotierenden Strömungsfeldes im fernen Nachlauf 
mit den Verdichtungsstößen, hervorgerufen durch die Leitflächenvorderkanten 
sowie durch den Übergang zwischen Lavaldüse und Messkammer, korrekt 
wiederzugeben. 
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A RANS-Simulationen 


Die hauptsächlich durch Zentrifugalkräfte erzeugte Strömungsstruktur des 
Nachlaufes hinter dem Modell mit 2 = 16° konnte von den RANS-Simulationen 
hingegen zumindest qualitativ abgebildet werden. 
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B Massenstrom in das nahe 
Nachlaufgebiet verursacht durch 
PIV-Messungen 


B.1 Messung des maximalen Massenstroms 


PB 7% 1,0: 10° Pa 


Tp ~ 295K 
AVB < 0,1 mô 


Ballon Schlauch Modell Vakuumbehälter 


Abbildung B.1: Skizze des Versuchsaufbaus zur Bestimmung des maximalen Massenstroms durch 
das Modell während der PIV-Messungen. 


Zur Bestimmung des maximalen Massenstroms in das Nachlaufgebiet, der 
während der PIV-Messungen aufgrund des Unterdruckes an der Basis 
gegenüber den atmosphärischen Bedingungen außerhalb des Windkanals 
entstand, wurde der in Abbildung B.1 skizzierte Versuchsaufbau verwendet. 
Luft stromte mit einem nahezu konstanten Druck pg ~ 1,0 - 10° Pa aus einem 
Ballon durch einen etwa 3 m langen Schlauch und das Modell in einen Vakuum- 
behälter mit dem Druck pyg < 1,5. 10? Pa. Der Schlauch und das Modell 
wurden in der Konfiguration verwendet, wie sie während der PIV-Messungen 
zum Einsatz gekommen waren. Aus der Volumenänderung des Ballons AVg 
über die Zeit wurde mit dem idealen Gasgesetz (Oswatitsch, 1952) und der 
Temperatur der Luft im Ballon Tg ~ 295 K der mittlere Maximalmassenstrom 
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durch das Modell zu ritmax = (4,65 + 0,10) - 1075 kg/s bestimmt. Da der Druck 
P» an der Basis während der Versuche höher war als pyg und die Strömung 
durch die zur Partikelerzeugung eingesetzten Zigaretten zusätzliche Reibungs- 
verluste erzeugte, war der tatsächlich auftretende Massenstrom mpyy in das 
Nachlaufgebiet während der PIV-Messungen nicht größer als der gemessene 
maximale Massenstrom Mmax- 


B.2 Einfluss des PIV-Massenstroms 
auf das Strömungsfeld 


B.2.1 Abschätzung mittels semi-empirischer Korrelation 


Von Danberg (1990) wurde eine semi-empirische Methode zur Bestimmung 
des Einflusses eines Massenstroms in das Rezirkulationsgebiet sowie dessen 
Temperatur auf den Basisdruck entwickelt (Danberg und Nietubicz, 1992). 
Moore und Hymer (2002) verglichen die Theorie von Danberg (1990) mit 
Messergebnissen für Überschallströmungen von Reid und Hastings (1959), 
Bowman und Clayden (1967) sowie Mathur und Dutton (1996) und konnten 
bis zu einem entdimensionalisierten Massenstrom von etwa 0,02 eine gute 
Übereinstimmung zeigen. Für die PIV-Messungen ergibt sich aus dem 
Verhältnis des Massenstroms in das Rezirkulationsgebiet mpry zu dem Massen- 
strom ApwUs/(nD?), den die Anströmung bezogen auf den Strömungsquer- 
schnitt des Modells hatte, ein entdimensionalisierter Massenstrom von weniger 
als 5,5 - 1075. Der entdimensionalisierte Massenstrom liegt somit innerhalb 
des von Moore und Hymer (2002) gezeigten Gültigkeitsbereiches der von 
Danberg (1990) entwickelten Methode, sodass der Einfluss auf den mittleren 
Basisdruck zu weniger als 0,1 % bestimmt werden konnte. 


B.2.2 DES-Simulationen mit zentralem Massenstrom 


Über einen lokalen Einfluss im Bereich der Einströmung sowie über einen 
Einfluss auf die Struktur des Nachlaufes kann mit der Theorie von Danberg 
(1990) allerdings keine Aussage getroffen werden. Deshalb wurden zusätz- 
liche DES-Simulationen auf dem mittleren Gitter m = 2 mit einem zentralen 
Massenstrom von misim. = 5 - 1075 kg/s durchgeführt. Dieser Massenstrom ist 


180 
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größer als der experimentell ermittelte maximale Massenstrom und damit in 
jedem Fall größer als der während den PIV-Messungen erzeugte Massenstrom. 
Abbildung B.2 vergleicht die simulierten mittleren Axialgeschwindigkeits- 
felder der Modelle mit Leitflächen mit und ohne zentralem Massenstrom in 
das nahe Nachlaufgebiet. 


Abbildung B.2a zeigt, dass der zentrale Massenstrom die Struktur der 
Nachlaufströmung für 4 = 0° nicht in dem Maße veränderte, wie dies für höhere 
Massenströme durch Mathur und Dutton (1996) experimentell gezeigt worden 
war. Im Vergleich zu dem Simulationsergebnis ohne zentralem Massenstrom 
ist das Rezirkulationsgebiet allerdings um 8% kürzer, was auf den um 5 % 
geringeren Basisdruck zurückzuführen ist. Die Simulationsergebnisse mit 
zentralem Massenstrom zeigen gegenüber der Simulation ohne Massenstrom 
in Abbildung B.3a ein erhöhtes Niveau der relativen turbulenten kinetischen 
Energie Koy ur Insbesondere der von Forsythe et al. (2002a) beschriebene 
unphysikalische Ubergang von der modellierten Turbulenz der Grenzschicht 
zur aufgelösten Turbulenz in der Scherschicht fand in den Simulationen mit 
zentralem Massenstrom weiter stromauf statt. Aufgrund dieses numerischen 
Effektes wurde das Fluid des Rezirkulationsgebietes stärker von der äußeren 
Strömung mitgenommen als in der Simulation ohne zentralem Massenstrom, 
sodass der Druck an der Basis abfiel (Sandberg und Fasel, 2006b). Obwohl sich 
die Struktur der Nachlaufströmung in Abbildung B.2a nicht ändert, zeigt die 
Veränderung des Übergangs zwischen modellierter und aufgelöster Turbulenz 
die Sensitivität der numerischen Simulationen bezüglich numerischer Effekte. 


Für die in den Abbildungen B.2b, B.2c, B.3b und B.3c gezeigten Nachlauf- 
strömungen hinter Modellen mit angestellten Leitflächen sind nur lokale 
Veränderungen der Strömungsfelder in der Nähe der Achse zu erkennen. 
Insbesondere das Einströmen des Gases in den Nachlauf bei x,r = O ist in 
den Abbildungen B.2 und B.3 für alle Modelle gut erkennbar. Insgesamt 
zeigt der Vergleich der numerischen Simulationen mit und ohne zentralem 
Massenstrom für die unterschiedlichen Leitflächenanstellwinkel A bis auf 
lokale Unterschiede im Bereich der zentralen Gaseinströmung eine gute 
Übereinstimmung der Nachlaufstrukturen. Eine grundlegende strukturelle 
Beeinflussung der Nachlaufstruktur durch den für die PIV-Methode notwen- 
digen Partikelmassenstrom konnte somit für Leitflächenanstellwinkel 2 von 
0°, 16° und 32°, basierend auf den numerischen Ergebnissen, ausgeschlossen 
werden. 
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B.2 Einfluss des PIV-Massenstroms auf das Strömungsfeld 


DES-Simulation, m = 2, gemittelt, ms; 


-2 
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* 
US ‚Sim. 
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Abbildung B.2: Einfluss des PIV-Massenstroms sim. = 5 - 1075 kg/s auf die simulierten (DES), 
azimutal und zeitlich gemittelten, relativen Axialgeschwindigkeiten im Nachlauf, eingezeichnet 
sind die zeitlich gemittelten Stromlinien, die Verdichtungsstöße hervorgerufen durch die Leitflä- 


chenvorderkanten (——) sowie deren Reflexionen von der Windkanalwand (— —), der Verdich- 
tungsstoß am Ende der Lavaldüse (+ ), der Expansionsfächer an der Modellhinterkante (— --), 
die Rekompression (— - —) und die Isotache UX = 0. 
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B.2 Einfluss des PIV-Massenstroms auf das Strömungsfeld 


DES-Simulation, m = 2, gemittelt, rsim. = 0 


DES-Simulation, m = 2, gemittelt, rsim. = 5 - 10 kg/s 


*” 
ko ut 
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Abbildung B.3: Einfluss des PIV-Massenstroms sim. = 5 - 1075 kg/s auf die simulierte (DES), 
azimutal und zeitlich gemittelte, relative, turbulente kinetische Energie kh u: im Nachlauf, 
eingezeichnet sind die zeitlich gemittelten Stromlinien, die Verdichtungsstöße hervorgerufen durch 
die Leitflächenvorderkanten (——) sowie deren Reflexionen von der Windkanalwand (— —), der 
Verdichtungsstoß am Ende der Lavaldüse (------: ), der Expansionsfächer an der Modellhinterkante 


(=), die Rekompression (— - —) und die Isotache UX = 0. 
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